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Using thermodynamic functions of gases derived from 
band spectra, calculations have been made of theoretical 
explosion pressures in mixtures of hydrogen and oxygen 
containing various inert gases which are compared with 
explosion pressures observed experimentally by different 
investigators. Details of methods of calculation are given. 
In carefully dried mixtures of hydrogen and oxygen with 
excess hydrogen, the observed pressures are too low. This 
effect is eliminated by the addition of small amounts of 
water vapor to the original mixture which, accorditg to an 
hypothesis advanced by Wohl and von Elbe, quenches the 
luminescence from OH radicals formed in the flame front. 
In mixtures containing excess oxygen or nitrogen the ob- 
served pressures are too high. The existence of a lag in the 
excitation of the vibrational energy levels of oxygen and 
nitrogen during the explosion period, as advanced by 
Wohl and Magat, is offered as explanation. The excitation 
lag is greater in oxygen than in nitrogen. Explosions of mix- 
tures of gaseous ozone diluted with oxygen do not show this 
effect, even though the explosion time for approximately 
the same explosion temperature hardly differs from that of a 
hydrogen-oxygen explosion which shows the excitation lag. 


This is explained on the assumption that there is no ap- 
preciable lag in the transfer of vibrational energy between 
like molecules. A similar excitation lag apparently does not 
exist in H,O and H:2. The phenomenon of time-variable 
heat capacities within the duration of an explosion may be 
used to explain the intense audible vibrations observed by 
the authors and others on exploding slow burning lean 
mixtures of hydrogen with air or oxygen within a certain 
concentration range. The possible bearing of this effect on 
the phenomenon of engine knock is mentioned. Explosions 
were made of mixtures of hydrogen and oxygen with helium 
as inert gas, which are free from both anomalous effects 
mentioned above and the energy of dissociation of the 
reaction 2H,O—=20H+H:2 was calculated from the data 
obtained with the aid of accurate thermodynamic functions 
of the gases involved. The average of the results of six 
explosions was found to be 126,000+2000 calories per 2 
moles of H,O. From this value one obtains 


H,O=H +0OH — (114,400+1000) calories 
OH=H+0O —(104,160+1000) calories. 


HE explosion method has been used con- 
siderably during the past 25 years for the 
determination of heat capacities and dissociation 
data without knowing precisely how much re- 
liance was to be placed on the use of measured 
explosion pressures for such purposes. Because of 
the recent exact evaluation of thermodynamic 
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functions of gases from band spectra one may 
calculate accurately the pressures which should 
be developed provided the explosion is adiabatic 
and provided complete dissociation equilibria 
and equilibria among the different quantum 
states are established. Since the experimental 
data on hydrogen-oxygen explosions are secured 
by agreement of the observations of various in- 
dependent investigators this procedure of com- 
paring calculated with observed pressures must 
bring to light any existing anomaly. 
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It has been usual to assume that the observed 
maximum pressure corresponds to a condition of 
equalized temperature throughout the vessel. 
In reality there is a temperature gradient rising 
inward from the wall to the point of ignition at 
the center. In the thermodynamic calculation of 
explosion pressures the effect of this temperature 
gradient must be corrected for. The pressure for 
the temperature gradient system is lower than 
that for the equalized temperature system by an 
amount (0.1 to 1 percent)* which varies with the 
composition of the explosive mixture. 

The writers have carried out the calculations 
for the existing reliable data®»® on ex- 
plosions of mixtures of hydrogen and oxygen 
containing various inert gases employing reliable 
spectroscopic data® on dissociation constants and 
energy content of the gases involved. 


METHOD OF CALCULATION 


The method of calculation for the different 
types of explosions is basically the same as that 
of Wohl ahd co-workers,’ with the addition that 
the theoretical explosion pressure is corrected 
for the effect of the temperature gradient in 
the exploded gas as outlined in a previous 


paper.* 


4 Lewis and von Elbe, J. Chem. Phys. 2, 659 (1934). 

5 Pier, Zeits. f. Elektrochemie 15, 536 (1909); Bjerrum, 
Zeits. f. physik. Chemie 79, 513 (1912); 81, 281 (1912). 
asa and von Elbe, Zeits. f. physik. Chemie B5, 241 
(932) and Magat, Zeits. f. physik. Chemie B19, 536 

8 Also supplementary explosions by the authors con- 
tained in this paper. 

® Free energy and energy content of water vapor. Gordon, 
J. Chem. Phys. 2, 65, 549 (1934). 

_ Free energy and dissociation constants of hydrogen. 
Giauque, J. Am. Chem. Soc. 52, 4816 (1930). See also 
Davis and Johnston, ibid. 56, 1045 (1934), for slight correc- 
tions above 2000°. 

Energy content of hydrogen. Davis and Johnston, 
J. Am. Chem. Soc. 56, 1045 (1934). 

Free energy, dissociation constants and energy content 
of oxygen. Johnston and Walker, J. Am. Chem. Soc. 55, 
172, 187 (1933). Lewis and von Elbe’s (ibid. 55, 507, 511, 
1933) values of the dissociation constants and energy con- 
tent include the contribution of the oxygen !A level. The 
free energy of Johnston and Walker should be corrected for 
the level. 

_ Constant of equilibrium of NO=1/2 N.+1/2 Os. 
Giauque and Clayton, J. Am. Chem. Soc. 55, 4875 (1933). 

Energy content of nitrogen. Johnston and Davis, J. Am. 
Chem. Soc. 56, 271 (1934). 
Free energy of OH. Johnston and Dawson, J. Am. Chem. 


- Soc. 55, 2744 (1933). 


Heat of formation of H,O from He and Oy. Rossini, 
Bur. Standards J. Research 6, 1 (1931). 
10 See reference 7 for other references. 


In the calculations which follow, the unknown 
quantity sought is Pe(calc.), the pressure for 
the unequalized temperature state actually exist- 
ing at the end of an explosion, which is to be 
compared with Pe(obs.), the observed explosion 
pressure. This is the most direct way of com- 
paring experiment with theory. 

Since H.O dissociates in two different ways, 
namely, 

I. 114,222 


II. 20H+H2—2H,0+ 126,000" 


and the dissociations are affected differently by 
hydrogen and oxygen, the method of calculation 
differs slightly with the type of gas mixture, 
that is, containing 


Case 1. Large excess of hydrogen. 

Case 2. Small excess of hydrogen or stoichiometric 
mixture. 

Case 3. Excess oxygen. 


Case 1. Explosions with large excess of hydrogen 
and with or without inert gases 


The inert gases may be nitrogen, argon and 
water vapor. At explosion temperatures nitrogen 
will neither dissociate nor react with hydrogen 
appreciably, but the dissociation of hydrogen 
must be considered. For the combustion of 1 
mole of hydrogen at room temperature, 


57,501 
= (1+mx.0)|CvD 1.0) 
Té-Ti 


Te Te 
CoD a» +MN:,A CuD a) (1) 
Ti Ti 


where 57,501 is the heat of combustion of 1 mole 
of ordinary hydrogen gas, at constant volume 
and room temperature. This value has been 
determined for 298.1°K, but may be assumed 
constant for a small temperature range. 7@ is 
the temperature which would exist at the end 
of an explosion if the temperature were equalized 
throughout the vessel without loss of heat to 
the wall. 77 is the temperature of the gas mixture 
before explosion. mu:0, mu:, etc., are the moles 
of the indicated gases in the burned mixture per 
mole of HO formed disregarding all dissocia- 
tions, that is, per 1/2 mole of Oz: originally 


11 See below. 
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present. CuD«i.0)|7,, etc., are the heat capacities 
between 7i and Té of the gases indicated in- 
cluding the contribution of dissociation to the 
heat effect. 

For large excess of hydrogen, dissociation J 
may be neglected. Thus 


Té 
+Du (2a) 


CvD 
Ti 


Cuca») +Du,, 


Tt 


a») (2b) 


Dy= re) /(T@— (3a) 


and 


(3b) 


where (AE°re)11 and (AE°rz)n: are the energies of 
dissociation of H:O and Hae, respectively, and 
x1 and xn: are the corresponding degrees of 
dissociation. 

Pe, the explosion pressure for the equalized 
temperature state is calculated from the gas law, 


(4) 


where APe is the excess of pressure due to in- 
crease in the number of molecules caused by 
dissociation 

APe=0.5(pout+pu), (5) 


where Pou and py are the partial pressures of 
OH and H existing in dissociation equilibrium. 

The dissociation equilibria are determined by 
the constants 


Kyu= 
Kou:= pu?/ pu, (7) 
pu.0, pu2, pou, and px are given by 
pu:0= (Té/T1) + pin:o ](1—x11), (8) 
pus= 
+ { 
+ pin.o}0.5x11 ], 
pou= (Té/Ti)[pin:/(1+mu:) + pin.o 
+ { Jxn:, (11) 


(9) 
(10) 


where pin: and pin:o are the partial pressures of 
these gases in the unburned mixture at 77. 

The value of Pe is obtained by the method of 
trial and error. A temperature 7@ is estimated 
and with the aid of the known values of heat 
capacities and dissociation constants,® x11 and xs 
and consequently Di: and Dn: are evaluated. 
If Eq. (1) is not fulfilled a new value of 7@ is 
chosen and the calculation repeated. 

The ratio Pe/Pe is determined separately from 
the relation‘ 


1 
Pe/Pe=Te(+a)/ Te(it+x)dn, (12) 


(1+) is the ratio of the number of molecules 
in the burned mixture in the dissociation equi- 
librium state and the undissociated state for the 
temperature 72. Therefore 


¥= APe/(Pe—APe) 


(13) 
The corresponding symbols for the unequalized 
temperature state are within the integral of Eq. 
(12) which is to be taken over all elementary 
fractions dn of the burned gas. m varies from 0 
to 1. One first obtains a curve of 2m|CvD|7*x Te 
as a function of m as outlined previously.‘ 
The area under this curve equals =m|CvD|}? 
x Té. This enables one to obtain Te and (1+) 
as functions of m and also Té and (1+2). For 
calculation of the ratio Pe/Pe it is preferable to 
use the value of T@ determined from the above 
area rather than by Eqs. (1) to (11), in order 
that slight inaccuracies in the Dm|CvD|7*xTe 
versus n curve shall not enter into the Pe/Pe 
ratio. In any case both values of Té are very 
close to each other. 


Case 2. Explosions with small excess of hydrogen 
or of stoichiometric mixture and with or 
without inert gases 


In this case one must consider both dissocia- 
tions I and II, and in the presence of nitrogen 
the formation of NO, namely, 2NO=N.2+0Oz2 
+42,800. It is advantageous to combine all 
dissociation terms into a single term ZmD, so 
that Eq. (1) now reads: 
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57,501 Té Té 
Té-—Ti Ti Ti 


Te 
+2imD. (14) 
Ti 


Let Mon, My, Mo and myo represent the 
number of moles of the indicated gases at 
equilibrium per mole of H.O formed for complete 
reaction, that is per 1/2 mole of O» originally 
present. 

To bring the gases to equilibrium first allow a 
sufficient amount of H,O to dissociate into He 
and O, and then let Hz and O2 combine accord- 
ing to 


until equilibrium is reached. Then 
(0.5/(Té— T1)) rz) 1(mout+2mo; 
+mo+myo) + (AE* re) + 72) 
(15) 


The following equilibria are to be considered 
in addition to Eqs. (6) and (7) 


Kyr= (16) 
Kvo:= po?/po:; (17) 
Kpxo= pno?/ (18) 


Mon, My, etc., are expressed by the following 
equations 


Mon, x, 0, No= (Pou, x, 0, (19) 


where pio: is the partial pressure of Oz in the 
unburned mixture at the temperature 77. 

Values of pu:, pu.o and pno are affected by 
the multiple dissociations and are given by 


pu:= (pin:—2pio:) Té/Ti+0.5pou 
(20) 
pu.0= (2pio:+ pin.o) Té/Ti— pou 
(21) 
pni= pin:Té/Ti—0.5pno. (22) 


Pe is obtained from Eq. (4), except that APe 
is now 


APe=0.5(pou+Pu) (23) 
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Case 3. Explosions with excess oxygen and with 
or without inert gases 


In this case the same dissociation equilibria as 
in Case 2 are considered. Eq. (19) changes to 


Mou, 0, No= (Pou, N, 0. no/ Pica») Tt/Té. (24) 


The term mo.| Cvo:| is substituted for the He 
term in Eq. (14). In Eq. (15) (AE°7ze); is multi- 
plied by (0.5mon+mu,+0.5my). 

The partial pressures of the various gases to 
be used in the equilibrium equations now read : 


pu.0= (pin: + pin.o) Té/Ti—0.5 pou 
—pu:—0.5pu, (25) 
por= (pio: —0.5pin:) Té/Ti—0.25 pon 
(26) 
px:= pin.Té/Ti—0.5pno. (27) 


Tables I to IV contain the results of the calcu- 
lations. 

For most groups of experiments the original 
values of Ji and Pi have been adjusted slightly 
to a common pressure of 760 mm Hg and a 
common temperature, in most cases 291°K. 
This necessitates a slight correction of the 
observed explosion pressure, which for small 
changes of Pi and Ji is accomplished by the 
relation 


APe=(760—Pi)Pe/Pi—[(291—Ti)/Ti] 
(28) 


The agreement between the calculated and 
observed explosion pressures of the first group 
in Table I is complete, the average deviation for 
15 ‘‘argon” experiments being +0.07 percent, 
that is, zero. The small systematic deviation 
(—0.61 percent) in the excess hydrogen explo- 
sions is accounted for by heat losses due especially 
to heat conducted through the ignition rod which 
is in contact with hot gas during the explosion. 
Such losses have been found experimentally by 
Wohl and Magat for these mixtures of high 
heat conductivity, but not for excess oxygen 
mixtures of low conductivity. 

The value used for the energy of the dissocia- 
tion reaction 20H +H2—2H:O is 126,000 calories, 
which was obtained from explosions (described 
later) of mixtures of hydrogen and oxygen con- 
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TABLE 1. Hydrogen-oxygen explosions with excess hydrogen and — and water vapor. Comparison of Peeatc. and Peors. 
Pi=760 mm Hg. 


Ti(H.0) 


Té Peécaic. Pecalc. Péobs. differ- 


mmHg mmHg mmHg 


nN 
w 


=) 


5926 —0.32 
(6406)* (+1.25) 


+0.07 


wn =e 


Average 


Average 


* Experimentally unreliable. 


TABLE II. Hydrogen-oxygen explosions (dry) with excess ee Comparison of Pecate, and Pe,),... Pi=760 mm Hg, 
Ti=291°K. 


= 


Explo- 
sion Té Té 


Author No. Cy 


Ti(H2) 


5 
Peobs. 


Pecale. 
mm Hg difference 


mm Hg 


Pecale. 
mm Hg 


Bjerrum 
Wohl von Elbe 


do. 
Wohl and Magat 
do. 
do. 


— 1.68 
—2.03 
—1.75 
—1.77 
—1.90 
—1.92 
—1.97 
—1.93 
—2.87 
—2.18 


Average —2.00 


taining a slight excess of oxygen, and helium and 
water vapor as inert gases. The results tabulated 
in Table I are not materially affected by an un- 
certainty in this energy value because, due to 
excess of hydrogen, the water dissociation is 
small. 

In excess hydrogen explosions Wohl and von 
Elbe® found that the absence of water vapor in 
the original mixture caused a marked lowering 
of the explosion pressure. This is confirmed by 
reconsideration of early experiments of Bjerrum® 
and later by the work of Wohl and Magat’ and 
Lewis and von Elbe.” The results are given in 


® Reference 4, and Table IV, this paper. 


Table II. The average deviation of 10 explosions 
is — 2.00 percent. 

Quite the opposite effect is observed in the 
type of explosion listed in Table III. Here the 
experimental pressures are too high, more pro- 
nounced in the excess oxygen (average deviation 
+1.22 percent for 6 explosions) than in the 
excess nitrogen explosions (average deviation 
+0.79 percent for 8 explosions). Since the results 
of different authors agree within the limits of 
error the effect appears to be well substantiated. 
The cause of the effect was attributed by Wohl 
and Magat to a lag in the excitation of the 
vibrational energy levels in oxygen and nitrogen 


Explo- 
sion Ti Té 
Author No. mA °K Cy II He ence 
24) H Pier 138 0.169 0.297 291 2041 0.043 0.163 5041 5009 5015 +0.12 
, do. 135 .166 .240 291 2047 .053 .189 5055 5023 5045 +0.44 
: do. 150 .140 -987 291 2153 048 .164 5270 5236 5274 +0.73 
H ¢ do. 149 .123 .771 291 2381 .180 385 5769 5731 5735 +0.07 
2 F do. 147 .126 .776 291 2392 .187 .395 5789 5751 5760 +0.16 
Iti i do. 148 .119 685 291 2440 -262 511 5897 5858 5853 —0.09 
do. 146 -115 .575 291 ‘2600 .926 6232 6188 6196 +0.13 
z do. 145 .114 .571 291 2613 .571 976 6247 6203 6210 +0.11 
5 do. 144 -113 546 291 2627 -604 1.03 6285 6240 6240 0.00 
t ; do. 155 1.033 401 326 2048 -050 0.120 4468 4449 4481 +0.72 
Oo 7 do. 152 0.877 .677 326 2259 101 .247 4852 4828 4823 —0.10 
do. 151 .830 .627 326 2325 .150 .337 4979 4953 4944 —0.18 
: 153 835.523. 3.095 2348 180 5032, «4994. —0.75 
on 
von Elbe 51 284 4.837 291 2469 441 892 5984 5945 
: do. 52 .103 .379 3.491 291 2672 1.060 1.730 6372 6327 
25) 
: Pier 166 0.101 4.057 0.015 0.045 5030 5017 5012 —0.10 
do. 167 3.398 -036 -118 5412 5398 5384 —0.26 
: do. 168 -086 2.683 .160 .300 5872 5852 5821 —0.53 
: Wohl & Magat 31 .121 3.392 .034 .110 5395 5381 5340 —0.76 
26) do. 37 .009 3.461 .050 -123 5458 5443 5420 —0.42 
} do. 28 oiae 3.007 .080 .177 5629 5613 5582 —0.55 
do. 32 oe 2.924 .100 .200 5695 5678 5625 —0.93 
27) do. 36 .029 3.019 .108 .212 5716 5699 5680 —0.31 
do. 29 .108 2.610 173 .310 5901 5881 5842 —0.66 
do. 38 .008 2.723 .183 .330 5939 5919 5845 —1.26 
cu- do. 30 025 2.657 .195 .350 5951 5940 5885 —0.93 
—0.61 
nal 
tly 
K : 
the 
the 110 3.260 2387 8.55 5.60 0.075 0.161 5595 5579 5485 
47 3.876 2202 8.38 5.53 .024 .069 5525 5512 5100 
49 3.924 2188 8.36 5.52 .024 .066 5194 5181 5090 
do. 43 3.678 2260 8.42 5.56 .037 .094 5343 5329 5235 
do. 44 3.650 2268 8.43 5.56 .040 .094 5357 5343 5242 
do. 45 2.726 2564 8.72 5.68 .180 .330 5941 5921 5808 
46 2.719 2566 8.72 5.68 "180 '330 5946 5926 5810 
28) 34 3.779 2231 8.40 5.54 .037 .081 5282 5268 5166 
27 3.500 2301 8.46 5.57 .048 242 5424 5410 5258 
33 3.324 2368 8.54 5.60 .069 .149 5556 5541 5420 
ind 
Up 
for 
nt, 
ion 
lo- 
ly 
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on. 
by 
igh 
ren 
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TABLE III. Hydrogen-oxygen explosions with excess oxygen and 


G. VON ELBE 


water vapor, and excess hydrogen and nitrogen and water vapor. 


Comparison of Pecaic, and Pes. Pi=760 mm Hg, Ti=291°K. 


Explo- 
sion Té Té Té Pecalc. Peale. fobs. % 
Author No. °K CPT =mD mmHg mmHg mmHg _ difference 
Wohl & Magat 71 2.817 0.090 0.039 2251 8.41 6.50 1.61 5246 5238 5314 +1.45 
Lewis & von Elbe 2.1 2.820 091 .044 2247 8.41 6.50 1.60 5237 5229 5277 +0.91 
Wohl & Magat 72 2.900 095 041 2221 8.39 6.49 1.52 5187 5180 5264 +1.62 
Lewis & von Elbe 2.2 2.897 .094 .045 2221 8.39 6.49 1.52 5187 5180 5235 +1.06 
Wohl & Magat 70 3.044 .096 042 2180 8.35 6.47 1.37 5106 5099 5140 +0.81 
Lewis & von Elbe y © 3.045 097 .046 2176 8.35 6.47 1.36 5097 5090 5165 +1.47 
Average +1.22 
Explo- % 
sion Ti Té Té Té Té Pecalc. Pecale. Pobs. differ- 
Author No MH, mn, °K °K Dy Dit, mmHg mmHg mmHg ence 
Pier 158 0.141 0.476 5.806 291 1589 7.62 5.26 0 0 5.70 3889 3885 3892 +0.18 
do 159 13 546 4.799 291 1744 7.83 5.33 0 0 5.79 4229 4224 4282 +1.38 
do. 160 105 517 3.718 291 1993 8.14 5.44 0.0185 0.068 90 4761 4755 4837 +1.73 
do. 161 096 443 3.186 291 2163 8.33 5.52 -061 168 5.97 5111 §102 5175 +1.43 
do. 162 .086 366 2.519 291 2402 8.56 5.62 -232 496 6.05 5584 5572 5614 +0.75 
do 163 .077 434 1.894 291 2595 8.76 5.70 471 806 6.11 5964 5949 5940 —0.15 
do. 157 1.003 952 3.966 326 1653 7.74 5.30 5.76 3593 3591 3610 +0.53 
do. 156 0.636 356 1.984 326 2347 8.55 5.61 0.18 0.40 6.05 4867 4860 4869 +0.19 


Average +0.79 


TABLE IV. Hydrogen-oxygen explosions (dry)— Miscellaneous mixtures. Comparison of Pecate, and Peobs, 


Explo- 
i Ti Pi 7s v8 


Peet, Fle, 
Té Té fcale. ‘cale. 4 obs.  differ- 


sion 
Author No. my, mo, my, my °K mmHg °K PH: mmHg mmHg ence 


Lewis & 
vonElbe 1.5 2.587 — 1.868 0.024 294.6 717.5 2013 
do. 1.868 .024 294.3 717.5 2735 8.91 
do. 12  — 0.326 3.127 040 295.8 717.0 2187 8.36 
do. 26 — 3452 0.048 — 292.7 757.0 2097 8 


5.45 5.92 4493 4484 84409 —1.67 


5.76 6.18 5799 55780 5806 +-0.45 
6.48 — 1.12 5.98 4782 4773 4762S —0.23 
6.44 1.06 5.95 4894 4888 4924 +0.74 


molecules. That such a lag exists has been shown 
by the work of a number of investigators.” 
Explosion experiments therefore would seem to 
offer a rather novel approach to this field of 
study. 

While the explosion times are too short to 
allow complete excitation of the vibrational 
states of oxygen and nitrogen they seem to be 
sufficient to do so for water vapor and hydrogen 
(Table I). 

The lag in the excitation of vibrational energy 
is not observed in explosions of ozone diluted 
with oxygen previously reported by the authors." 
That the absence of the effect in such explosions 
cannot be explained on the basis of longer explo- 

18 Kneser, Ann. d. Physik 11, 761, 777 (1931); 16, 377 
(1933); Phys. Zeits. 32, 179 (1931); Zeits. f. Physik 77, 
649 (1932); Nature 129, 797 (1932); J. Acous. Soc. Am. 5, 
122 (1933). Zener, Phys. Rev. 38, 27 (1931); wong’! Na- 
ture 129, 200 (1932); Eucken, Miicke and Becker, Natur- 


wiss. 20, 85 (1932); Sharratt and Griffiths, Proc. Roy. Soc. 


A147, 292 (1934). 
4 Lewis and von Elbe, J. Am. Chem. Soc. 55, 511 (1933). 


sion times is shown by comparison in Table V 
of explosion 2.3 (Table III) with a typical ozone 
explosion at nearly the same explosion tempera- 
ture 7@. The last column gives the time elapsed 
between ignition and establishment of maximum 
pressure. It is conceivable that highly excited 
oxygen molecules liberated in the decomposition 
of ozone may readily transfer vibrational energy 
to molecules of their own type (oxygen mole- 
cules), whereas the transfer of such energy from 
freshly formed water molecules to oxygen mole- 
cules is more difficult. If the above argument is 
correct, barring unforeseen complications in ex- 
perimental manipulation, one might also expect 
disappearance of the excitation lag effect in 
nitrogen by exploding mixtures of N2xO and H:2 
diluted with Ne. Such explosions have not yet 
been reported. 

An explanation of the water-vapor effect was 
proposed by Wohl and von Elbe, who suggested 
that OH radicals formed in the flame front 
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TABLE V. Comparison of explosion times of ozone and 
hydrogen-oxygen explosions. 


Explosion Té Time 
Explosion No. (°K) (sec.) 
Ozone 3.3 2122 0.0146 
H.—O» 2.3 2176 .0128 


during the combustion of hydrogen and oxygen 
radiate energy by luminescence; that the radia- 
tion may be quenched by small amounts of 
water vapor initially present in the explosive 
mixture but not by the ‘newly-formed highly 
energized water molecules nor by hydrogen.” 
The question concerning the efficiency of quench- 
ing the luminescence by other gases, such as Oz 
and Ne, cannot be answered from the present 
data. The results of the dry explosions of Lewis 
and von Elbe in Table IV merely suggest a 
superposition of the excitation lag effect on the 
water-vapor effect. The latter predominates in 
explosion 1.5 because the nitrogen content repre- 
sents only’ about 3/7 of the total inert gas. 
The former predominates in explosion 1.3 because 
nitrogen represents practically all the inert gas; 
in addition, the explosion time. is very short, 
thus contributing to the excitation lag. The 
difference between explosions 2.6 and 1.2 is due 
to the fact that the excitation lag in Oz is greater 
than in Ne, the former more effectively masking 
the water-vapor effect. 

The experimental results of explosions 1.5, 1.3 
and 1.2 (corresponding to 59.96 percent, 31.99 
percent, 19.99 percent He in air, respectively) 
have been verified by duplicate explosions in a 
larger vessel.‘ 


DETERMINATION OF THE ENERGY OF DISSOCIA-~ 


TION 20H+H2—2H:0 


If any thermochemical data are to be derived 
from explosions of this kind they must obviously 
be free from the anomalies observed in explosions 
in Tables II, III and IV. The first group of 
explosions in Table I meets these requirements. 
If the slight excess of hydrogen is replaced by a 
slight excess of oxygen the dissociation of H,O 
into OH and Hz becomes of such magnitude as 


4 Concerning criticism of Garner and Tawada’s (Trans. 
non» Apr 26, 36 (1930)) objection to this explanation 
an 


see Wo Magat, reference 7. 


to render such explosions suitable for the de- 
termination of the energy of dissociation. 

This value was determined by Bonhoeffer and 
Reichardt” with a large limit of uncertainty, 
namely, 128,000+5000 calories per 2 moles of 
H,O, by measuring the intensity of light absorp- 
tion in heated water vapor. Suitable explosions 
were carried out by Wohl and Magat, but the 
results of their three published experiments 
differ markedly. Nevertheless they estimate the 
value to be 124,000 calories. 

We have made six explosions using helium 
instead of argon as the inert gas. Previous com- 
parative experiments‘ in the same vessel employ- 
ing identical explosive mixtures of hydrogen and 
oxygen and either helium or argon or various 
mixtures of ‘the two inert gases gave identical 
maximum pressures. Considering the difference 
in heat conductivity of the mixtures, this agree- 
ment indicates the absence of any noticeable heat 
loss through the ignition rod. (The latter was vi 
very small diameter and was made of highly non- 
conductive material.) 

The method of calculating the energy of 
dissociation is that given in Case 3 above, except 
that instead of solving for Pe one solves for 
which determines the value of K 11. 
The results are given in Table VI. 

From this new value the following heats of 
reaction are obtained. The limits of error lie 
practically entirely in the above determination, 
since the values of the Oz dissociation (117,300 
cal.), of the Hz dissociation (102,800 cal.), and 
of the heat of combustion of Hz to HO (57,111 
cal.) are known accurately. 


H,0 =H+OH — (114,400+1000) calories 
OH =H+0-—(104,160+1000) calories. 


Thus it is evident that the separation of the 
first H from H,O requires more energy than the 
separation of the second H. 


GaAs VIBRATIONS IN EXPLOSIONS WITH EXCESS 
OXYGEN AND NITROGEN 


An interesting effect observed in explosions 
containing excess O2 or Ne but never excess He 
or A, which we believe is connected with the 


16 Bonhoeffer and Reichardt, Zeits. f. physik. Chemie 
A139, 75 (1928). 
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TABLE VI. Energy of dissociation 20H +H:=2H,0. 


Expt. Pi 
No. 


Ti Peobs. Pe Té AE% 


mm Hg mm Hg 


Average 


excitation lag, is that when slow-burning mix- 
tures within certain composition limits are ex- 
ploded, intense audible vibrations, which are 
clearly recorded on the pressure diagram, are 
set up.” These vibrations commence long before 
maximum pressure!’ and continue for some time 
thereafter. The limiting percentage of hydrogen 
at which the vibrations disappear was found to 
be about the same in a 12-inch as in an 18-inch 
vessel (the amplitude of the vibrations is some- 
what larger in the larger sphere) and the same 
whether Nez or Os was used as the inert gas. 
The vibrations are more intense in O, than in 
explosions. 

The vibrations have nothing to do with the 
violence of the explosion, because they do not 
appear in richer mixtures which burn more 
rapidly and with greater intensity. The effect is 
shown in the pressure records in Fig. 1. 

Moreover the effect is not due to a resonance 
phenomenon of the vessel because its appearance 
is independent of the size and shape of the vessel. 
Thus Campbell, Littler and Whitworth’® ob- 
served similar vibrations in hydrogen explosions 
with excess oxygen in a cylinder in the same 
composition range as in the present experiments 
in a spherical vessel. Like effects were found in 
explosions of CH, or CH, with excess oxygen. 
Morgan” reported similar vibrations with certain 
mixtures of coal gas or CH, and air exploded in 
a cylinder. Maxwell and Wheeler®™ obtained 

‘7 Also unpublished observations by G. von Elbe and 
Wohl and Magat, reference 7. 

18 This is in agreement with the finding of Wawrzinick 
(Automobiltechnische Zeits., Feb. and March, 1933; 
N. A. C. A. Tech. Mem. No. 711) that the sound started 
before the flame reached the end of the combustion cham- 
ber, the sound being a manifestation of the pressure 


fluctuations inside the bomb, since both had the same 


frequency. 
one Cc ampbell, Littler and Whitworth, J. Chem. Soc. 339, 


20 Morga n, Phil. Mag. 53, 1161 (1927). 
_ 2 Snnaileed Wheeler, Ind. Eng. Chem. 20, 1041 (1928). 


audible vibrations with pentane and air and toa 
smaller degree benzene and air in a cylinder, 
a cube, and a sphere. In all vessels the maximum 
effect appeared with a mixture containing 3.5 
percent pentane. 

Nor are the vibrations due to a natural fre- 
quency of the diaphragm, because the same 
effect is found with different types of diaphragms, 
corrugated and plain; the frequencies of the 
vibrations vary decidedly with the composition 
of the explosive mixture, and do not correspond 
to the natural frequency of the diaphragm. 

It is worthy of note that the vibrations are 
peculiarly connected with the presence of oxygen 
and nitrogen in the exploded gas, both of which 
show the excitation lag, and that the vibrations 
are more intense in the presence of oxygen, which 
has the larger excitation-lag effect. 

The origin of the vibrations may be proposed 
as follows: 

When an elementary fraction of gas burns, the 
Oz and Nz molecules do not take up their normal 
heat capacities immediately. The rate at which 
they acquire normal heat capacities depends 
upon time and principally temperature (probably 
exponential with temperature). If the combustion 
temperature approaches some higher critical 
value, the amount of gas which accumulates in 
the subnormal heat capacity state is small; 
that is, the normal heat capacity is established 
very soon after the passage of the flame front 
through an elementary fraction. The reverse is 
true if the combustion temperature is low, be- 
cause the rate of attainment of normal heat 
capacity is slow. Under the latter conditions the 
burned gas is overexpanded. If this gas is now 
subjected to rather sudden adiabatic compression 
its temperature will rise and correspondingly its 
rate of attainment of normal heat capacity, as a 
consequence of which the mass of gas suffers 
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still further contraction in volume. In a spherical 
vessel this is to be expected (and is found) at 
the stage in the process where the rise in pressure 
is rapid, that is, where a relatively large amount 
of gas is burned in a short time. The mass of 
burned gas in the center undergoes rapid shrink- 
age in volume, which gives rise to a rapid mass 
movement of the gases inward. This may suffice 
to set up a pressure wave which, because of little 
damping, continues for some time. 

In Fig. 1 the temperatures in the burned gases 
for any corresponding position of. the flame front 
are about 450°C higher in the 20 percent He in 
air mixture than in the 15 percent He in air 
mixture. The temperatures of the gases in any 
layer just after the flame has passed it ranges 
from about 1730°C in the central layers to about 
1900°C in the layers near the wall for the 20 
percent H2 mixture; correspondingly, these tem- 
peratures are 1280°C to 1450°C for the 15 per- 
cent Hz mixture. The rate of attainment of 
normal heat capacity in the O2 and Ne molecules 
is much faster in the former explosion, and one 
may conceive that in this case subnormal heat 
capacity exists only in the layers of burned gas 
immediately behind the flame front. Thus the 
excitation lag which is still noticeable at the 
maximum pressures of such mixtures (compare 
Tables III and IV) is confined to the outer or 
last burned layers of gas while the remainder 
of the gas is practically fully excited. These 


conditions are not favorable for initiation of the 
gas vibrations described above. The small excita- 
tion lag effect still observable in the high- 
temperature explosions (compare Ne explosions, 
Table III) is probably due to the very short 
combustion time. 


RELATIONSHIP TO THE PROBLEM OF 
ENGINE KNOCK 


That audible vibrations are concomitant with 
knocking in a gas engine of the Otto type, 
varying in intensity depending upon the severity 
of the knock, has been known for some time.” 
In the light of the above hypothesis on the origin 
of the gas vibrations described, the phenomenon 
of time-variable heat capacities may be given 
serious consideration in connection with studies 
of engine knock. We believe that it plays a role 
in many phenomena accompanying the knock, 
such as loss in power, increased radiation and 
gas vibrations. It would seem that the function 
of antiknock compounds is to guide the course of 
chemical reaction in the yet unburned charge to- 
ward preventing sudden inflammation of the latter 
and thus to reduce to a minimum both the setting 
up of a pressure impulse and the excitation-lag 
effect in the nitrogen gas which dilutes the 
charge. 


22 See for example Withrow and Rassweiler, Ind. Eng. 
Chem. 24, 528 (1932); Automobile Eng. August, 1934. 
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Solutions containing heavy ions in water were examined 
with Cu Ka radiation, care being taken to avoid effects of 
general radiation. The diffraction: patterns of saturated 
solutions of Th(NO3), and UO2(NOs)2 consist mainly in a 
pronounced ring, shifting to smaller angles and getting 
fainter with increasing dilution. This indicates a more or 
less regular arrangement of the ions in the liquid, for which 
the name ‘‘superarrangement”’ will be used to distinguish 
it from the arrangement with smaller periodicity of the 
water molecules. Both arrangements are of course of the 
“Jiquid-type.” The above results show that the “super- 
arrangement” expands and gets more diffuse with increas- 
ing dilution, as might be expected. Another kind of pattern 
is exhibited by solutions of AgNO;, Pb(NO;)2, Ba(NOz)2 
and lead acetate. Here the intensity does not show pro- 


(Received November 2, 1934) 


nounced maxima and minima. More especially it does 
not decrease markedly when the diffraction angle tends to 
zero. An explanation is offered, assuming the existence of a 
large fraction of undissociated molecules, with a gaseous 
distribution. A third class is formed by solutions of Lil, 
RbBr, LiBr and other alkali-halides. These yield a broad 
ring the position of which is independent of concentration. 
It is shown that these facts may be accounted for by 
ascribing this ring mainly to the interference between the 
scattering from a heavy ion and that from the water mole- 
cules around it. To clear up some of the above points a few 
experiments are made with solutions of analogous light 
ions. Finally a general theory on diffraction by solutions is 
formulated and applied to some of the experimental cases. 


§1. 


HERE are reasons to expect that the distri- 

bution of ions in solution is less at random 

than in a gas of corresponding density (or in a 

solution of neutral molecules). I have been occu- 

pied for some time in seeing if some information 

may be gained on this subject by the method of 
x-ray diffraction.': 

From a theoretical point of view it would per- 
haps be preferable to investigate dilute solutions. 
The Debye-Hiickel theory may then be expected 
to be at least a close approximation to the facts. 
Now this theory gives a more or less definite 
picture of the arrangement of the ions. The influ- 
ence of this arrangement on the diffraction pat- 
tern may be easily calculated.! The experimental 
verification would, however, require very accu- 
rate (preferably absolute) measurements as the 


relatively strong scattering from the solvent is 


superposed on the effect in question. 

For this and other reasons I have paid more 
attention to strong solutions (1-10 normal) of 
heavy ions.* Though at first the experimental re- 
sults appeared rather puzzling, I believe that I 
have secured a sufficient agreement with theory 


1 Prins, Zeits. f. Physik 56, 617 (1929). 

2 Prins. Zeits f. Physik 71, 445 (1931). 

_ %Prins, Naturwiss. 19, 435 (1931) (in particular work of 
Krishnamurti quoted there). 

4See for — solutions references 1 and 3 and G. W. 


Stewart, J. C 


em. Phys. 2, 147 (1934). 
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to be able to draw some essential conclusions. 
For different reasons the interpretation is more 
difficult than with simple liquids: (1) the scatter- 
ing from the solvent is always more or less pres- 
ent, (2) a large number of a priori possibilities 
exists for the ‘‘structure of strong solutions,’’ and 
(3) only a very restricted number of heavy elec- 
trolytes possess a sufficient solubility to make 


' their study possible. Because of these complica- 


tions it has proved to be practically impossible 
to work back from the diffraction pattern to the 
structure in the same way as for simple liquids® 
(see, however, §8). 


§2. 


To get a first orientation about the problem we 
shall for the moment suppose that only one kind 
of the ions in solution scatters appreciably 


5 F. Zernike and J. A. Prins. Zeits. f. Physik 41, 184, 1927 
(see also work of other authors quoted there). The method 
indicated in this paper and developed more fully in (1) is the 
same as that applied afterwards by Debye and Menke to 
the same problems (Phys. Zeits. 31, 797 (1930) and else- 
where). It should be observed, however, that in a more 
primitive form analogous ideas occur in older papers of 
different authors; in the first place Keesom, Debye and 
Raman (reference 7). I avoid the word ‘“‘cybotaxis,” though 
it aims at the same conception, because it is often used in a 
sense, with which I cannot entirely agree. For instance the 
assumption of cybotactical ‘“‘groups” with “optimum size” 
and “expanding anisotropically” (Phys. Rev. 37, 9 (1931)) 
appears to me not be to justified by the subject. With this 
I do not mean to detract from the main idea underlying the 
word cybotaxis and the many important experimental 
contributions connected with it. 
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stronger than the water molecules and that there 
is no large difference between the size of this ion 
and that of the water’ molecules (e.g., CsOH or 
Lil). In this case it appears profitable to 
separate the scattering power of the ion into a 
part equal to that of water, say A, and a “‘surplus 
scattering power’’ B. It is clear that the scattering 
power A is in this case distributed over the 
whole volume in much the same way as if it were 
filled with water alone. The amplitude which the 
scattering from this distribution causes at a 
definite point of the diffraction pattern (at a 
certain moment) will be called A, that from the 
rest B. The observed intensity is then pro- 
portional to the mean value of the square of the 
amplitude: 


I= (A+B)?=A?+2AB+B?. (1) 


The first term in the right-hand sum simply 
represents the diffraction pattern of water 


(modified in the sense of Fig. 6). The second term — 


may be considered as due to the interference of 
the scattering from the heavy ion with that from 
the surrounding water molecules. This term is 
easily seen to be of more importance than the last 
term if B is small. As follows from Table II (§8) 
and Fig. 7 this is the case with the alkali-halides 
(§6). On the other hand the last term is more 
important when B is large. This is the case with 
the salts of Th, U and Pb. (§§4, 5). 

The last term B?, of course, represents the 
diffraction pattern due to the surplus scattering 
power alone and it is through this term that the 
arrangement of the heavy solute particles can be 
detected, this being the principal object of this 
investigation. Regarding this arrangement two 
extreme cases may be considered. 

(1) The ions (or molecules) of the solute may 
be distributed in a random* way (as in a gas). 

(2) Under the influence of the mutual forces 
they may tend to form a more or less regular 
arrangement. This arrangement will then be 
superposed on the “‘liquid-arrangement”’ of the 
solvent molecules in much the same way as the 
“superlattice’’ in some alloys is superposed on 


*Even in a gas the distribution is not completely at 
random as the molecules cannot interpenetrate. This 
‘‘geometricalf packing effect” is, however, of minor impor- 
tance throughout this invesiigation, as the ‘‘dilution” is in 
general at least 1: 10. 


X-RAY DIFFRACTION 


IN SOLUTIONS 


Fic. 1 Diffraction pat- 
tern of a gaseous distribu- 
tion (1) and of a “‘liquid- 
arrangement”’ (2). 


Fic. 2. Experimental 
arrangement for x-ray re- 
flection at heavy liquids. 


the finer-meshed “‘lattice’’ constituted by the 
packing of all the metal atoms. This suggests the 
name “‘superarrangement”’ for the arrangement of 
the solute. It will of course presumably be of the 


liquid-type. 


Now in the first case the diffraction pattern 
will be of the type represented by curve 1 in Fig. 
1, while in the second case it is of type 2. Curve 1 
is simply identical with the scattering curve for 
the isolated molecules (or ions), curve 2 is more 
complicated. It differs essentially from 1 in 
possessing a minimum at zero angle and a 
maximum (‘‘amorphous ring”’) at some definite 
larger angle. Both the maximum and the mini- 
mum will be the more pronounced the more 
regular the arrangement is.* Moreover the maxi- 
mum should shift to smaller angles with in- 
creasing dilution (proportionally to its cubic or 
perhaps square root). 

It may be added that the first and second term 
in (1) give rise to curves of type 2, the maximum, 
however, not shifting with concentration. These 
notions will suffice to guide us in the description 
of the experiments. A more complete theory is 
given i. §8. 


§3. EXPERIMENTAL METHOD 


Because of the high absorption coefficients of 
heavy atoms the application of the usual Debye- 
Scherrer method, in which the radiation passes 
through the substance, would require rather 
thin films of liquid. It is possible to make these 
by enclosing the liquid between mica windows? 
which are eventually drawn together by its 
capillary action. Some of the following results 
have been verified by this method. However as it 
is cumbersome to get the most favorable 
thickness (except for light liquids, Fig. 6) the 
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arrangement® shown in Fig. 2 was in general 
preferred. 

The x-rays from a copper target strike the 
liquid surface under a glancing angle of about 3° 
and the reflected pattern is photographed on a 
cylindrical film. The slit (1.5 mm long, 0.4 mm 
wide) is covered by 10u Cu foil, the film by a 
homogeneous 12u Ni foil. These foils serve to 
monochromatize the radiation. The most im- 
portant factor in this respect is, however, that 
the tube is operated at a tension of not more than 
15 kv. The exposition time is 3-10 hours. The 
concentration of the solution is kept roughly 
constant during this time. Sometimes the solution 
was stirred continuously so as to make sure that 
the surface layer was not different in constitution 


from the bulk of the solution. It should be* 


observed that the radiation penetrates into the 
liquid over a distance of at least some microns; 
there are many reasons to assume that so thick a 
surface layer has the same structure as the 
liquid in bulk. 

In this article the results are shown as pho- 
tometer curves. The top line* always indicates 
the galvanometer zero point (100 percent 
blackening). Two scales are given, one for the 
“glancing angle” (i.e., half the diffraction angle) 
g and one for the variable s, defined by.** 


(2) 


It should be observed that with the experi- 
mental arrangement used the blackening at large 
angles is shifted a little to smaller angles as the 
radiation from the first half of the liquid cannot 
reach the film in this case (Fig. 2). This effect has 
roughly been corrected for in the curves by a 
corresponding transformation of the scales, but 
the blackening at, say, s>2 is still considerably 


s=4nsin ¢/\=22/d, where 2d sin 


distorted. In order to locate a large angle maxi- . 


mum more accurately, in some cases the tube 
with attached spectrograph was set at the 
corresponding glancing angle.* No correction is 
then necessary. In general no high precision is 
claimed for the numerical results in this article. 

The absorption in the liquid does not play a 


6 J. A. Prins, Physica 6, 315 (1926). D. Coster and J. A. 
Prirts, J. de Physique 9, 153 (1928). 

* The lower line has ro special meaning. 

** A factor 10° is systematically omitted in the s-values; 
in other words s is expressed in A. 
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large role in distorting the blackening. The 
smaller the glancing angle of emergence, the 
stronger the absorption but also the larger the 
effective surface (see Fig. 2). Only at very small 
angles the relative intensity observed may be a 
little too low, not only because of the absorption, 
but also because the screen (5), which cuts off the 
direct (and totally reflected) ray, may also cut off 
part of the diffracted radiation. 


§4. Resu.ts. Crass I 


We shall first discuss the results for thorium 
and uranyl nitrates, which for our purpose may 
be grouped together as class I. 

As is apparent from Fig. 3, the principal peak 
present in the case of the saturated soiution 
moves to smaller angles and gets less pronounced 
with increasing dilution. This points to the 
conception that the Th**** ions are arranged ina 
liquid superarrangement embedded in the liquid 


arrangement of the water molecules, the inter- 


ionic distances in the superarrangement in- 
creasing and at the same time getting less sharply 
defined with increasing dilution. 

The position of the peaks may also be calcu- 
lated or at least estimated theoretically from the 
number of Th ions per cm* solution. For this 
purpose a knowledge of the density of the 
solutions would, however, be needed and no 
sufficient data seem to exist. I have made a rough 
estimation by extrapolating the volume con- 
traction for analogous ions and found a sufficient 
agreement with the experimental positions. This 
corroborates the above conception. 

It should be pointed out that this conception is 
quite the opposite of the existence of colloidal or 
“subcolloidal”’ crystals of thorium nitrate or 
hydroxide dispersed in the medium. The above 
facts are incompatible with this alternative 
picture of the “structure of solutions’? which 
seems to be held by many. The difference is of 
course that our superarrangement is homo- 
geneously distributed throughout the solution. 
The same seems to me to be true for solutions in 
general ;* for different cases the superarrange- 
ment may of course range from a pronounced 
“‘liquid”’ to a pronounced gaseous distribution. 


_ * Exception must be made for the case of a mixture of two 
liquids at a temperature close to that of complete inter- 
miscibility. 
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Fic. 3. Photometer curves of class I. Diffraction by 
thorium-nitrate solutions. Lower curve, Th (Nos), : 6H,O, 
solid, taken with broad slit. Middle curve, Th(NOs3)4 : H2O 
= 1:15, saturated solution. Upper curve, Th(NOs3)4 : 
=1 : 40, diluted solution. 


These considerations do not exclude the 
possibility that a small fraction of the Th ions 
may be present as colloidal crystals, e.g., of 
Th(OH), and that these may cause the solution 
to show the Tyndall phenomenon. 


Analogous results as with Th(NO3), were found 
inaslightly less pronounced way with UO2(NOs)e. 
In this case the saturated solution is not so 
viscous as in the previous case and the presence of 
colloidal particles is still less probable. 


§5. Crass II 


A more diffuse kind of diffraction pattern is 
shown by saturated and diluted solutions of 
Pb(NOs3)2, Pb acetate, AgNO; and Ba(NOsz)s. In 
this case no sharp maxima are present and no 
marked decrease of intensity is to be seen when 
zero diffraction angle is approached. These 
circumstances make it difficult to investigate the 
influence of concentration on the position of the 
maxima. At any rate, however, this influence 
does not appear to be strong. Typical photometer 
curves are reproduced in Fig. 4. The slight 
decrease at small angles in these curves may be an 
experimental defect (see §3). 

As the diffraction patterns show only weak 
characteristics, they do not afford a good foothold 
for a theoretical explanation. As is shown more 
fully in §9, the difference between this class and 
the previous one may be explained by assuming 


DIFFRACTION 


IN SOLUTIONS 
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Fic. 4. Diffraction by solutions of class II. Upper curve, 
H,O=1: 10. Lower curve, Pb(NOs)2: 
=1: 30. 


that at least a large fraction of the heavy 
constituents are more or less distributed at 
random in the present case. It is reasonabie to 
suppose these heavy constituents to be undis- 
sociated molecules. 

It seems quite probable a priori that undis- 
sociated molecules should be distributed in a 
more ‘‘gaseous’’ way than ions. I feel obliged 
however, to point out that in reality a phe- 
nomenon may be present here which is of a more 
general nature than the difference between 
charged and neutral particles. I mean a difference 
which may crudely be designated by distinguish- 
ing between ‘hygroscopic’ and ‘‘non-hygro- 
scopic” particles. It is at least remarkable that 
the salts of class I crystallize with water of 
crystallization, while those of class II crystallize 
without it. This way of looking at the difference 
in question is not opposite to the previous one. 
On the contrary: it appears quite natural that 
ions should be ‘‘hygroscopic,” while in general 
neutral salt molecules should be less so. Viewed in 
this way, it might be only a matter of termi- 
nology. I think it better, however, to insist on the 
idea of hygroscopy (or hydrophily) because this 
idea may also be applied to some neutral mole- 
cules. Let us first formulate the idea in question: 

“Hygroscopic”’ particles in general crystallize 
with a definite or with an indefinite amount of 
water, their saturated solutions are very viscous 
and may easily be undercooled. These properties 
are linked up in the present argument with a 
tendency of the particles to form a superarrange- 
ment, where every particle is permanently sur- 
rounded by water molecules onlw. 
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As examples of neutral molecules to which this 
idea might apply I suggest fructose (‘‘hygro- 
scopic’) and glucose (not or at least less ‘‘hygro- 
scopic’). This seems to me to find confirmation in 
the x-ray patterns of the solutions.* The well- 
known difference between ‘“‘emulsoidal’’ and 
“‘suspensoidal” sols may be of exactly the same 
kind. 


§6. Crass III 


The diffraction patterns of the alkali-halides I 
have investigated differ radically from the two 
preceding types. They consist in a diffuse ring at 
a considerable angle not shifting appreciably with 
changing concentration (the contrasts, however, 
diminishing a little on dilution). The diameter of 
the ring shows only small variations for different 
substances, being larger for the bromides than 
for the iodides (Fig. 5) and slightly larger for 
RbBr than for NH,Br and LiBr. 

These facts find a natural explanation if we 
assume the ring to be largely due to the inter- 
ference of the scattering from a heavy ion 
(I-, Br~, Rb*) with that from the surrounding 
water molecules. As has been pointed out in §2, 
we may expect this effect, due to the ‘‘mixed 
term” in (1), to be relatively more important 
here than in the previous cases, where the surplus 
scattering power is larger. For the same reason 
the absence of an “‘inner ring”’ in the present case 
does not prove that all tendency to a super- 
arrangement is absent, but only that it cannot be 
very strong.* 

The following numerical data corroborate the 
above explanation. From theoretical consider- 
ations we may expect the diffraction angle of the 
maximum 2¢ to be related to the distance a 
between the center of the heavy ion and that of a 
water molecule touching it by the relation: 


A= 2(fa) sin ¢. 


* It is possible for example that in Fig. 5 the left-hand 
rt of the curve for Li I is partly due to such an inner ring. 
use of its diffuseness, however, it would be hardly 
possible to observe an eventual shift with concentration. 
On the other hand in class I (and perhaps in class II too) 
there seems also to be a large angle maximum independent 
of concentration. With TH(NOD, it is situated at about 
s =3.0. It is reasonable to explain this maximum as due to 
the ‘mixed term” too. In all cases a large fraction of the 
intensity at s=2 is of course due to the term A? in (1) i.e., 
to the water alone. 
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Fic. 5. Diffraction patterns, class II. Solutions of alkali- 
halides. Upper curve: Lil : HXO=1:15 molar. Lower 
curve: RbBr : H.O=1 : 9 molar. 


Here f is a factor for which by slightly different 
ways of reasoning any value between 0.8 and 1 
may be found.’ We shall here somewhat arbi- 
trarily put it equal to 0.85. The preceding 
formula may then be written: 


a=h/1.7 sin ¢ 


or by referring to (2) 
a=7.4/s. (3) 


The values for a calculated by means of this 
formula from our experimental s-values are given 
in Table I (third column). They may be com- 


TABLE I. Distances in A between centers of tons and of 
surrounding water molecules. 


Ion s observed a from (3) a from ionic radii 
3.6=2.2+1.4 
Br- 

Rb* 

Tht+++ 


pared with the values calculated from the well- 
known ionic radii, derived from crystallographic 
data (water molecule r= 1.4A). 

The numerical agreement between column 3 
and 4 is better than our rough measurements and 
calculations give a right to expect and might even 


7 W. H. Keesom and J. de Smedt, Proc. Amst. Acad. 25, 
118 (1922); C. V. Raman and K. R. Ramanathan, Proc. 
Ind. Ass. for Cultiv. Science 8 (2), 127 (1923). P. Debye, 
Phys. Zeits. 28, 135 (1927). See also references 6 and 1. 
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be slightly fortuitous. But the parallelism be- 
tween the two columns can hardiy be accidental 
and is a strong support in favor of the interpre- 
tation of the results advocated in this paragraph. 

The preceding diffraction patterns of class ITI 
have in common with those of class I that the 
intensity decreases markedly when the diffraction 
angle tends to zero. It cannot be assured experi- 
mentally that the intensity tends to zero and 
indeed it follows from theory (§$§$8, 9) that only a 
certain not very small limiting value should be 
approached (Fig. 8). Nevertheless the experi- 
mental decrease may be said to be clearly marked. 

This feature gets more and more lost when the 
following series is successively investigated: 
(class II1), CsNOs;, SrBre, BaCl., Ba(NOz)s, 
Ba(NOs;)2, Ba acetate, (class II). This suggests 
that the salts between classes III and IT must be 
considered as transition cases from class III to 
class II, in good accord with their position in the 
periodic system. 


§7. Licgut Ions AND WATER 


To make sure that the small angle intensity 
observed with Pb(NOsz)e etc. ($5) is not due to 
the presence of the NO;~ ‘ion alone I have also 
photographed the diffraction pattern of a NaNO; 
solution with a mole ratio of 1:10. No small 
angle scattering of any importance was observed. 
The pattern resembled that of an equivalent 
NaCl solution. This pattern may be interpreted 
as mainly due to the arrangement which the 
water molecules assume in the electric field of the 
ions, these being too light to contribute ap- 
preciably to the diffraction pattern themselves. 

This conception of a rearrangement of water by 
ions is due to Bernal and Fowler* who made use 
of experimental results of Meyer.’ According to 
this conception in pure water the H.O molecules 
are arranged in a ‘“‘tetrahedral’’ and more 
especially in a ‘‘quartz-like”’ liquid arrangement. 
This arrangement is due to the peculiar electrical 
constitution of the H2O molecules consisting of 
two positive and two negative poles arranged 
‘etrahedrally in the (almost spherical) surface of 
the molecule. This constitution causes every 


tae Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 
1933). 
*H. H. Meyer, Ann. d. Physik 5, 701 (1930). 


DIFFRACTION 


IN SOLUTIONS 


Fic. 6. Debye-Scherrer photographs taken with Cu Ka 
radiation penetrating a liquid layer 0.5 mm thick between 
mica windows 3u thick. Upper figure, pure water. Just out- 
side the main ring a weaker second ring is present. Lower 
figure, LiOH : H,O=1 : 10 (saturated). The two rings have 
fused to one of slightly larger diameter than the main ring 
of pure water. 


molecule to be surrounded by only four neighbor- 
ing molecules which join it each with one of their 
poles (+ to — of course). This tetrahedrical 
arrangement is, however, wrenched loose by the 
strong electric field of the ions, if these are 
introduced. A more or less close packing, where 
every molecule has more than four neighbors, is 
the result. This not only accounts for the well- 
known volume contraction observed on adding 
an electrolyte to water but also for the changing 
of the peculiar x-ray diffraction pattern of water 
to the more universal pattern of a close-packed 
liquid (see Fig. 6). 

Fig. 6 represents some experiments I have 
made with solutions of NaOH and LiOH to 
confirm and extend the analogous results of 
Meyer, who used NaCl and LiCl! solutions. To 
these the objection might be raised that the 
scattering by the Cl- ion contributed appreciably 
to the modified pattern. This objection may be 
shown not to be very serious, as has already been 
observed by Meyer, but I have thought it would 
be more convincing to show the effect with 
LiOH. In this case the objection loses all 
significance as the scattering by the Li* ion is 
negligible and that of the OH~ ion may well be 
identified with that of a H2O molecule. 

Incidentally it may be mentioned that I have 
the impression that 5 normal solutions of LiOH 
and NaOH do not produce yuite so large a modi- 


| 
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fication of the pattern as equivalent solutions 
of LiCl and NaCl. 


$8. THEORY 


Without an exact theoretical treatment any 
discussion would become too vague. We shall 
make use of the theory given by Zernike and 
Prins® for ordinary liquids, as it is very well 
suited for the present case of solutions. In this 
theory the arrangement of the molecules and ions 
is expressed by a distribution law g(r) for the 
distances between the molecular (or ionic) centra. 
As was shown by the authors it is then possible to 
calculate the x-ray pattern from the distribution 
function and vice versa.’ 

For the sake of simplicity we shall begin by 
supposing the water to contain only one kind of 
heavy ions, or stated more accurately: We sup- 
pose a certain volume to contain N particles of 
scattering power (see below) A and n particles of 
scattering power a; the particles of the first kind 
will be called (water) molecules, those of the 
second kind ions. The scattering by the volume in 
question will then at some definite point of the 
diffraction pattern give rise to an amplitude: 


where the separate terms indicate the amplitudes 
of the elementary waves from the different 
particles. Now what is wanted is the intensity, in 
other words the mean square of the above 
quantity: 


This works out to a sum of all possible products, 
occurring in three types: 


Va. 
NN Nn nn 


The mean value, e.g., of the first sum may now be 
found by first taking together all terms having 
A; as first factor, then taking the mean of this 
aggregate keeping A; constant and finally taking 
the mean over A,. This gives 


=A?+ a (sr) /sr J. 
0 


PRINS 


Here A denotes the mean* amplitude of the 
elementary wave scattered by one molecule. For 
spherically symmetrical molecules A?= A * and so 
we may designate both quantities with A?. This 
then is the exact definition of the ‘‘scattering 
power” A introduced above. As regards the 
integral in the right-hand expression it is 
equivalent to the mean sum over / (excluding 
/=k) figuring in the left-hand expression. The sin 
(sr)/sr and the distribution function g(r) enter 
into it in virtue of the phase differences existing 
between A; and A; for all possible distances and 
diffraction angles. Finally the second summation 
(over k) now yields a factor N (or n in the third 
sum) as all indices have disappeared. So finally: 


sin 


1=wayi+ dr4rr°g(r) 
0 


sin (sr) 
+2NAa dr4rr*g’(r) 
0 sr 


0 


or, designating the three integrals, respectively, 
by G(s), G’(s), G’’(s): 


I= 
+na*[1+G'"(s)]. (5) 


In these formulae primed letters refer to the 
distribution of ions with respect to a water 
molecule; doubly primed letters refer to the 
distribution of ions with respect to one of them. 

If the case of a 1:10 Lil solution is con- 
sidered as an example and the presence of the 
Lit ions is disregarded, then one may put 
n= N/10, a=5A (see Table II). This renders the 
second term more important than the third, 
which explains the remarkable features of the 
diffraction patterns of class II, mentioned in §4. 
On the other hand with higher values of a the 
last term gets more important. Jf the distribution 
function g” corresponds to a superarrangement 
this last term may be held responsible for the 
shifting peak found with class I (§4). Jf on the 
other hand the distribution expressed by g’’ is 
more or less a gaseous one, the presence of the 


*In (4) A means the maximal amplitude. This causes 
3 NA? to appear where here in (5) NA? appears. 
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H,0 Br r NO; 


Li? 

Fic. 7. Dimensions and scattering power of some ions. 
If the shaded volume is supposed to contain an electron 
density 10 times that of water it will give the right surplus 
scattering power. 


last term accounts largely for the features of class 
III, at least in the case of lead salts. In the case of 
AgNO; the preponderance of the third term is 
dubious and a closer examination (§9) is needed 
to check the agreement between theory and 
experiment. 

In all cases our principal aim is to get a 
knowledge about g’’. This may be derived from 
an accurate determination of 1+G”. This object 
may be facilitated by comparing the diffraction 
pattern for a heavy electrolyte such as CsOH 
with the one for an analogous light electrolyte 
such as NaOH. If the two diffraction patterns, 
reduced to the same value of N, are subtracted, 
we should expect the first term to disappear, 
while the second term might perhaps be esti- 
mated with sufficient accuracy to find the last 
term. This method was alluded to in the last 
part of §1. 

In Fig. 7 and in Table II numerical data on the 


TABLE II. Dimensions and scattering power of some ions 
and of H.0. 


H:O Br IJI- NOs Rb* Ag* Tht+++ 
Radius in A 40 1.95 2.20 1.3-2.6 1.49 1.13 1.32 1.10 
35 53° 32 35 45 74 80 


1. 
Effective number ofelectrons 10 
Surplusnumber ofelectrons 0 7(15) 14(25) 2(11) 23 40 66 75 
0.58 0.73 0.5 0.85 1.03 1.22 1.27 


Radius broken circle Fig. 7 0 


scattering powers and dimensions of the ions, 
discussed in this article, are exhibited. The 
scattering power is represented by the “effective 
number of electrons” which may be found from 
the total number of electrons by eventually 
taking into account the anomalous scattering of 
K and L electrons. For some purposes it is more 
useful to know the ‘“‘surplus scattering power,” 
obtained by diminishing the previous quantity by 
the amount of electrons that would be needed to 


fill the volume of the ion with the same electron 
density as is present in a water molecule.* It is 
easy to transcribe formulae (4) and (5) under the 
assumption that a denotes the surplus scattering 
power and N the ¢éotal number of particles. The 
first term should then almost completely coincide 
with the /ower diffraction pattern exhibited in 
Fig. 6. 


§9. ONE-DIMENSIONAL MODELS 


We shall assume in this paragraph that no 
forces exist between any particles in the liquid 
except of course insofar as they keep it together 
in a definite volume. It is then possible in 
principle to calculate the diffraction patterns 
from (4), as the g-functions are completely 
defined by this assumption and by the dimen- 
sions of the particles and the available free space. 
In practice it is, however, too difficult to perform 
this calculation in the three-dimensional case. 
The analogous calculation is, however, easily 
made rigorously in the one-dimensional case and 
the result may be expected to provide a rather 
good picture of the corresponding three-dimen- 
sional case. We shall show in this way that the 
diffraction patterns of class II and III may be 
understood without introducing new principles. 

Consider a straight line of length L con- 
taining a large number aN of “molecules” of 
length a and scattering power A and a large 
number BN(a+6=1) of “ions’’ of length 6 and 
scattering power B. The total free space L—aNa 
— Nb divided by N will be called the mean free 
space per particle /. On the assumptions stated 
above the arrangement is a matter of geometrical 
probability. By an argument analogous!’ to 
that used by Zernike and Prins® the following 
expression may be derived for the diffraction 
pattern of the arrangement (considered as section 
of a plane “grating’’ at normal incidence) : 


I/N=real part cf aA*+ 6B? 
2(aA 6) 
1—ils — (ae***+ 


* Instead of this density the electron density of close- 
cked water in bulk should be taken if the ion is relatively 


rge. The corresponding figures have been added in paren- 
theses in Table II. The truth lies somewhere between these 
numbers and those preceding them. 

10 The mathematical treatment may be simplified by 
introducing generating functions for the intermolecular 
distances. 
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Fic. 8. Theoretical diffraction pattern for one-dimen- 
sional I~ solution (1:10). Broken curve represents ex- 
pression (6). Full curve is derived from it by taking into 
account the ‘‘atomic scattering powers’’ for H.O and I- 
from the Thomas-Fermi model. 


In this expression s=2z sin (~=diffraction 
angle), this quantity is the analog of the variable 
s from (1) in the three-dimensional case. 

To get the case of a 1 : 10 I- solution we put: 
a=0.9, B=0.1, b=2a, =a, B=5A, a+/=2.8A. 
The result is shown in Fig. 8, which should be 
compared* with Fig. 5. 

To discuss expression (6) we may separate it in 
four terms: (1) The term aA?+ 8B? which will be 


called the ‘‘gas-term.”’ It is the most important: 


term for large values of s. (2) The first term 
arising from working out the quadratic numer- 
ator contains the factor A? and will be called 
the “water-term.”’ (3) The second term arising in 
the same way and containing AB will be called 
the ‘‘mixed term.” (4) The third term arising 
in the same way and containing B? will be called 
the “ionic term.” 

In our case the mixed term is slightly more 
prominent than the water term and both terms 
are more prominent than the ionic term. If the 
mixed term alone were present the position of the 
maximum would roughly correspond to a 


* For different reasons the s-scales cannot be expected to 
correspond exactly. A satisfactory agreement may, how- 
ever, be established in this respect, too, by a more careful 
elaboration of the correspondence. 


“grating constant’’ a+/+-)+/ in other words to 
the mean distance between the centers of an ion 
and a neighboring water molecule. With in- 
creasing dilution the water term gets more 
predominating and the maximum slowly shifts to 
the position corresponding to the distance be- 
tween adjoining water molecules. This shift is, 
however, so slow that it is no wonder that it has 
not been noted experimentally (§6) with con- 
centrations varying at most from 1:5 to 1 : 20. 
With the other alkali-halides the difference be- 
tween the dimensions of the ion and the water 
molecule is too small to give rise to an ap- 
preciable shift. 

The same model may more or less serve for the 
AgNOs3 and Pb(NOs)2 solutions when the salt 
molecule is substituted for the I ion. It has in 
both cases roughly the same dimensions as the 
I ion. The principal difference consists in the 
following two points: 

I. The scattering power of the salt molecule is 
larger than that of the I ion. , 

II. The “molecular scattering curve’’ for the 
salt molecule shows a steeper decrease beginning 
at a smaller angle than the “ionic scattering 
curve”’ of I-. 

Circumstance I renders the gas term 1 (see 
above) relatively more prominent. This makes 
the minimum at small angles less marked. 
Circumstance II has an analogous effect. A 
rough calculation shows that even with AgNO; 
(and a fortiori with Pb salts) the practical 
disappearance of the minimum found experi- 
mentally (Fig. 4) may be well explained in this 
way. 

The shifting peak of class I can of course not 
be explained by the model of this paragraph. It 
is due to special forces (§5) causing a super- 
arrangement. On the other hand we may now be 
pretty sure that an analogous superarrangement 
is not present to an appreciable degree with class 
II (and that with class III it cannot be very 
strong) as the experimental results appear to be 
explained without introducing it. 

At,this stage I believe that only very accurate 
(and if possible absolute) measurements could 
provide more precise information. It is certainly 
disappointing that the elements 85 and 87 are not 
available for this purpose. 
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The energies of a number of hydrocarbon molecules have 
been calculated by the Heitler-London-Pauling-Slater 
method with two objects in view, first, to see if this simple 
theory leads to the correct energies of formation, and 
second, to determine what significance can be attached to 
the empirical additivity rule for bond energies. Our equa- 
tions involve three parameters, which roughly speaking, 
determine the gross energies of a C—H bond, a C—C bond, 
and a C=C bond. These have been evaluated from the ob- 


served energies of methane, ethane and ethylene. The 
energies calculated for other hydrocarbons then agree sur- 
prisingly well with the observed values. Nevertheless, the 
validity of an empirical bond additivity rule seems to be 
purely fortuitous. The calculations further indicate that 
the theoretical significance ordinarily ascribed te the so- 
called “observed resonance energies’ is a very question- 
able one. 


AN VLECK,?:*: Penney® * and their 
collaborators have applied the Heitler- 
London-Pauling-Slater theory with considerable 
success in explaining the structure of polyatomic 
molecules. While the theory has thus shown its 
adequacy for such semiqualitative considera- 
tions, it would indeed be optimistic to expect of 
it an equally good account of the total energies, 
that is, the energies liberated in forming mole- 
cules from widely separated atoms in their 
normal states. Our calculation of the total 
energies of a number of hydrocarbon molecules 
thus affords a much more severe test of the 
approximations involved in the H-L-P-S pro- 
cedure. The calculation is additionally interesting 
since through it one may hope to obtain a 
theoretical explanation of the empirical bond 
additivity rule. 

Let us suppose that a hydrogen orbit, a, is 
paired to a carbon orbit, a. The total interaction 
energy between the hydrogen and carbon atoms 
is then 

Wuc=3/2Jaat+Suc, (1) 


Suc= Quc— 1/2 


i=s,o, 7, 


where 


Here the J’s are the ordinary exchange integrals, 


1A preliminary account was given at the Ann Arbor 
Meeting of the American Physical Society, June, 1934. 
(Phys. Rev. 46, 335A (1934)). 

* National Research Fellow. 

2 J. H. Van Vleck, J. Chem. Phys. 1, 177 (1933). 

3 J. H. Van Vleck, J. Chem. Phys. 1, 219 (1933). 

4 J. H. Van Vleck, J. Chem. Phys. 2, 20 (1934). 

5 J. H. Van Vleck, J. Chem. Phys. 2, 297 (1934). 

6 W.-G. Penney, Proc. Roy. Soc. A144, 166 (1934). 

7W. G. Penney, Proc. Roy. Soc. A146, 223 (1934). 
(1934) G. Penney, Proc. Phys. Soc. (London) 46, 333 
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and Qyc is the Coulomb energy. The steric 
term, Syc, is a function only of the internuclear 
distance and is quite independent of the choice 
of carbon orbits,’ i.e., of the angles between the 
bonds. Advantage has been taken of this in- 
variance to write Syc in terms of unhybridized 
orbits. The individual J’s, on the other hand, do 
depend on the bond angles. If the hydrogen 
orbit is not paired to any orbit on the carbon 
atom the first term in (1) drops out and only 
the steric term is left. In a similar way, the 
interaction between two carbon atoms is 


Wee'= 3/2 Jivt+Sco’, (2) 
pairs 


where Soo =Qeco—1/2 


j=8, 0, 7, 


J ij. 

To the approximation of pure electron pair 
bonding, the total energy of the molecule is 
W= > +L Wewe;+ W.(C;) 

ij i<j 
+) M(H.H)). 


i<i 


(3) 


The last sum in (3) represents the hydrogen- 
hydrogen repulsions. W,(C;) is the energy of 
the valence state of the ith carbon atom, the 
normal state of the carbon atom being taken as 
the zero of energy. Because of the influence of 
the L—S structure of the carbon atom W,(C;) 
is not necessarily the same for all carbon atoms 
in the molecule, but depends on the angles 


* This is a simple consequence of the unitary character of 
the transformation from one set of orbits to another. A 
proof is given by Penney, refererice 7. 
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between the bonds made by the atom in question. 
In the cases we will consider, W,(C;) lies between 
7 v and 7.5 v. For any particular molecular 
model, the third and fourth sums of (3) can 
readily be evaluated by means of formulas 
given by Van Vleck.’® The first and second sums 
of (3) contain terms of two. types: interactions 
between adjacent atoms, i.e., between atoms 
which are bonded to each other, and steric 
interactions between distant, unbonded atoms. 
It will be seen from (1) and (2) that the steric 
terms consist of sums of terms of opposite sign 
which largely annul each other. The steric 
terms, and Wec:, between distant atoms 
are certainly small, and we shall provisionally 
suppose that they are so small that they may be 
neglected. If we obtain good agreement between 
the calculated and observed energies, we can 
regard this supposition as substantiated. The 
terms Wyc between adjacent atoms can be 
calculated by using Van Vleck’s‘ values of the 
exchange integrals N,,, Nss, Nes, and by de- 
termining Syc from the observed energy of 


methane. We shall neglect the small differences © 


in internuclear distances in different molecules. 
The Wcec: terms between adjacent carbon atoms 
appear at first glance much more troublesome, 
since they involve numerous exchange integrals, 
Cagys in Penney’s notation,® whose values are 
not reliably known. The problem is greatly 
simplified, however, by the following considera- 
tions, which we shall illustrate by reference to 
Penney’s model of the ethylene molecule.® In 
this model Wyc is the same for all four hydrogen 
atoms, and is a function of the angle w (see 
insert in Fig. 1. Of course for ethylene the lower 
C atom is also bonded to two H atoms in the 
same manner as the upper one). Wee can be 
written 


Wee'= We-c'+3/2a, (4) 


where a is the exchange integral J,, in the 
notation of (2), or —C,+** in Penney’s notation. 
The first term in (4) is the bonding energy of 
the principal carbon-carbon bond (the one di- 
rected along the C—C line), the second term 
represents the energy of the second carbon bond 
(the r—-7 bond), which is independent of w and 


10 Reference 4. Note corrections given in reference 5. 
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w 

Fic. 1. In computing We_c: it has been supposed that 


the hybridization ratios of the two orbits in the C—C bond 
are equal, an obvious condition in order that Wo_c- be a 


minimum. 


of no interest for our present purposes. In Fig. 1 
we have plotted Wyc as a function of w. It will 
be seen that the curve is very flat over a wide 
range of values of w. Hence if we suppose that 
We-c has a sharp minimum in this range, the 
carbon-hydrogen interaction plays very little 
part in determining the form of the molecule: 
the hydrogen bonds simply arrange themselves, 
in directions which bring the carbon single bond 
energy to its minimum. The stability of the 
benzene ring immediately suggests that the 
minimum of We_c should fall at w=30°; we 
shall assume this to be so in our further calcula- 
tions. One may test this assumption by using 
the values of the various exchange integrals, 
Capys, Suggested by Penney." The variation of 
We-c with w found in this way is shown in 
Fig. 1. We see that with this choice of exchange 
integrals, We-c does indeed’ have a sharp 
minimum at 30°. 

For Van Vleck’s @ type bonds” the curves 
obtained for Wyc and We_c are again very 


1 Reference 7. The set of values labeled (2) has been em- 
ployed since Dr. Penney informs the writer that these are 
probably the best estimate. 

12 See Fig. 2 and Eq. (33), reference 3. In this case the C 
atom is bonded to three H atoms and another C atom. If 
the C—C line is taken as the polar direction, the azimuths 
of the three H atoms are 0°, 120°, 240°, respectively, rather 
than 0°, 180°, as shown in Fig. 1, and @ is the angle CCH. 
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similar in appearance to those shown in Fig. 1. 
It thus appears that for all molecular models 
involving bonds of the w and @ types, maximum 
stability is attained when the various bonding 
angles are such that Wge_c, for each carbon- 
carbon bond, takes on its minimum value.” 
In order to find the energy of the principal 
carbon-carbon bond we accordingly need not 
know the values of the individual carbon-carbon 
exchange integrals, Cag,s but only that of the 
particular linear combination of them involved 
in We_c: at its minimum. Hereafter, Wce_c: is to 
be understood as representing the minimum 
value of this quantity, unless otherwise stated. 
We shall determine Wce_c from the observed 
energy of ethane. The only other carbon-carbon 
exchange integral which plays an important role 
in our calculations is a, which, as (4) shows, 
fixes the energy of the second bond in a carbon- 
carbon double bond. a will be evaluated from 
the observed difference between the energies of 
ethane and ethylene. 

Before we proceed to the explicit calculations, 


one troublesome point must be examined. The ~ 


observed molecular energy is the sum of the 
“electronic” energy and the energy of nuclear 
motion. Our calculations should refer only to 
the electronic part of the energy. However, it 
is not usually possible to actually compute the 
nuclear energy. This can be done approximately 
in a few cases. Zahn" estimates that the nuclear 
energy of methane is 1.25 v, that of ethane is 
1.97 v, and each additional CHe group increases 
the nuclear energy by about 0.7 v. In a similar 
way, using the data of Mecke” and of Eucken 
and Parts!* and Sutherland,'’ one finds the value 
1.4 v for ethylene and 0.7 v for acetylene. It will 
be seen that these figures can be represented 
fairly well by a rule of the bond additivity type 
which assigns a nuclear energy of about 0.3 v to 
each C—H bond, and an energy of about 0.1 v 


3 If, with w type bonding, the carbon atom is bonded to 
two other carbons and one hydrogen, rather than to two 
hydrogens and one carbon, the influence of the H—C 
interaction is considerably more important than indicated 
by Fig. 1. Fortunately, in this case the H—C interaction, 
as well as the C—C interaction, favors the 30° configura- 
tion. See reference 7. 

4C,. T. Zahn, J. Chem. Phys. 2, 671 (1934). 

1% R. Mecke, Zeits. f. shell Chemie B17, 1 (1932). 
usa and Parts, Zeits. f. physik. Chemie B20, 184 

17 Quoted by Penney, reference 8. 


to each carbon-carbon bond. We must suppose 
that this representation of the nuclear energies 
is sufficiently accurate for our purposes. The 
nuclear energies then need not be taken into 
account explicitly: their only effect will be to 
increase the “observed” value of Syc by 0.3 v, 
and of Sco’ by 0.1 v. 

From the observed energy of methane, W(CH,) 
= —16.3 v, we find, by a calculation similar to 
that given by Van Vleck,‘ that Wyc for methane 
is —6.47 v. Once this value is determined it is 
a simple matter to compute Wyc for any other 
arrangement of carbon orbits. 

For ethane we shall use the model proposed by 
Penney.® In order to make We_c a minimum. 
we must take @=112°, this being the value of 6 
which gives the same hybridization ratio in the 
C—C bond orbits as does w= 30°. The energy of 
ethane is then 


W(C2He) = 6Wyc(0= 112°) + 
+2W,(C)+ M(HiH)). 


i<j 


(5) 


The values of the various terms are: Wyc(@ 
=112°)=—6.42 v, W,(C)=7.06 v, >> M(H.H;) 
=3.72 v. In computing >> M(H;:H;) we have 
taken the carbon-hydrogen internuclear distance 
to be 1.08A, and the carbon-carbon distance to 
be 1.37A. From the observed energy, W(C:H¢) 
= — 27.57 we find Wo-c'= — 6.89 v. 

The difference in energy between ethane and 
ethylene is 


W(C2H¢) W(C2H,) = 112°) 
—4Wyc(w= 30°) —3/2a+2LW,(C) —W,(C’) ] 
+2 M(H.H,;) J. 


The primed terms here refer to ethylene. Us- 
ing Wyc(w=30°)=—6.59 v, W,(C’)=7.20 v, 
[> M(H.H,) =1.09 v, W(C:H.)— W(C:H,) = 
—5.73 v, we find a=—2.72 v. It should be 
noted that this value of a@ is entirely compatible 
with Penney’s calculation of the twisting fre- 
quency of ethylene,® since his calculation de- 
termines not a, but a—a’, where a’= —Crras 
(a’=J,, in the notation of (2)). Pauling and 
Wheland!* give a= —1.5 v, but, as we shall see, 
@ 


18 L. Pauling and G. W. Wheland, J. Chem. Phys. 1, 362 
(1933). 
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the argument which leads to their value is a 
very dubious one. 

It is now a simple matter to compute the 
energies of hydrocarbons for which simple elec- 
tron pair bonding is a reasonable approximation, 
i.e., which do not involve resonance effects. 
In our models of molecules of the ethane and 
ethylene series, the carbon atoms are arranged 
in a staggered chain. The staggering is of course 
necessary to preserve the 120° angles between 
the carbon-carbon bonds. For example, for 
propane, the bonding of the two end carbon 
atoms is of the @ type, with @=112°, while the 
bonding of the center carbon atom is of Van 
Vleck’s class (a) type,!® with (180°—2w), the 
angle between the carbon-carbon bonds, equal to 
120°. The energy of propane is thus”? 


W(CsHs) = 6 112°) +2 39.3°) 
+2W2(Ci1) + W.(C2) M(HH)), 
with 
Wauc(w= 39.3°) = —6.377, W,(C2)=7.10 v, 
> M(H:H;)=4.90 v, 


and the remaining terms as given in (5). This 
gives W(C3;Hs) = — 38.92 v. The observed value 
is — 38.9 v. The results of the calculations for a 
number of molecules are given in Table I. 

The approximation of simple electron pair 
bonding is inadequate for molecules such as 
acetylene and benzene which involve multiple 
a—am bonds. Here more accurate methods must 
be employed. For acetylene, we use a linear 
model. If the axis of figure be the z axis, each 
carbon atom has a 7 orbit directed along the x 
and y axes. The interaction between these four 
orbits must be considered as a four-electron 
problem, rather than as two electron pairs with a 
small interaction term. The energy can be written 
W(C2H2)= Wot Wi, where Wo=2Wyct 
+2W,.(C)+ M(HiH:2), and W, is the energy due 
to the four a orbits. For this model Wyc= 
—6.61 v, W.(C)=7.51 v and M(H;,Hz2) is neg- 
ligible. Unlike the other cases we consider, for 
acetylene Wce_c does not have its minimum 
value. In order to estimate its deviation from 


19 See Fig. 1 and Eqs. (32) and ‘34), refe:ence 3. 
20 The value w=39.3° is determined from the constraint 
between the angles given by Eq. (32), reference 3. 


the minimum we have used Penney’s values of 
the carbon-carbon exchange integrals. The devia- 
tion is found to be 0.33 v, thus We_¢- = — 6.56 v. 
The familiar four-electron formula?! gives 


tata’ (6) 


where a! = — =J ex’, and is the intra- 
atomic exchange integral between the 7 orbits. 
The ordinary four-electron formula comprises 
only the term under the radical sign in (6); 
the necessity of adding the remaining terms can 
be seen as follows. The term under the radical 
sign represents the complete contribution to the 
energy from the four 7 orbits, hence W, should 
contain no contribution from these orbits. How- 
ever, if we examine W, in detail we see from (2) 
that Scc contains, for each choice of a pair of 
m@ orbits, one on each carbon atom, either the 
term — 3a or the term —}a’, depending on the 
relative orientation of the orbits. The total 
contribution to Scc of our four 7 orbits is thus 
—a-—a’. Similarly, from Van Vleck’s formula 
for W.,(C), we see that the contribution of our 
orbits to W,(C) is The term ata’ 
+K,, must thus be added to the radical in (6) 
in order to correct for these steric terms which 
appear in W>. A simpler example of the same 
point is afforded by noting from (2) or (4) that 
the difference in energy between a carbon single 
bond and a carbon double bond is 3a, not a as 
one might at first suppose. The energy of the 
ma—m bond is indeed a, but it must be re- 
membered that if this bond were lacking, the 
interaction of the two z orbits would contribute 
a term —4a to the energy, and this no matter 
what the choice of carbon orbits. The difference 
in energy due to the presence of the r—7 bond 
is thus $a. A similar interpretation holds for the 
added terms in (6). 

According to Van Vleck,’ K,.,=0.56 v. The 
exchange integral a’ can be obtained from the 
difference a—a’=—0.72 v deduced by Penney 
from the observed twisting frequency of ethylene. 
Using these values we find W(C2H2)= — 14.9 v, 
compared to the observed value, — 15.2 v. Since 
most of the stability is due to Wi, Wo amounting 


21 F, London, Zeits. f. Electrochemie 35, 552 (1929); J. C. 
Slater, Phys. Rev. 38, 1109 (1931). 
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to only —4.76 v, the agreement between the 
calculated and observed values is a very nice 
confirmation of the value we have found for a. 

The interaction of the 7 orbits in benzene 
must be treated as a six electron problem. 
The energy is W= W.+ W,, with 


Wo= 6Wacl(w= 30°) +6We_ec: +6W,(C) 
+2 M(HiH)),’ 


and”? 
(7) 


Here 6 and y are the values of — C,,,, at inter- 
nuclear distances corresponding to the distance 
between alternate carbon atoms, and the distance 
between diagonally opposite atoms, respectively. 
The steric correction term 3a+38+3y/2 is 
included in W;. An estimate of the magnitude 
of 6 and y can be obtained from Bartlett's” 
theoretical calculation of C,,,,. While the abso- 
lute magnitude of C,,,, is not given correctly 
by this calculation, it is not unreasonable to 
suppose that the ratios of the values of C,+r, 
at different internuclear distances are fairly well 
represented. Using Slater’s™ value of the screen- 
ing constant for the mz orbits, one concludes that 
8 and y are quite small and we have accordingly 
neglected them in computing the value of 
W(C.He) given in Table I. It should be empha- 
sized that in computing the energies of molecules 


TABLE I. Energies of hydrocarbons. 


Calc. 
(16.3)v 


Obs. 


16.3v 
27.6 
38.9 
50.2 
21.8 
33.1 
61.8 
67.1 


15.2 
52.2 
56.7 
82.5 


Methane 

Ethane 

Propane 

Butane 

Ethylene 

Propylene 

1,2,3,4 tetrahydrobenzol 
Cyclohexane 


Acetylene 
Benzene 

1,2 dihydrobenzol 
Naphthalene 


2 F, Seitz and A. Sherman, J. Chem. Phys. 2, 11 (1934); 
R. Serber, J. Chem. Phys. 2, 697 (1934). 

3 J. H. Bartlett, Jr., Phys. Rev. 37, 507 (1931). 

*4 J.C. Slater, Phys. Rev. 36, 57 (1930). 

°° The experimental values have been taken from the 
following sources: F. D. Rossini, Bur. Standards J. Re- 
search 13, 21 (1934); L. Pauling, J. Am. Chem. Soc. 54, 
3570 (1932); and references 14, 28, 30. 


such as benzene it is not sufficient to consider 
merely the interaction between 7 orbits. It is 
imperative also to include the forces in the plane 
of the benzene ring, as we have done. The im- 
portance of bond forces in the plane of the ring 
has been shown by Penney ;’ the importance of 
the steric terms which appear in W, has been 
discussed above in connection with acetylene. 
The errors which attach to incomplete treatments 
which consider only z orbits will be illustrated 
in the discussion of results. 

On taking account of the interaction of the 
four 7 orbits in 1,2 dihydrobenzol, we find 
W,=3.23a. For napthalene we have used Sher- 
man’s* result for the resonance energy. This 
gives W,=9.54a. 


DISCUSSION OF RESULTS 


The agreement between the calculated and 
observed energies is remarkably good, much 
better than one has any right to expect from so 
crude a theory. It should be pointed out that 
inclusion, in the H-L-P-S treatment, of non- 
orthogonality and higher order permutation in- 
tegrals generally seems to spoil the results of 
the simple theory.”’ Any agreement obtained by 
use of the simple theory must consequently be 
regarded as rather accidental. 

It is not at all clear why molecules for which 
(3) is a good approximation should so nearly 
obey a bond additivity law. The variation in 
Wuc from molecule to molecule is over } v, 
while W,(C) varies by almost 3 v. In addition, 
the hydrogen-hydrogen repulsion energy is often 
quite large, and it by no means obeys an 
additivity rule. The situation is not clarified by 
considering net, rather than gross, bonding 
energies. For example, the net C—H _ bonding 
energy in ethylene is 4.45 v, in ethane 4.04 v, 
and in cyclohexane 3.95 v. In computing these 
net bonding energies, the energy of the hydrogen- 
hydrogen repulsions has been divided equally 
among the C—H bonds. The energy of the C—C 
bond in ethane is 3.36 v, and that of the C=C 
bond in ethylene is 3.77 v. The empirical bonding 
energies given by Pauling and Sherman” are 


ad {; Sherman, J. Chem. Phys. 2, 488 (1934). 

27H. M. James, J. Chem. Phys. 2, 794 (1934); A. S. 
Coolidge and H. M. James, ibid. 2, 811 (1934). 
(1939) Pauling and J. Sherman, J. Chem. Phys. 1, 606 


CH, 
(27.6) 
C;Hs 38.9 
CyHio 50.2 
CoH, (21.8) 
C3H¢ 33.1 
CoH io 61.9 
CH» 
14.9 
52.2 
56.3 
CyHs 81.4 
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(C—H)=4.07 v, (C—C)=3.15 v, (C=C) 
=5.56 v. The wide divergences between these 
values force one to the conclusion that the 
validity of an empirical bond additivity law is 
purely fortuitous. The bond additivity rule does, 
in fact, overstabilize cyclohexane and 1,2,3,4 
tetrahydrobenzol by } v and } v, respectively. 
The H-L-P-S method thus gives a somewhat 
better account of the energies of these molecules 
than the additivity rule, although, of course, 
the accuracy of our treatment is hardly sufficient 
to justify the comparison. 

Although no general bond additivity rule 
should be expected to hold for molecules whose 
energies are given by (3), the theory does predict 
an additivity relation between the energies of 
the higher members of the paraffin series, 
C,,Hense. This follows from the fact that the 
bonding of each added CHe group is of the same 
type, namely Van Vleck’s class (a) type, while 
for n>6 the repulsions between the two added 
H atoms and the H atoms on the two end C 
atoms become negligible. This agrees qualita- 
tively with the results of Rossini,» who finds 
an additivity relation for the paraffins when 
n>5, but slight deviations for the lower members 
of the series. 

_ Before discussing benzene, a few words re- 
garding Pauling, Wheland and Sherman’s”® theo- 
retical formulas for resonance energies may not 
be amiss, since the precise significance of their 
method of calculation is by no means clear. 
Thus Hiickel*® has argued that if benzene were 
considered to have three single carbon-carbon 
‘bonds and three double bonds the energy would 
be W’=W,+3a, whereas solution of the six 
electron problem would give W= W 5+2.605a. 
The resonance energy W—W’ is consequently 
asserted to be —0.395a, and not 1.105a, as 
stated by Pauling and Wheland. But this argu- 
ment is evidently not correct. In fact, we see 
from (2) or (4) that W’=W +4$a, while (7) 


29 Pauling and Wheland, reference 18; Pauling and 
——. J. Chem. Phys. 1, 679 (1933); Sherman, reference 
2 


oF, Hiickel, Aromatic and Unsaturated Molecules, 
cara, Conference on Physics, London, October, 
1934, 


gives, on putting B=y=0, W=W,+5.605a. 
The difference W—W’ has just the value found 
by Pauling and Wheland. One can readily show 
that their method will always give the correct 
resonance energy. The fallacy in the argument 
to the contrary arises from the fact that it does 
not properly take into account the steric terms 
which appear in 

One must carefully distinguish between the 
theoretical resonance energy and the so-called 
“observed resonance energy’”’ used by Hiickel 
and by Pauling, Wheland and Sherman. The 
former is the difference between the theoretical 
energy obtained by treating the interaction of 
the orbits as an electron problem, and the 
energy calculated under the assumption of simple 
electron pair bonding. The latter is taken to be 
the difference between the observed energy and 
the energy computed from the empirical bond 
additivity law: the name “resonance energy” has 
been applied to this quantity under the sup- 
position that it is equal to the theoretical 
resonance energy. However, the theory as we 
have developed it indicates a serious breakdown 
of the bond additivity law as used for molecules 
such as benzene. For benzene, we find, under 
the assumption of simple electron pairing, W’ 
= —49.2 v, while the bond additivity law gives 
W’= —50.6 v. For 1,2 dihydrobenzol the cor- 
responding figures are —51.6 v and —56.3 v, 
and for napthalene —75.9 v and —79.2 v. This 
failure of the bond additivity law results from 
the very low value of the theoretical C=C 
bonding energy for such molecules, 3.77 v, 
compared to the empirical value, 5.56 v. The 
theoretical resonance energies for the three 
molecules are —3.0 v, —4.7 v, —5.5 v, respec- 
tively, while the ‘“‘observed resonance energies” 
are —1.6 v, —0.4 v,3! —3.2 v. The ‘‘observed 
resonance energies’ appear to have no theoretical 
significance whatsoever. It consequently seems 
very difficult to attach any meaning to calcula- 
tions, such as those of Pauling, Wheland and 
Sherman, which refer to these quantities. 


41 It should be noted that the Pauling-Wheland-Sherman 
procedure predicts a resonance energy, —2.6 v, in definite 
disagreement with this value. 
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An Experimental Investigation of the Effect of Pressure on Phase Equilibria of 
Sodium Tungstate and of Related Thermodynamic Properties 


R. W. GorANSON AND F. C. KRrAcEK, Geophysical Laboratory, Carnegie Institution of Washington 
(Received December 7, 1934) 


In order to clarify some uncertainties concerning the t=588.8+0.025(p—40) with a heat change of 105 joules/ 
physical properties of sodium tungstate the effect of pres- gram and volume change of 0.035 cm*/gram. The triple 
sure on phase equilibria to 1000 bars, the heats of inversion point lies at =588.8° and to 80 bars. The II—III 
and of melting, and the specific volume (density) at 30° _ inversion is very similar in behavior to the III] inversion; 
were determined. Sodium tungstate is trimorphous at here ¢=587.6+0.029 » on the heating curve. In order to 
atmospheric pressure, form I being stable from 695.5° represent specific volume asa function of temperature from 
(melting point) to 588.8°, form II stable from 588.8° to 20° to 1300° the data of Jaeger and of Austin were also 
just below 587.6°, and form III stable below this tempera- _utilized. The published density of sodium tungstate was 
ture. The I= II inversion is estimated to take place at found to be in error and therefore a redetermination was 
588.8° with a heat change of 14 joules/gram but no volume made. The value obtained by us, 5.13+0.01, is 20 percent 
change (independent of pressure). Because of the large higher than that given by F. W. Clarke (J. L. Davis). The 
amount of undercooling in the I-III inversion, the II—I low-temperature modification expands from 0.1959 cm*/ 
transition could be obtained metastably to 613 bars. gram at 30° to 0.203, cm*/gram at 587.6°. The salt inverts 
The III=I inversion takes place promptly and repro- at this temperature, attaining the volume 0.238, cm*/gram, 
ducibly on heating but not so on cooling, the temperature and finally, it again increases in volume at the melting 
change being dependent in large measure on the past history point, 695.5°, to 0.256; cm*/gram. 
of the charge. Heating curve values lead to the expression 


INTRODUCTION powdered sodium tungstate which, however, was 
not completely anhydrous.' The last traces of 
water were removed by melting the salt. 

Three crystalline forms of (anhydrous) sodium Boeke? states that the low-temperature modi- 
tungstate are known to exist but the exact fication of sodium tungstate is optically isotropic. 
relationship between these phases at atmospheric The material was found to have a weak bire- 
pressure is still somewhat in doubt, no work fringence, however, and examination of x-ray 
having been done at other pressures. The tem- powder photographs shows that the symmetry 
perature-pressure-phase equilibria relations of js not that of the cubic system.’ It has a mean 
these three solid phases and of the liquid phase refractive index (A=547my) of 1.678. 
were therefore determined to 1000 bars (10° 
dynes/cm2) pressure. In order to utilize these © TEMPERATURE-PRESSURE PHASE EQUILIBRIA 
data to obtain other thermodynamic properties The three crystalline modifications of sodium 
of sodium tungstate the heats of inversion and tungstate will be designated I, II and III, to 
of melting and the density at 30°C were re- denote the high, intermediate, and the low 
determined. From these data, together with temperature forms, respectively. The liquid 
Jaeger’s liquid density data and Austin’s thermal phase will be designated as L. 
expansion data from 20° to 500°C, the specific 


volume could be evaluated as a function of Previous work . 
temperature from 20° to 1300°C. The results recorded in the literature‘ show 


that on cooling two breaks are obtained in the 


Scope of the investigation 


Description of the material used 


4 ge 1 states that the dihydrate can be dehydrated at 100°. 
A = grade of sodium Cungnate wae fe It was found, however, that material dried at 125° retained 
crystallized either by cooling a hot saturated up to 0.5 percent of water. 
solution or by precipitating it by the addition of (1906) E. Boeke, Zeits. f. anorg. allgem. Chemie 50, 359 
95 percent ethyl alcohol to a nearly saturated ‘ woe indichend to Me. C. J. Ksanda of this Laboratory 
;_ for the x-ray photographs. 

— The resulting sodium tungstate di ‘ The literature to 1929 is given in Gmelin’s Handbuch 
ydrate was dried at 125°, yielding very finely der anorg. Chemie, Wolfram system No. 52, 212 (1933). 
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curve : (a) asmall heat effect in the range between 
590° and 576°, and (b) a large heat effect in the 
range between 583° and 560°. On heating there 
appears to be only a single arrest at about 588°. 
Van Klooster and Germs found that if cooling 
was arrested at 584° there was on subsequent 
reheating an arrest at 591°. Kracek® made further 
measurements and found that the small heat 
effect took place at 586° on cooling and that it 
could be reversed readily, provided the tem- 
perature did not fall low enough for the large 
heat effect to begin, taking place at 592° on 
reheating. The large heat effect was observed to 
take place with considerable undercooling, and 
reheating afterwards gave an apparently single 
arrest at 588°. This behavior is almost an exact 
parallel to that of potassium nitrate, and the 
modifications were accordingly numbered in the 
Smithsonian Tables to agree with Bridgman’s 
nomenclature for the latter substance,* assuming 
the intermediate form to be metastable at atmos- 
pheric pressure. Boeke observed with a heating 
microscope that the melt crystallizes in birefrin- 
gent dendritic crystals which transform at 590° to 
a still more birefringent form and then invert at 
‘584° to become optically isotropic. On reheating, 
the phenomena are said to occur in reverse 
order. Boeke’s observations would indicate that 
the intermediate form of sodium tungstate has a 
temperature range of stable existence at atmos- 
pheric pressure instead of coming in metastably 
as was supposed by Kracek. 

The results of this paper establish, among 
other things, the validity of Boeke’s observa- 
tions, and also show that, with proper technique, 
both arrests can be observed on heating as well 
as on cooling. 


Experimental procedure 


The pressure investigation was carried out 
in a wacer-jacketed steel bomb containing a 
platinum-rhodium resistance furnace. A detailed 
description of the bomb has been given in an 
earlier paper.’ Pressures were measured with a 
resistance gauge consisting of a hundred ohm 
coil of manganin wire immersed in the pressure 


5 Kracek, Smithsonian Physical Tables, 8th Edition, p. 
261 (1933). 

6 P, W. Bridgman, Proc. Am. Acad. 51, 579 (1916); F.C. 
Kracek, Chem. 34, 225 (1930). 

7 Roy W. Goranson, Am. J. Sci. 22, 481 (1931). 


fluid but kept at room temperature (25°). The 
coil is connected to the apparatus for measuring 
the resistance through a pressure-tight elec- 
trically insulated packing. The sensitivity in 
measuring changes in resistance (about 0.2 ohm 
per kilobar) was such that one scale division 
corresponds to about one bar. The scale for 
this electrical pressure gauge was calibrated by 
comparison with an absolute gauge of the free- 
piston type. The pressure medium used was 
nitrogen gas. 

Temperatures were measured with Pt-PtRh 
thermocouples. The leads are carried out of the 
bomb through pressure-tight electrically insu- 
lated packings in the bomb lid. The scale for 
these thermocouples was obtained by calibration 
at several fixed points and interpolated from 
the table of L. H. Adams.* 

The differential temperature, i.e., the differ- 
ence in temperature between the inert body, in 
this case silica glass, and the sodium tungstate 
charge, was measured with PtRh-AuPd-PtRh 
thermocouples. 

The inversion and melting temperatures were 
obtained from differential heating and cooling 
curves. 

Some of the measurements at atmospheric 
pressure were made in an ordinary platinum 
resistance furnace in connection with the calo- 
rimetric data. 


Experimental results 


The accuracy of the determination of the 
temperature-pressure-phase equilibrium curves 
for the III=—II and III=I inversions is limited 
by the rather wide passive regions within which 
the rate of transition is not measurable. The 
upper bound (heating curve values) of these 
regions is reproducible to within 0.5°, taking 
place promptly but not isothermally. The lower 
bound (cooling curve values), although prompt 
when initiated, is nevertheless variable and not 
readily reproducible because it depends in a 
large measure on the past history of the sodium 
tungstate charge. It was found that if the powder 
was first sintered or melted to a cake it would 
react more or less as a single crystal, i.e., the 
individual crystalline grains, by this procedure, 
seem to have been brought into contact intimate 


*‘L. H. Adams, Int. Crit. Tab. 1, 59 (1926). 
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Fic. 1. Typical differential heating (A, G) and cooling 
curves (B to F inclusive) through the two inversions. 


enough to initiate the reaction in their neighbors. 
Moreover, this procedure of first sintering or 
fusing to a cake increases the thermal con- 
ductivity and thereby lessens temperature gradi- 
ents in the charge; the unfused precipitate, 
dehydrated at 125°, is an extremely fine powder 
and a good thermal insulator, and such a charge 
would therefore have high temperature gradients 
when heating or cooling. Differential cooling 
curves yield a lower bound averaging about 6° 
below the upper bound, and undercooling of as 
much as 15° has been observed. By repeated 
reversals of the heating and cooling curves 
before the transition was complete the under- 
cooling could be apparently largely eliminated, 
but even with this technique the maximum 
value obtained at atmospheric pressure for the 
inversion temperature on the cooling curve was 
still 5° below the inversion temperature on the 
heating curve ; at 420 bars pressure this technique 
was able to bring the cooling curve value to 
within 1° of the heating curve value. The 
essential characteristics of the differential heating 
and cooling curves do not vary appreciably with 
pressure, but the curves do become more repro- 
ducible at high pressures, perhaps as a result of 
increased efficiency in heat transference. Such 
curves are reproduced in Fig. 1; they represent 
runs made at one bar pressure, but the shapes 
of the transition curves are typical of any of 
the pressure runs. The curves of Fig. 1 were all 
obtained with the same sample; cooling curves 
B to F inclusive demonstrate the variability in 


TABLE I. Zemperature-pressure relations of the III=II 
and III=1 inversions. ‘‘H’’ denotes the beginning of in- 
version on heating and ‘‘C”’ the beginning on cooling. The 
asterisk value is calculated (as shown later in the text). 


Pressure 
in bars 


Pressure 


Temperature 
in bars H Cc 


Temperature 
H Cc 


598.5 
603.9 
606.0 
605.5 
611.3 
613.0 
614.0 


587.6* 
587.5 
538.8 
589.3 
592.4 
594.6 
595.5 
599.2 
599.2 


599.4 
598.8 


605.2 
606.4 


581.1 
582.6 
582.7 
583.6 


592.5 
597.7 


values, indicated by “‘e,”’ obtained for the IIT>III 
inversion. 

The experimental results for the III=II and 
the III=I inversions are tabulated in Table I. 

The II-I inversion, when it follows immedi- 
ately upon the III-II inversion, is not readily 
detectable on the differential heating curve be- 
cause it takes place only one and a half degrees 
above the beginning of the III-II inversion at 
atmospheric pressure and with only about one- 
seventh the heat effect of the latier. With careful 
technique the break can be obtained (Fig. 1, 
points “‘b’”’ on curves A and G) and this break 
can be easily accentuated by reversing the 
cooling curve before the I[-III transition has 
been completed. This inversion is prompt but 
not isothermal; superheating and undercooling 
of about one-half degree each have been ob- 
served for this transition. Because of the large 
amount of undercooling in the I>III transition 
it was found that ordinarily the transitions would 
take place metastably in the order I>II-III 
up to 613 bars pressure. It was also possible to 
reverse the cooling curve before the II-III 
transition occurred and so to obtain the meta- 
stable III transition up to 300 bars pressure. 
Typical records of such reversal curves are 
presented in Fig. 2 and were made in the order 
as numbered. The metastable prolongation of 
the I=II pressure-inversion curve was therefore 
determined to 613 bars pressure. The experi- 
mental results for the II=—I inversion are 
collected in Table II. 

The temperature-pressure results given in 
Table I for the III-I (heating curve values) 
transition can be represented, within the accu- 
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Fic. 2. Typical differential heating and cooling curves of 
the I=II inversion made in the order as numbered without 
going through the II=III inversion. 


racy of the experiments, by the slope dt/dp 
= 0.025 or, if no account be taken of stable and 
metastable states and extrapolated to one bar, 
by the equation 


= 587.8 +0.025p, 


where ¢ denotes temperature in degrees C, and 
p pressure in bars. 

For reasons already given, the cooling curve 
values (I-III) were not used in the calculation, 
especially as in a number of the experiments the 
cooling curve values represent the metastable 
II-III rather than the stable I-III transition. 


TABLE II. Temperature-pressure relations of the II=—I 
inversion. “‘H’’ and “‘C” denote the beginning of inversion 
on the differential heating and cooling curves, respectively. 


Pressure 


Temperature 
in bars Cc 


The temperature-pressure curve for the II=I 
inversion, as can be seen by inspection of Table 
II, is independent of pressure, within the experi- 
mental error; that is, 


tu—1= 588.8+0.000p. 


The location of the triple point can be evaluated 
only within limits from the data available. 
The intersection of the III-I and the II=I 
curves yields p= 40 bars at 588.8°. This must be 
regarded as the lower pressure limit. The true 
location of this point depends, obviously, upon 
the location of the true equilibrium temperature- 
pressure curve for the III=I inversion, and, 
hence, lies at a somewhat higher pressure. Now, 
since at 420 bars the III=I band was found to 
be only one degree wide, and if we suppose that 
this width would persist to lower pressures if 
the I-III transition could take place without 
the metastable intervention of the sequence 
I>II-III, it may be estimated from the data 
that the upper limiting value of pressure for the 
triple point is 80 bars at 588.8°. Hence, 


Na.WO,(I, II, III): = 588.8° 
p=40 to 80 bars. 
Melting curve 


The temperature of melting is taken as the end 
of melting and the temperature of freezing as 
the beginning of crystallization or as the maxi- 
mum where undercooling occurs. The experi- 
mental values are presented in Table III. 

TABLE III. Pressure melting of Na2W0O,. ‘‘H” denotes 
the end of melting and “C” the beginning of crystalliza- 


tion. The lower four values denote a single sequence and are 
used in calculating the slope of the curve. 


Temperature 
H Cc Mean 


695.5 


695.0 695.0 
705.4 705.5 
712.0 711.7 
716.6 716.1 


Pressure 


The pressure-melting relations, where tem- 
perature represents the mean of the melting and 
freezing temperatures, was represented as a least 
square quadratic passed through the points and 
is given by the following equation: 
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Fic. 3. The temperature-pressure phase diagram of sodium 
tungstate. ~ 


where ¢ is in °C and p in bars. 

The temperature-pressure phase diagram for 
sodium tungstate is given in Fig. 3. Circles with 
a darkened half-segment represent the I=II 
inversion where ¢ is the mean of heating and 
cooling curve values except above 300 bars where 
only cooling curve values could be obtained. 
The dotted portion of the curve denotes the 
metastable prolongation of the I=II inversion. 
Open circles represent the IIIoII and III-I 
inversion (heating) ; plusses and crosses represent 
the values obtained from cooling curves, the 
crosses denoting the I-III inversion and plusses 
the II-III inversion. Dark circles represent the 
means of the melting and freezing temperatures. 


CALORIMETRIC DATA 


The heats of melting and of the two in- 
versions were measured by the differential heat- 
ing and cooling curve method described in 
earlier papers. The values are based on time- 
differential temperature areas for which calibra- 
tion of the thermal equivalent of the apparatus 
was made to 500°. No basic calorimetric data 
suitable for use with our apparatus (with the 
possible exception of the value of the heat of 
freezing of NaCl, which, however, volatilizes so 
copiously at its melting point as to give very 
erratic results) exist for the temperature range 
above 500°, and hence it was necessary to base 
our determinations on an extrapolation of the 


°F. C. Kracek, J. Phys. Chem. 34, 225 (1930); R. W. 
Goranson and F. C. Kracek, ibid. 36, 913 (1932). 


thermal equivalent curve beyond this tem- 
perature. 
The results obtained were as follows: 


Liquid=Na2WO,I, joules/gram 
Na2WO,I=Na.WO,II, 14 “ 
Na2WO,II=Na.WO,III, 105 “ 


Hare’ has measured the heat of the combined 
transitions (I-III) by a rather different radiation 
method which depends upon the absolute evalua- 
tion of the radiation constants of the apparatus 
for each experiment. He obtained 136 joules/ 
gram. It should be noted in this connection, 
however, that unfused, ignited powder obtained 
by dehydration of the dihydrate gave for the 
inversions the values 18 and 123 joules/gram, 
published by Kracek in the Smithsonian Tables. 
It was found subsequently that the same filling 
of the apparatus always gave a larger value 
with the unfused powder than after the salt had 
been fused. The discrepancy probably is due to 
the low thermal conductivity of the original 
powder. The values given in the tabulation above 
are for material after fusion. The experimental 
error is estimated to be +5 percent. 


SPECIFIC VOLUME DATA 


Changes of specific volume at the transitions and 
on melting 
From the calorimetric data and the slopes of 
the transition curves the changes in the specific 
volume at the transitions can be calculated. 
The results of such calculations are given in 
Table IV. 


TABLE IV. Thermodynamic properties of the transitions at 
atmospheric pressure. 


dp/dt AQ Av 
Transition deg/bar joules/gram cm*/gram 


d 0.029 105 0.035 
Il >I 588. i+ 0.00 
I —L 81 0.018, 


Temp. at 


For the III-I transition, which is not stable 
at atmospheric pressure, dt/dp=0.025, AQ=119 
joules/gram, and Av=0.035,) cm'’/gram. The 
temperature of the III-II inversion at atmos- 
pheric pressure is calculated to be 587.6°, which 
agrees with the experimentally determined value. 


10 A. Hare, Phil. Mag. 48, 418 (1924). 
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Specific volume of the liquid 


Jaeger™ has measured the density of liquid 
sodium tungstate at three temperatures. His 
values are: 3.685 at 917°, 3.502 at 1128°, and 
3.356 at 1330°. A recalculation of these points 
was made in terms of specific volume yielding 
the equation 


v=0.2714+0.703(¢—917)10~ 
— 1.423(t—917)710- 


where v denotes the specific volume of the liquid 
which, from this equation, is 0.257 at the 
melting point. 


Specific volume of the crystals at room tem- 
perature 


The only published value of the density of 
anhydrous sodium tungstate at room tempera- 
ture is that by Clarke and Davis,” who give 
4.1833 at 18.5° and 4.1743 at 20.5°. From these 
data and those of Jaeger the change in specific 
volume from solid at 20° to liquid at 695° is 
calculated to be 0.017. This result is inconsistent 
with the data of Table IV from which the 
specific volume change due to the transitions and 
melting alone is calculated to be 0.0532. Accord- 
ingly, the density at room temperature was re- 
determined. 

The density results on the very fine- 
grained powder resulting from dehydration of 
Na2WO,.2H.0 clustered about the value 4.82, 
with deviations of about 0.01, while those on the 

fused material were in the neighborhood of 5.10, 
and depended on the grain size, the coarsest size 
used (10-24 mesh) yielding the highest value. 
This variability resulted from an imperfect re- 
moval of the air trapped in the sample. Con- 
sequently the ordinary technique was modified 
by squeezing toluene (the displacement fluid) 
into the pores of the grains under a pressure of 
6000 bars for about one-half to one hour; this 
additional treatment increased the measured 
density by about two parts in 500. Accordingly, 
the density of sodium tungstate is taken to be 
5.134+0.01 at 30°. The corresponding specific 
volume is 0.195+0.0004 cm/gram. 


uF, M. Jaeger, Zeits. f. anorg. Chemie 101, 183 (1917). 
377) W. Clarke (J. L. Davis), Am. J. Sci. (3), 14, 283 
1877). 
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Fic. 4. The specific volume of sodium tungstate as a func- 
tion of temperature. 


Thermal expansion of sodium tungstate 


The linear thermal expansion of sodium tung- 
state has been determined by J. B. Austin™ 
from 20° to the inversion temperature. We have 
calculated a least square quadratic for the 
volume expansion of modification III from his 
data to 500° on the assumption that his samples 
represent aggregates of randomly oriented crys- 
talline grains. The equation is 


V20[_53.194(¢— 20) 10-8 


which, if the measured value of the specific 
volume at 30° be inserted, becomes 


v=0.1947-+ 10.4(¢—20)10-*+ 20)7107°, 
v denoting specific volume. 


Specific volume as a function of temperature 


By combining the various data we can finally 
express the specific volume of.sodium tungstate 
as a function of temperature from 20° to 1300° 
at atmospheric pressure. In this temperature 
interval there are two discontinuities, the III II 
transition at 587.6° with Av=0.0359 and the 
I=L transition at 695.5° with Av=0.0182. These 
results are represented graphically in Fig. 4. 


‘8 We are indebted to Dr. Austin for allowing us to use his 
results before publication, 
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The Compressibility of and an Equation of State for Gaseous Ethane 
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The compressibility of gaseous ethane (C2H¢) has been determined from 0.5 to 5.0 moles per 
liter and from 25° to 250°C with an estimated accuracy of 0.2 percent or better. Above 250°C 
ethane decomposes at an appreciable rate when confined in steel. An equation of state is fitted to 
the data. The average deviation of the calculated from the observed pressure is 0.13 atmosphere 
or 0.21 percent. The equation of state extrapolates to the critical pressure (at 7 moles/liter and 
32.1°C) with an error of 10 percent and calculates the weight of a normal liter to 0.2 percent. 
The vapor pressures of liquid ethane at 0° and 25°C are given. 


HE compressibility of the hydrocarbons is 

of practical as well as theoretical interest. 

The only hydrocarbon for which extensive mod- 

ern compressibility data are available is methane, 

which has been investigated by Keyes and Burks.! 

We have studied the compressibility of gaseous 

ethane from 25° to 250°C and from 0.5 to 5.0 
moles per liter. 

The apparatus and the method of procedure 
have been described in detail elsewhere.? A 
weighed quantity of ethane (approximately 3 
grams or 0.1 mole) was confined’ in a chrome- 
vanadium steel bomb (of 200 cc volume) which 
was placed in a thermostat regulated to +0.001°C, 
the temperature being measured with a platinum 
resistance thermometer. The density of the gas 
was varied and measured by the introduction or 
removal of known volumes of mercury from a 
calibrated mercury injector thermostated at 30°C. 


Pressures were measured oa a dead-weight piston 


gauge. The effect of the temperature of the bomb 
and of pressure on the volume measuring system 
was determined experimentally in a blank run 
in which there was no ethane in the bomb. 

The pressures were measured along isometrics 
as well as along isotherms. At each temperature 
the method of procedure was to work from the 
smallest to the greatest gas density and then to 
check the pressure at the smallest density. At 
275° it was found that the pressure had decreased 
during the course of the measurements, and the 


* Contribution No. 335. 
'F. G. Keyes and H. G. Burks, J. Am. Chem. Soc. 49, 
1403 (1927). 
1034) Beattie, Proc. Am. Acad. Arts and Sci. 69, 389 
{ 


pressure corresponding to a given density de- 
creased slowly but steadily with time. Remeas- 
urement of the vapor pressure at 25°C showed 
an increase in pressure over the previously deter- 
mined value and also an increase in pressure with 
decreasing vapor space in the bomb, but no ap- 
preciable change with time. This indicates that 
at about 275° in steel, ethane cracks to give, as 
one product, hydrogen which rather rapidly dif- 
fuses through the walls of the bomb at this tem- 
perature. When the bomb was opened, a thin film 
of an oily tar was found on the inner walls. 

Temperatures are given on the international 
platinum resistance thermometer scale. The 
strain-free mica-cross platinum resistance ther- 
mometer was calibrated at the ice-point (0°), the 
steam point (100°) and the sulfur boiling point 
(444.6°) at the beginning and at the end of the 
measurements. Special hypsometers, which per- 
mit very exact calibration at a constant pressure 
of one normal atmosphere, were used. The maxi- 
mum difference between temperatures computed 
from the two calibrations was 0.001°, this maxi- 
mum occurring at 250°C. Each day the ice-point 
of the thermometer was determined and correc- 
tion made for the variation, which amounted at 
most to 0.002°C during the course of the investi- 
gation. The precision of the determination of 
temperature with one platinum resistance ther- 
mometer is about 0.002°C, but the uncertainty 
of the readings of a single thermometer relative 
to the international platinum resistance scale 
increases above 100°C to about 0.01° to 0.02°C 
at 250°C. 

Pressures are given in normal atmospheres. 
The, pressure gauge was calibrated at the vapor 
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TABLE I. Values of the constants of the Beattie-Bridgeman equation of state for gaseous ethane (C2H.). 


A=A,(1—a/V) B=B,(1—6/V) e=c/VT? 


Ao Bo 


Units: normal atmospheres, liters per mole, °K (T°K =i°C +273.13). 
0.05861 0.09400 0.01915 


Amagat units: normal atmospheres, V =1 at 0°C and 1 atm., °K 


3.69658 X 10-3 11.9320 X 10-8 2.6402 X 10-8 4.2344 X 10-3 0.8627 X 10-3 40.543 X 108 


Berlin units: meters of mercury, V =1 at O°C and 1 meter of mercury, °K 


15.7971 X 10-3 3.4847 X 10-3 5.5888 X 10-3 1.1386 X 10-3 53.510 X 108 


3.70799 X 10-3 


TABLE II. Comparison of the pressures calculated from the equation of state with the observed pressures for gaseous ethane. 
For each temperature the first line gives the observed pressure and the second line gives the observed minus the calculated 
pressure. The equation used for the calculated pressures is given in Table I. 


Density, ' eel 0.5 1.0 1.5 2.0 


Temp., 


2.5 3.5 
Pressures, normal atmospheres 


32.84 
—0.25 


38.41 
—0.11 


43.80 
—0.02 


49.03 
—0.01 


54.25 
+0.06 


59.31 
+0.02 


27.34 
—0.19 


31.24 
—0.11 


35.08 
—0.01 


38.79 
+0.00 


42.49 
+0.03 


11.11 
—0.05 


12.24 
—0.02 


13.34 
—0.02 


14.43 
—0.03 


15.52 
—0.03 


16.60 
—0.03 


17.67 
—0.05 


18.71 
—0.09 


19.80 
—0.09 


20.89 
—0.08 


0.049 
0.293 


20.14 
—0.14 


22.58 
—0.07 


24.95 
—0.05 


27.28 
—0.05 


29.64 
+0.00 


31.89 
—0.05 


34.17 
—0.06 


36.44 
—0.07 


38.64 
—0.15 


40.87 
—0.19 


0.083 
0.276 


25 obs. 


obs.-calc. 
50 obs. 
obs.-calc. 


obs. 
obs.-calc. 


obs. 
obs.-calc. 


obs. 
obs.-calc. 


obs. 
obs.-calc. 


46.08 
—0.01 


49.66 
—0.05 


53.25 
—0.05 


56.73 
—0.15 


60.23 
—0.22 


0.082 
0.198 


obs. 


64.36 
obs.-calc. +0.02 


69.38 
4-0.01 


74.27 
—0.10 


79.16 


obs. 
obs.-calc. 


obs. 
obs.-calc. 


obs. 
obs.-calc. 


Average deviation (atm.) 
Average % deviation 


—0.18 


0.078 
0.166 


36.88 39.66 41.35 
—0.25 —0.16 +0.00 


44.27 49.02 52.87 
—0.08 +0.00 +0.15 


51.37 57.98 63.78 
+0.00 +0.05 +0.11 


58.29 66.68 74.36 
+0.04 +0.06 +0.05 


65.13 75.26 84.81 
+0.11 +0.12 +0.08 


71.78 83.61 94,99 
+0.07 +0.07 +0.02 


78.40 91.93 105.13 
+0.07 +0.09 +0.06 


100.19 115.20 
+0.13 +0.14 


108.24 125.02 
+0.02 +0.06 


116.30 134.84 
—0.03 +0.05 


0.073 0.072 
0.110 0.090 


55.97 
+0.35 


68.89 
+0.10 


81.50" 
—0.04 


93.91 
—0.08 


106.06 
—0.13 


118.16 
—0.05 


130.18 
+0.09 


141.90 
+0.06 


153.59 
+0.10 


0.111 
0.140 


58.49 
+0.57 


73.52 
+0.05 


88.25 
—0.24 


102.73 
—0.37 


117.00 
—0.40 


131.20 
—0.26 


145.30 
—0.02 


159.08 
+0.05 


172.87 
+0.26 


0.247 
0.267 


60.56 
+0.75 


77.79 
—0.12 


94.73 
—0.61 


111.46 
—0.80 


127.98 
—0.80 


144.41 
—0.58 


160.77 
—0.18 


176.73 
+0.01 


192.77 
+0.44 


0.477 
0.459 


0.083 
0.151 


Total average deviation (atm.), 0.131; total average % deviation, 0.213 


pressure of liquid carbon dioxide at 0°C for which 
the value 34.401 normal atmospheres’ is used. 
The gauge was recalibrated after the completion 
of the measurements and showed the normal in- 
crease in gauge constant (normal mm of mercury 
per gram load) with time (0.01 percent per year), 
for which correction was made. The uncertainty 
in the determination of pressure is less than 0.05 
percent. 


30. C. Bridgeman, J. Am. Chem. Soc. 49, 1174 (1927). 


Gas volumes at each temperature are expressed 
in liters per mole and are referred to the specific 
volume of liquid mercury at the temperature in 
question and 1 atmosphere. The density of mer- 
cury at 0°C was taken to be 13.5955 grams per 
milliliter and the effect of temperature on the 
volume of mercury at 1 atmosphere was obtained 
from the equation given in International Critical 
Tables.‘ At 0.5 mole per liter gas density, the gas 


4 International Critical Tables, Vol. 11, p. 457. McGraw- 
Hill Book Co. Inc., New York, 1927. 
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volume in the bomb was about 200 cc; at 5.0 
moles per liter it was 20 cc. The mercury com- 
pressor can be calibrated to +0.004 cc at one 
atmosphere pressure and the uncertainties of the 
corrections applied for the effect of pressure and 
temperature on the volume measuring system are 
about 0.02 cc. Yet the errors in the measurement 
of volume are larger than these figures indicate, 
being about 0.1 percent. At each measurement 
we waited from 1 to 2 hours for equilibrium and 
usually no change in pressure was observed after 
the first half hour. 

The mass of ethane was determined by weigh- 
ing in a stainless steel bomb against an almost 
identical counterpoise. Proper air-buoyancy cor- 
rections were made. The accuracy of the deter- 
mination of mass was 0.01 percent. 

The final tabulated pressures are believed to 
have an accuracy of 0.2 percent or better. 

The ethane was obtained from the Linde Air 
Products Research Laboratory, and guaranteed 
to be at least 99.9 percent pure. In order to 
remove any non-condensible gas which might 
have been present, it was distilled by us three 
times from one bulb to another by use of liquid 
air, the first and last portions being discarded. 
After each distillation the ethane was frozen, and 
pumped by means of a mercury diffusion pump. 
It was then loaded into the weighing bottle, 
weighed, and transferred to the compressibility 
bomb. That the ethane was of extraordinary 
purity was indicated by the fact that at 0° the 
vapor pressure of the liquid remained constant 
within our experimental error (0.001 atmosphere 
for relative measurements at constant tempera- 
ture) while the vapor volume was decreased from 
20 to 5 cc; and at 25° the vapor pressuré remained 
constant within the experimental error while the 
vapor volume was decreased from 40 to 5 cc. 

The constants of the Beattie-Bridgeman equa- 
tion of state’ were determined, from the data. 
They are given in Table I. In Table II are given 
the observed pressures and the deviation of the 
pressures calculated from the equation of state 
from the observed values. The average deviation 
is 0.13 atmosphere or 0.21 percent. Although the 
equation of state gives quite a satisfactory rep- 
resentation of our data, there is one trend in the 


°J. A. Beattie and O. C. Bridgeman, Proc. Am. Acad. 
Arts and Sci. 63, 229 (1928). 


data that is not reproduced by the equation. The 
equation requires that the curvature of the iso- 
metrics at any temperature increase (numeri- 
cally) with density. The curvature of the experi- 
mental isometrics first increase with increasing 
density (to about 3.5 moles per liter) and then 
begin to decrease. In other words the constant c 
of the equation of state should slowly decrease 
with increasing density in order that the equation 
represent this trend in the ethane data. 


VAPOR PRESSURES 


In the course of our investigation we measured 
the vapor pressure of liquid ethane at two tem- 
peratures. The results are given in Table III. 


TABLE III. Vapor pressure of liquid ethane. 


Temp.°C Vapor pressure, normal 
atmospheres 
LC.TS 


Our values 
23.56 


41.5 


23.57 
41.37 


WEIGHT OF A NORMAL LITER 
In Table IV are given the value of the weight 


TABLE IV. Weight of a liter of ethane at 0°C and 1 normal 
atmosphere. 


1.3535 grams 
1.3566 grams 


Calculated from the equation of state 


of a liter of ethane at 0°C and 1 normal atmos- 
phere calculated from the equation of state, and 
the observed value. This is not a satisfactory 
agreement; the difference may be due to the 
equation which in general gives a low value for 
the weight of a normal liter for compressible 
gases (for the gases so far studied the deviation 
is smaller than in the case of ethane); or to the 
purity of the ethane used in the measured weights 
of a normal liter, the suspected impurity (COz) 
giving a high value. 


CALCULATION OF THE CRITICAL PRESSURE 


The critical data for ethane’ are given as 48.8 
atmospheres, 0.21 gram per cc (7 moles/liter), 
and 32.1°C, the critical density being uncertain. 


6 International Critical Tables, Vol. III, p. 230. 
7 International Critical Tables, Vol. III, p. 3. 
8 International Critical Tables, Vol. III, p. 248. 
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If we substitute the critical temperature and 
density into the equation of state, we calculate 
44.3 atmospheres for the critical pressure, 10 per- 
cent lower than the observed value. This should 
not be considered a bad disagreement since it 
involves an extrapolation of 40 percent from our 
highest measured density. It should be pointed 
out that temperature extrapolations are in gen- 
eral fairly safe but density extrapolations are 
dangerous. This effect is one that is true of equa- 


tions of state in general and it is also true of the 
experimental data. Since the isometrics are so 


nearly linear, the experimental data can safely 


be extrapolated over quite a temperature range; 
but in the region near the critical, the isotherms, 
no matter what variables are plotted, are too 
curved to permit extrapolation. That is, from the 
data in the range 0.5 to 5.0 moles per liter, we 
cannot predict the data in the region of 7 moles 
per liter. 
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The Critical Constants of Propane 


James A. BEATTIE, NOLAND POFFENBERGER AND CANFIELD HaDLock, Research Laboratory of Physical Chemistry Massa- 
chusetts Institute of Technology* 


(Received November 17, 1934) 


The critical constants of propane (C;Hs) are: t-=96.81+0.01°C, p.=42.01+0.02 normal at- 
mospheres, v, = 0.195 liter per mole (4.43 cc per gram), d, = 5.13 moles per liter (0.226 gram per cc). 
The uncertainty in the critical volume and density is 1 percent. 


N the course of our measurements of the com- 
pressibility of gaseous propane,! we investi- 

gated a number of isotherms in the critical region 
and located the critical point. The method used 
was the same as that employed in the measure- 
ments of the compressibility of propane and of 
ethane,” and has been described elsewhere.* 

Pressures are given in normal atmospheres; 
temperatures in degrees centigrade on the inter- 
national platinum resistance scale, and volumes 
in liters per mole. In the region of the critical 
point of propane the uncertainties in the meas- 
urement of pressure do not exceed 0.03 percent; 
those in the location of temperature do not ex- 
ceed 0.005°C; and those in the measurement of 
volume do not exceed 0.05 percent. 

In the determination of the critical data, the 
purity of the substance is of primary importance. 


* Contribution No. 339. 
1J. A. Beattie, N. Poffenberger and C. Hadlock, to be 
published shortly. 
2J. A. Beattie, C. Hadlock and N. Poffenberger, J. 
Chem. Phys. 3, 93 (1935). 
( te Beattie, Proc. Am. Acad. Arts and Sci. 69, 389 
1934). 


The propane was obtained from the Research 
Laboratories of the Linde Air Products Company 
at Buffalo and was stated by Doctor L. I. Dana 
to be 99.9 percent pure or better. That the pro- 
pane was of extreme purity is indicated by the 
measurements of vapor pressure at 50° and 75°C 
given in Table I. It will be noticed that at 50°C 


TABLE I. The effect of varying the vapor volume on the vapor 
pressure of propane at 50° and 75°C. 


50°C 75°C 
Vapor Vapor Vapor Vapor 
volume pressure volume pressure 
cc normal atm. normal atm. 


16.897 28.112 
16.899 d 28.113 
16.898 28.113 
16.898 28.113 
16.898 28.113 
16.897 28.113 
16.897 28.112 
16.898 

16.898 

16.901 


the vapor pressure remained constant within the 
accuracy of our measurements (+0.001 atmos- 
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CRITICAL CONSTANTS OF PROPANE 


TABLE II. Isotherms of propane (C;Hs) in the critical region. 


Temp., °C 96.80 
Density Volume 


moles/liter liters/mole 


96.81 


96.82 96.90 97.00 97.10 97.20 


Pressure, normal atmospheres 


42.0115 


42.0120 
42.0130* 
42.0135 
42.0155 
42.0190 


42.0215 
42.0290 


41.9300 


41.9290 


41.9305 
41.9310 
41.9355 
41.9425 
41.9865 


42.0120 
42.0200 
42.0290 
42.0785 


41.3965 
41.3970 


42.0360 


42.0675 
42.0730 
42.0765 


42.0995 


42.1360 
42.1445 
42.1515 
42.1525 
42.1560 


42.1580 
42.1595 
42.1620 
42.1655 


42.1725 


42.1695 


42.2060 
42.2150 
42.2245 
42.2280 


42.2320 


42.2370 
42.2400 


42.2335 


42.2775 
42.2895 
42.3010 
42.3040 


42.3090 


42.3150 
42.3190 


42.0145 
42.0175 
42.0190 
42.0195 


42.0205 
42.0220 
42.0235 
42.0240 
42.0260 


42.0295 
42.0365 


42.0800 
42.0855 
42.0920 


42.2535 42.3345 


42.1105 
42.1675 


42.1960 
42.2575 


42.2810 42.3640 


* Critical Point 


phere) when the vapor volume was decreased 
from 93 to 18 cc; and that the vapor pressure 
increased 0.002 to 0.003 atmosphere when the 
vapor volume was further decreased to 0.15 cc. 
That is, during a 600-fold decrease in the vapor 
volume the increase in vapor pressure was prac- 
tically within our experimental error. At 75°C 
the vapor pressure remained constant while the 
vapor volume was decreased from 50 to 7 cc. 

In Table II are given the compressibility data 
on propane in the critical region, the isotherms 
between 96.7 and 96.9 being plotted in Fig. 1. 
The pressures are given to 0.0005 atmosphere 
since relative values are consistent to about 0.001 
atmosphere. In the two phase region the densities 
listed are the average densities of the whole 
system. 

From Fig. 1 it is evident that the critical tem- 
perature is 96.81-+0.01°C and the critical pres- 
sure 42.01++0.02 normal atmospheres. The deter- 
mination of the critical density cannot be made 
with so great an accuracy as those of temperature 
and pressure. The top of the ‘‘steam dome’”’ and 
the inflection point with horizontal tangent of 
the 96.81°C isotherm agree quite well for a criti- 


Pressure (normal atmospheres) 
& 


96.70 


0.18 0.22 0.24 


0.20 
Volume (liters per mole) 


Fic. 1. Isotherms of propane in the critical region. The 
radius of each circle is 0.002 atmosphere. 


cal volume of 0.195 liter per mole or a critical 
density of 5.13 moles per liter. The critical vol- 
ume seems to have been located to +1 percent. 
International Critical Tables* gives t-=95.6°C, 
p-=43 atmospheres. 


4 International Critical Tables, Vol. III, p. 248. McGraw- 
Hill Book Co. Inc., New York, 1928. 
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4.201 0.2380 41.9030 41.9660 
4.341 "2304 41.9175 
4.490 12227 41.3975 41.9280 41.9940 
4.651 "2150 41.9300 41.9990 42.0075 
4.823 ‘2073 41.9295 42.0020 42.0095 
4.895 "2043 42.0030 42.0110 
: 
5.293 "1889 42.0060 
5.380 "1859 42.0075 
= 
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The Relation Between Vector and Bond Eigenfunction Methods for Spin Degeneracy 


RICHARD S. BEAR! AND HENRY EyrinG, Department of Chemistry, Princeton University 


(Received October 29, 1934) 


A method previously developed for the calculation of 
matrix components for problems in which spin degeneracy 
only is important is shown to be equivalent to an operator 
method derivable from the vector approach of Dirac. 
Complete rules for the simple construction of matrices for 
all possible cases which can arise under these conditions of 
spin degeneracy are given. An especially simple way of 


determining secular equations for states of the identical 
representation is described, and the method is applied to 
the calculation of the energy of an infinite linear extension 
of sodium atoms. An inspection method of resolving a 
crossed six bond eigenfunction in terms of the uncrossed 
set is given. 


INTRODUCTION 


ECENTLY Van Vleck? and Serber*® have 

extended the vector method of Dirac,* 
applying it to the calculation of matrix coeffi- 
cients of exchange integrals. Van Vleck was 
chiefly concerned in demonstrating the value of 
the vector approach in single-atom systems, but 
indicated its utility in molecular problems. A 
more generally useful approach to the latter, 
however, is that employing bond eigenfunctions.* 
This paper considers relations between the vector 
and bond eigenfunction methods for configura- 
tions of orbits in which spin degeneracy only is 
important. In addition certain improvements in 
the mechanics of calculation with bond eigen- 
functions are described and applied to a con- 
sideration of the energy of an infinite linear 
extension of sodium atoms. 


THE VECTOR AND BonD EIGENFUNCTION 
METHODS 


In considering a system of WN electrons dis- 
tributed among an equal number of orbits 
a, b, «++ in the many centered field of a molecule 
or group of atoms one proceeds as follows. 
A complete set of antisymmetric eigenfunctions 
is constructed for the configuration in a manner 


1 National Research Fellow in Chemistry at Princeton. 
2 Van Vleck, Phys. Rev. 45, 405 (1934). 
3 Serber, Phys. Rev. 45, 461 (1934). 
4 Dirac, Proc. Roy. Soc. A123, 714 (1929). 
5 (a) Slater, Phys. Rev. 38, 1109 (1931); 
( a Kimball and Eyring, J. Am. Chem. Soc. 54, 3876 
1932); 
(c) Eyring and Kimball, J. Chem. Phys. 1, 239 (1933); 
(d) Pauling J. Chem. Phys. 1, 280 (1933). 
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elaborated by Slater.6 One of these can be 
written 


1 
2X (—1)?P(aa(1), 68(2), 


(N!)} P 
in which aa(1) expresses the orbit-spin eigen- 
function occupied by electron 1, b8(1) that of 
electron 2, and so on. The sum of the P’s per- 
mutes the electrons in NV! ways among the orbit- 
spin combination orbitals, and (—1)? is positive 
for even, negative for odd permutations. There 
are 2" eigenfunctions of the type ¢; obtained by 
taking all possible combinations of the two spin 
eigenfunctions a or 6 for each electron. The ¢’s 
formed in this manner are all independent, and 
though in general not orthogonal, they form a 
convenient set in terms of which we suppose any 
state arising when perturbations are introduced 
can be expanded. Because of the method of 
formation, the ¢’s are already eigenfunctions of 
the operator S, for spin momentum quantized 
along an optional z axis. A convenient way of 
expressing all y’s with a particular eigenvalue 
for S, is the following: 


(—1)?P(aa(1), b8(2), ++) 


Yr 


(2) 


Here ¢; is a particular one of the antisymmetric 
eigenfunctions with the required eigenvalue of S:. 
The operator Q, permutes the spins of g; with 
respect to the orbits in such a way as to give 
rise to a new function ¢,. 


6 Slater, Phys. Rev. 34, 1293 (1929). 
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Frequently it is useful to consider another set 
of independent functions instead of the g’s. 
These will be linear combinations of y’s which 
have the same eigenvalue of S,: 


Lengi= LenQien (3) 


where c,; is the coefficient of g; in the expansion 
of y,. 

If one now applies the perturbations of electron 
interactions, H, new states arise which can be 
expanded in terms of the set of y’s. Perturbation 
theory then shows that the matrix of T—W, for 
the set of y’s can be used as a secular determinant 


whose solutions for W, the energy, when equated . 


to zero will be the required energy levels of the 
perturbed system. The only matrix components 
which will be different from zero are those 
between y's with the same eigenvalues for S,. 
Hence we can write nonvanishing components as 
follows: 


[ver 
(4) 


Consider the integral 


(5) 


Ort, = 


The integrals of the type /(Qu¢1*)Q:(Q.¢1)dr 
are unity if Q.g=Q:Q.¢, otherwise they are 
zero, providing we neglect non-orthogonality of 
the orbitals as we shall do here. 

The integral, Q,:, :, defined in (5) is the 
coefficient of an expression (Kgi—WaAg;) in 
the expansion of (4), where Kg; is an exchange 
integral between spin-orbit products which differ 
in permutation of their spins by the permutation 
Qi. Ag; is the corresponding integral with unity 
instead of H between the spin-orbit products. 
Because we neglect non-orthogonality of our 
orbitals unless 7=1, ie., unless Qi is 
the identical permutation. The particular iden- 
tical permutation Q, gives rise to the quantity 
(J—W), where J is the familiar sum of coulombic 
integrals. The usual designation of exchange 
integrals, however, is with respect to a permuta- 
tion, Pi, of the orbits in relation to the spins, 


such that for the corresponding integrals we have 
Prt, i= i, the positive or negative sign 
occurring when the permutation is even or odd, 
respectively. We have, therefore, the expansion 
of (4): 


f —Wdr=(J—W)P rt, 1 
+2 (6) 


The expression (6) is similar to Eq. (14) of 
Serber.’ We can also write (6) in, the form 


i<j 


where the subscript.ij indicates an exchange 
between orbits 7 and j. The quantities J,,; and 
P,., 4; are the matrix components of unity and P;; 
between the states y, and y;. K,;; is the exchange 
integral between the orbitals z and j. 

Bond eigenfunctions have the important prop- 
erty of being eigenfunctions of S*, the operator 
for the square of the total spin angular mo- 
mentum, as well as for S,. They are formed 
from the antisymmetric functions (¢’s) which are 
eigenfunctions with a particular eigenvalue for S,. 


where 6;; is +1 for any ¢, which has @ on the 
orbit 7 and 8 on j; it is —1 when 8 is oni and a 
on j; or it is 0 when both 7 and 7 are associated 
with the same spins. Convenient methods of 
determining an independent set of bond eigen- 
functions have been developed.’ 

In order to work directly with these y’s which 
we call bond eigenfunctions it is necessary to 
have an operator for P;;. Diract has shown that 
(in units) 


3(14+45;-5;). 
(ss+s;)?- (s;?+5;?) 
Sits;=Siz, =F, 


we have 254°5;= 


(8) 
Since 


and 


Here s; and s; are the operators for spin vectors 


7 Rumer, Nachr. Ges. Wiss. Gottingen, M. P. KI., p. 
337 (1932). 
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of electrons 7 and j, and S;;? is the operator for 
total spin angular momentum associated with 
the two electrons. As a result of substitution in 
(8) we find: 


(1 (9) 

Methods have been developed (reference 5c) 
for operating with S? on any antisymmetric 
function gy. These are readily adapted to S;;? in 
the following way. As in the case of S?, 


+S,, gt (10) 
Here all operations are to be carried out on the 
electrons 7 and j only. The coefficient 4/27 is 
omitted on S,, i;, since wé are dealing in h/27 
units. The z’s multiplying S,, ;; are, of course, 
the square root of —1. 

The effect of the operation (1—.S;;7) on any ¢ 
is to interchange the spins of electrons i and j 
and multiply the result by —1, as is easily 
verified. Consequently, the effect of P;; in the 
integral (5) is to pick out from y, and y; the 
number of times that ¢’s of the two differ only 
by a single exchange of spins on the electrons 
1 and j. The negative sign corresponds to the fact 
that single-exchange integrals occur with signs 
opposite to that of the coulombic one. These 
results, of course, were to be expected, but the 
considerations just stated constitute a simple 
verification of the vector method of determining 
matrix coefficients. 

One could, of course, proceed to calculate 
coefficients of exchange integrals occurring in 


matrix components between bond eigenfunctions. 


by the laborious process expressed in (5) and (9). 
However, complete and simpler rules have been 
found already from entirely equivalent reasoning. 
The reader is referred to previous papers®® * ° 
where methods have been given for determining 
the coefficients from “‘sequences,’”’ each of which 
is related to the two bond eigenfunctions of 
a particular matrix component. The vector 
method, however, does lead to the interesting 
comparison with the sequences which follows. 


a ia Frost and Turkevich, J. Chem. Phys. 1, 777 
® Stearn, Lindsley and Eyring, J. Chem. Phys. 2, 410 
(1934). 
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The vector form of the operator P;; suggests 
rewriting (7) as 


W)ydr= (J- W-3 


i<j 
—2 D0 Kij(si-s;) re. 


i<j 


(11) 


We now write the sequence in an equivalent 
although formally somewhat different form than 
was used in previous work. Some simple con- 
ventions permit us to write the matrix component 
between two bond eigenfunctions of any multi- 
plicity. Unless the two eigenfunctions are of the 
same multiplicity, i.e., have the same number of 
bonds the matrix component is zero. Where 
they do have the same number of bonds we 
classify the orbitals in the following way. If two 
orbitals are bound together in one eigenfunction 
we put them into two different groups. Now if 
in the second eigenfunction either or both of these 
orbitals are bound to new orbitals the new 
orbitals are assigned to the group not containing 
the bond partner. This process is continued until 
the assignment of orbitals starts repeating itself 
or until a new orbital found is not bound to a 
partner. If the first case (repetition) arises the 
two groups are said to form a cycle and in the 
second case they are said to form a chain. After 
completing a cycle or chain the process is then 
repeated with some orbital not yet appearing 
until all orbitals in the two bonds are classified 
as cycles or chains. Orbitals not bound in either 
bond eigenfunction are counted as chains of one 
link. Any such matrix component can then be 
written as 


w- 1/2 K,; 


m cycles 
+1/2(Ki;"” —K,;'""))} 
+6u2 


u even 
chains 


chains 


(12) 


The symbol 6,; is, of course, 1 if w=7, and 0 if 
u+%i where u is the number of even chains. 
For the unprimed exchange integrals, K;;, the 
z orbital is to be taken from one group in a par- 
ticular cycle or chain and j is to be taken from 
the other group. In K;,;’ both orbital z and j are 
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to be from the same group. For K;;’’ the 7 and 7 
must be from different chains and if one is from 
the first group of a chain the other must be from 
the second group. For K;;’"’ the 7 and 7 must be 
from different chains and either both from first 
or both from second groups. It is immaterial 
which group in a cycle or even chain is counted 
as first or second. For an odd chain an orbit 
unbound in either eigenfunction belongs to the 
first group. This determines the second group. 
Having chosen the second groups we use them 
to determine p in the following way. The ex- 
ponent # is the total number of times that second 
group orbitals are written as the first letter of a 
bond in either eigenfunction. J is the sum of the 
coulombic integrals; W is the energy and LKii 
< 


is the sum of all possible single exchange in- 
tegrals. The exponent m is the number of cycles. 
The justification of (12) is readily obtained 
along the lines developed in reference (5c) for 
singlet states. Since for singlet states there are 
no chains (12) then reduces to (—1)?2"[M 
+ > 3/2(Ki;—Kzi;’) (13) where the expression 


m cycles 
has been simplified by making the substitution 
M=Q-W-1/2>K;x;. 
<j 
We now give two matrix components as 


examples of (12). First {(ab cd eg f hijkl mn) 
X W) (ac bd ef g hij kl m n)dr=(ad/bc, fg/e, 
h, m, m, (i/j)(k/l)) 
Secondly (ab cd eg fhijklmn)(H—W)(ac bd 
ef ghijklmn)dr=(ad/bc, fg/e, h, k,l, m, n) 
= 
—ad—bc —fg} +1/2{eh+ek+el+em+en—fk—fl 
—fm—fh —fn — gk — gl — gm — gn —kl —km— gh 
In the 
bond eigenfunctions two letters written with no 
space between are understood to be bound 
together. A letter separated by spaces from both 
neighbors is not bonded. In the first case a is 
bound to 6 in the first eigenfunction so they are 
assigned to different groups in the sequence. 
Now in the second eigenfunction } is bound to d 
so they are assigned to different groups and 
finally in the first eigenfunction d is bound to c. 
Continuing we simply repeat our letters so that 
we have the cycle ad/bc. Similarly we have the 
odd chains fg/e and h and m, etc., and finally 
the two even chains (i/j) and (k/l). The paren- 
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theses are used to distinguish even chains from 
cycles. The application of (12) to our expression 
(sequence) which classifies our orbitals into 
cycles and chains leads immediately to the 
expanded matrix component. In the expanded 
form two letters written together signify the 
corresponding exchange integral between the two 
orbitals. 

The simple ratios between the coefficients of 
W and the coefficients of the various primed and 
unprimed exchange integrals Ki; in (12) could 
in fact have been deduced from (11). Re- 
membering that S,, S, and S, have the forms 


spectively, when written as the matrices oper- 
ating on the” eigenfunctions a and 6 of S,. 
We obtain for the corresponding operators 
S2i°S2;, SyitSy; and their sum S;-S; 
operating on ai8;, a:8;, Bia;, 8:8; the matrices 


0 0 0-1 1000 
0010 0-1 2 0 

(: 1 0 (: 2-1 
-1 00 0 0001 


Our problem is the comparison of matrix com- 
ponents for twice the last operator with those for 
the unit matrix for bond eigenfunctions. The sub- 
scripts 7 and j on the spin eigenfunctions a and 6 
are to indicate that these are the spin eigen- 
functions associated with the ith and jth orbitals. 
As we have seen bond eigenfunctions are linear 
combinations of the antisymmetric configuration 
eigenfunctions g; of Eq. (1) with coefficients of 
+1 or 0. Because of the symmetry of the energy 
operator /7 with respect to interchange of elec- 
trons we have the well-known equality 


ocr 


ff 
f (N) (13) 


where Q,(aa(1)b8(2)-+-+) is simply the diagonal 
element of ¢,. To calculate the coefficient of the 


1© Pauli, Zeits. f. Physik 43, 601 (1927). 
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exchange integral ab in (13) we need simply 
replace H by —25S,-S2 with a similar procedure 
for any other coefficient. To calculate the matrix 
component of unity H is replaced by 1. Now the 
matrix components for bond eigenfunctions is 
clearly a sum of terms like (13). An examination 
yields the following results for the various cases. 
The coefficients of K;; for cycles arises from the 
sum of matrix components of the operator 
—2S;- 5S; for four types of integrals (13) occurring 
equally often. We give the spin eigenfunctions 
associated with the ith and jth orbital for both 
¢u and yg, for the four types (1) ai8;, ai8;; 
(2) —Bia;; (3) —Bia;, a8; and (4) —Bia;, 
—B;a;. The only contributions to J,, arises from 
the first and last of the four pairs so that the 
ratio —2(S;-S;)-1+I,+ is readily seen to be 3/2. 
For K;; in odd chains two types of integrals are 
different from zero. Here as in all other cases all 
integrals occurring at all occur equally often, 
i.e., 2” times. The spin eigenfunctions for ¢, and 
are a:8;, and so that 
again For K;;' in a cycle 
four types occur a;a;, aia;; BiB;, BiB;; Bia;; 
Bia;, a8; giving a ratio —3/2. For K;;’ in an 
odd chain there are the two types aia;, — aia; 
and —£,a;, a8; again giving the ratio —3/2. 
For K;;'’ between two odd chains we get Bia;, 
Bia;, giving the ratio 1/2; for K;;’’’ between two 
odd chains a,a;, asa; giving the ratio —1/2. 
Finally for K;;’ between two even chains we 
have giving zero for the matrix 
component of J,; and 1 for the matrix component 
of —2S;-S;; and for we have —Bja;, — 
giving the values 0 and —1 for the respective 
matrix components. The coefficients in all cases 
are of course multiplied by the common factor 
(—1)?2™. 

An additional feature of the present method of 
calculating matrix coefficients is the ease with 
which multiple exchange coefficients can be calcu- 
lated. Any multiple permutation of orbits can 
always be expressed as the product of single 
exchanges ; for example P;;,=P::P;;. The coeffi- 
cient of Ki;,— WaA;;, is most conveniently found 
by taking advantage of the observation, similar 
to that used by Heitler," that P;; exchanges the 
orbits 7 and 7 in the bond eigenfunction upon 


" Heitler, Zeits. f. Physik 79, 143 (1932). 
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which it operates, changing the sign of the 
function as well. One can then use P;;= —1/2 
X (1+4s,-5;) to determine the multiple exchange 
coefficient after the usual method of sequences. 
In the second matrix component given consider 
Pave=PeaPas. Then P.s(ac bd ef ghijklmn)= 
— (bc ad ef gh ij kl m n) determines the sequence 
with (ab cd eg fhijkimn) to be now —(ac/bd, 
1/j, fg/e, h, k, l, m, n), from which —1/2(1 
+4s.-5a)r is seen to be 2(—1)42?. Hence in this 
instance the coefficient multiplying (Kase — W Aasc) 
is 8. More complex changes in the sequence arise 
when the exchanges involve interchanges of 
orbits coming from different cycles or from a 
cycle and a chain, but the procedure is the same. 


SINGLET ENERGY CALCULATIONS FOR 
SYMMETRICAL SYSTEMS 


Most stable chemical systems are in singlet 
states, and energy levels of this multiplicity 
assume particular importance therefore. By com- 
bining bond eigenfunctions into new ones in a 
way indicated by group theory from the atomic 
configuration, it is possible to factor the singlet 
secular equation into subdeterminants of smaller 
degree (for examples see reference 8). Usually, 
however, one is most interested in the lowest 
energy level, and then much simpler methods, 
involving no very great use of group theory, are 
available. Usually also the lowest level results 
from solution of the determinant related to the 
combinations of bond eigenfunctions belonging to 
the identical representation. 

Seitz and Sherman” have indicated a fairly 
simple way to determine eigenfunctions of the 
identical representation. An equivalent and 
somewhat simpler method is described by Stearn, 
Lindsley and Eyring.’ The expression 

R 


where xz is the group theory character belonging 
to the symmetry operation R for a given one- 
dimensional representation and where y is any 
bond eigenfunction, can be used to determine the 
required number of #’s of this symmetry. For 
the particular identical representation all xr’s are 
unity, and the determination of #’s is extremely 
simple: Ry. The procedure is to determine 
R 


12 F. Seitz and Albert Sherman, J. Chem. Phys. 2, 11 
(1934). 
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VECTOR AND BOND EIGENFUNCTION METHODS 


all possible ®’s, one from each of the independent 
set of singlet bond eigenfunctions. Some of the 
resulting eigenfunctions (#’s) will be the same, 
except possibly for sign, while others though 
different may not all be independent. It is a 
simple matter, however, to pick out from these 
an independent set of #’s. 

This procedure can be further improved by 
drawing the independent bond eigenfunctions on 
a figure which has the symmetry of the problem 
under consideration. When this is done it be- 
comes apparent immediately that certain y’s, 
with the same pattern of bonds, differ from each 
other by permitted symmetry operations. These, 
then, when given the proper sign and summed 
will form part of a particular #, though not 
necessarily all of it. Certain other eigenfunctions 
with the same bond pattern will also enter the 
® function, and these additional ones will corre- 
spond to “crossed-bond” functions when drawn 
on the circle of orbits used to determine the in- 
dependent set. These can be resolved in terms of 
the independent y’s, and from all #’s so found, 
equal in number to the number of different 
patterns formed by the independént set when 
drawn on the figure with the symmetry of the 
problem, the largest number of independent 
functions (#’s) is selected for the construction of 
the secular equation. 


RULES FOR RESOLVING CROSSED BOND 
EIGENFUNCTIONS 


Simple empirical rules for the uncrossing of 
bond patterns containing as many as five bonds 
have been given." It is possible in the simple 
way which follows to extend these methods to 
twelve-electron eigenfunctions. For these six- 
bond cases one traces along bonds of the crossed 
eigenfunction A (never traversing the same bond 
twice) the shortest connections between the 
extreme points of the bonds in the uncrossed 
eigenfunction B. If this is possible the expansion 
of A contains B at least once. If three cr more 
paths pass through a common point the coeffi- 
cient of B is p—2+g. What is meant by the 
number of paths p passing through a common 


8 Eyring and Sun, J. Chem. Phys. 2, 299 (1934). 


. 1. Examples of the methods of tracing bond paths 
on crossed-bond eigenfunctions (A) in order to determine 
the number of times (n=p—2+gq) a given uncrossed 
eigenfunction (B) is contained therein. The value of p—2 
is to be taken as unity for values of p< 3. Filled in areas are 
those bounded by unused sections of bonds and are to be 
regarded as essentially points in determining p and g. Pairs 
of A eigenfunctions connected by plus signs are those ob- 
tained by uncrossing one of the crosses of the function 
immediately preceding the pair. 


point needs to be defined. Frequently parts of 
the bond lines in A are not traversed in drawing 
the required shortest paths. All points connected 
by untraversed lines are for this purpose to be 
counted as one point. Also any path which 
passes through a point and which with the 
perimeter circle separates other paths common 
to the point into two groups is omitted in de- 
termining » for the point. This discarding is 
continued until there are no more such paths. 
Since q is zero in cases of five or fewer bonds, 
it will also be zero when the six bonds of A can 
be divided into two or more non-intersecting 
groups. Otherwise g is determined by the number 
of pairs of diametrically opposite bonds on the 
perimeter of B which when drawn on A have 


- the following properties. For a given pair to be 


counted in the determination of g each member 
of the pair must have its path as drawn on A 
terminate along two bond lines which cross; 
the paths for members of such a pair must 
touch at some point, either directly or by the 
“shrinking” of untraversed lines as discussed in 
the previous rules; and the pair members must 
not be separated by the traced path of another 
bond, for one or two such pairs g is 1 or 2, 
respectively. For three such pairs q is 4. Self- 
explanatory examples are given in Fig. 1. 
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ENERGIES FOR HIGHER MULTIPLICITIES 


The method of using bond eigenfunctions 
drawn on a figure with the symmetry of the 
actual problem at hand combined with the above 
rapid method of resolving crossed eigenfunctions 
thus becomes a simple one for the determination 
of singlet levels of the identical representation. 
A similar procedure is valid for levels of higher 
multiplicity and the identical representation. 
The method for selecting independent sets of 
bond eigenfunctions for higher multiplicities 
follows.'* Between two of the orbits in the usual 
circular arrangement one inserts an additional 
orbit. One then draws non-crossing bonds within 
the circle, leaving the extra orbit and v—1 
others unbonded, where v is the multiplicity. 
If it is possible in a particular eigenfunction to 
draw lines between the extra orbit and all other 
unbonded ones without crossing any of the 
bonds, the eigenfunction is a member of the com- 
plete set. After selecting in this manner a 
complete set, the extra orbit is, of course, dis- 
carded. The ‘“‘pattern” method of determining 
’s as combinations of members of this set can 
then be applied, but in general the levels thus 
determined may not contain the lowest level of 
the system. 


USEFULNESS OF SPACIAL SUPERPOSITION 
PATTERNS 


It is possible to use the patterns to determine 
&’s belonging to other representations as well if 
the group theory characters are known. The 
method is then still superior to an experimental 
choice of y’s from which to build the ®’s by use 
of the symmetry operations and characters. In 
fact one immediately sees that the ©’s will always 
be composed of bond eigenfunctions of the same 
pattern with coefficients depending upon the 
representation. 

Having selected an independent set of ’s for 
the problem in hand, one now finds that the use 
of the patterns is of further value in the calcula- 
tion of matrix components. For this purpose it is 
more convenient to use the %’s expressed in 
terms of bond eigenfunctions which are of the 


4 Rumer, Teller and Weyl, Nach. Ges. Wiss. Gottingen, 
M. P. KI. p. 519 (1932). 
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same pattern, irrespective of whether they are 
members of the independent set or not. One then, 
in a manner introduced by Pauling®4 for orbits 
arranged in a circle, superimposes bond eigen- 
functions upon each other in order to determine 
matrix components as follows. Suppose ®, and , 
are the two of the ®’s for which a component is 
to be determined and y,, and ws are typical 
bond eigenfunctions composing the two, re- 
spectively. Then the matrix component /®, 
xX W)® = S Dulrubru( — W)> 
is the sum of ones for which coefficients can be 
found by the sequence method. The super- 
position patterns of y,, and Wy are simply 
diagrammatic ways of writing the sequences, the 
parts of which will be recognized as follows. 
Certain “‘islands’’ composed of a broken closed 
line with an even number of bonds are cycles. 
Unclosed lines of no bonds or of an even or odd 
number of bonds are chains. First group orbitals 
are alternate ones starting with the first orbital 
in an odd chain and any arbitrary orbital in an 
even chain or cycle (island). The other orbitals 
belong to second groups. As long as one is 
dealing with “canonical’’ eigenfunctions, i.e., 
those formed according to the previous definition 
with orbits of each bond written in alphabetical 
order, it is unnecessary to consider the quantity 
(—1)? in (12) as other than positive. The 
coefficients of J, W, and the exchange integrals 
are then directly determinable for each super- 
position pattern from (12). The advantages of 
the spatial superposition patterns thus employed 
are that (1) any one of the y,.’s, apart from the 
coefficients c,, and Ci, gives the same types of 
superposition patterns with all of the y,,’s as 
does any other y,,, hence the same coefficients of 
J, W, and the exchange integrals; and (2) it is 
possible to allow directly for several exchange 
integrals being equal as may be indicated by the 
spacial arrangement of orbits. These factors 
result in reduction of the labor of calculating 
matrix components. 

The application of the above methods of 
obtaining eigenfunctions and matrix components 
for the identical representation from figures 
drawn according to the symmetry of the problem 
is particularly simple when the atoms are 
arranged in one or two-dimensional configura- 


ti 
: 
b 
ti 
tk 
th 
al 
di 
de 
di 
né 
th 
to 
in 
sin 
va 
atc 
Th 
ene 
atc 
the 
sio 
cir 
ato 
gre 
sim 
inti 
lim 
| ene 
hov 
sho 
in t 
I: 
are 
moc 
non 
ator 
obv 
will 
four 
The 


VECTOR AND BOND EIGENFUNCTION METHODS 


tions, such as a straight line, circle, square, etc. 
Then it is possible to draw an independent set of 
bond eigenfunctions directly on this configura- 
tion, and no symmetry operation on any one 
of these will give rise to eigenfunctions outside of 
the complete set. The number of different ®’s is 
then equal to the number of patterns that arise, 
and all of them are independent. When three- 
dimensional figures are necessary no way of 
determining the independent bond functions 
directly is possible because of the two-value 
nature of the spin eigenfunctions. In such cases 
&’s will contain bond eigenfunctions other than 
those of the independent set, and these will need 
to be resolved in order to pick out independent 


ENERGIES OF FILAMENTS OF SopIUM ATOMS 


In order to apply the foregoing considerations 
in a simple way we have calculated the lowest 
singlet level for the identical representation for 
various even-numbered configurations of sodium 
atoms arranged on straight lines or in circles. 
The object was to find a limiting value for the 
energy per atom approached as the number of 
atoms increases. This energy would then be 
that of a single atom in an infinite linear exten- 
sion of atoms, a one-dimensional lattice. The 
circle as well should show the same energy per 
atom in the limit, and since the circle has the 
greater symmetry it was hoped that the greater 
simplification of the secular equations thereby 
introduced would allow closer approach to the 
limiting value by permitting determination of 
energies for larger numbers of atoms. Actually, 
however, the straight line of atoms more quickly 
shows a trend which indicates the desired energy 
in the limit. 

In Fig. 2 the independent sets for six atoms 
are shown as drawn on linear and circular 
models. The method in the linear case is to draw 
non-crossing bonds on the line connecting the 
atoms or to one side of that line. It becomes 
obvious immediately that a quadratic equation 
will give the energies for the circle, while one of 
fourth degree arises for the linear configuration. 
The #’s in the circular case are ®,=I+II, 
$,=I]I+IV+V; and for the linear model 6, =I, 
$, =I], 6; =III+IV, and =V. 


/N 


Fic. 3. 


Pauling®* has given superposition patterns for 
the circular configuration, and in Fig. 3 are 
given several for the linear set. The examples 
have been chosen to illustrate the existence 
of apparently different superposition patterns 
which, however, are easily recognizable as giving 
the same expansion in terms of the coulombic 
and exchange integrals. All of the examples 
given yield F=2'[M+3/2(5a+3y+6) —3/2(48 
+25) ], where a, 8, y, 6 and are the exchange 
integrals between orbits or atoms separated, 
respectively, by 1, 2, 3,4 and 5 times the distance 
between adjacent atoms. The matrix component 
JS%,(H— W) = F and JS W)®ydr 
=2F. The other matrix components are de- 
termined similarly. 

In this manner the secular equations for 
circular and linear configurations up to 12 and 8 
orbits, respectively, have been set up. These are 
of twelfth degree for the largest circle considered 
and of tenth degree for 8 atoms on a line. 
Beyond these numbers of atoms the determinants 
are entirely too large to be handled with reason- 
able labor. We shall not repeat the matrix com- 
ponents here, since the determinants are of 
limited application. 

For estimation of coulombic and exchange 
integrals we have employed the method pre- 
viously used in similar investigations.“ * The 
dissociation energy of the Naz molecule, D, is 


% Taylor, Eyring and Sherman, J. Chem. Phys. 1, 68 
(1933). 
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Fic. 4. Energies of various sodium atom configurations. 
Coulombic contributions. — — Exchange contri- 


butions. © atoms arranged in circle. A atoms arranged on 
line. Curves without circles or triangles after Slater. N is 
the number of atoms. 


17.5 kg cal., D’=17.7 kg cal., wo=159 cm“, 
and =3.15A. From these constants a Morse 
curve was constructed. For the coulombic bind- 
ing between two atoms at a given distance from 
each other we have employed Rosen’s value’® of 
28.3 percent of the total binding energy as de- 
termined from the Morse curve. The corre- 


sponding exchange integral then is 71.7 percent | 


of the binding energy. Energies have been calcu- 
lated only for configurations in which the dis- 
tance between adjacent atoms is 7p for the Nae 
molecule. This is sufficiently close to the equi- 
librium configuration in each case. 


TABLE I. 
Circle Line 
J(circle) W-J Ww W/N J W-J Ww WIN 
2 17.7 8.85 17.7 8.85 
4 25.6 11.2 36.8 9.20 16.5 19.1 35.6 8.90 
6 40.1 14.5 54.6 9.10 28.0 25.4 53.4 8.90 
8 50.4 21.3 71.7 8.96 39.5 31.8 71.3 8.91 
10 62.6 29.5 92.1 9.21 
12 73.0 38.0 111.0 9.25 


The results are given in Table I and plotted 
in Fig. 4. It will be noticed that the coulombic 
and exchange energies for the linear configura- 
tions increase by practically constant amounts 
for each pair of atoms added. Hence the total 
energy per atom for these cases rapidly ap- 
proaches the constant value of about 8.91 kg cal., 
which is the desired limiting energy value. The 


1 Rosen, Phys. Rev. 38, 255 (1931). 
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circular configurations possess greater coulombic 
energy than do the linear ones. On account of 
the unusual importance of coulombic energy, the 
circular arrangements of the atoms appear to be 
the stabler ones. 

The straight lines for coulombic and exchange 
energy plotted against N for the linear con- 
figurations enable one to express the total energy 
by the following equation : 


W=5.0+5.7(N—2)+3.2(N+2). (14) 


The first two terms arise from coulombic origins, 
those of the third from exchange contributions. 
It becomes apparent from (14) that the energy 
per atom for infinite N is composed of 5.7 kg cal. 
from coulombic sources and 3.2 kg cal. from 
electron exchanges. 

These results can be compared with those of 
Slater,” who found for an infinite linear lattice 
of identical atoms to the approximation of in- 
cluding only interactions between adjacent 
atoms: 


W=(N—1)J'+0.290(.N —3)K, (15) 


where J’ is the coulombic attraction between 
adjacent atoms and K the exchange integral 
between adjacent orbits. By using the values 5.0 
and 12.7 for J’ and K, respectively, the results of 
(15) have been calculated for small values of NV 
and plotted in Fig. 3. Though it appears that 
Slater’s equation is not very good for small 
values of N, it may be written 


W=5.0(N—1)+3.7(N—3), 


and it then becomes evident that radical differ- 
ences between (14) and (15b) disappear for 
large N. At such large values of N the Slater 
equation indicates somewhat less coulombic 
energy and a little more exchange energy than 
the method used here gives, but these differences 
nearly cancel to give the limiting value for the 
total energy per atom of 8.7 kg cal. as compared 
to the 8.9 determined by including interactions 
between atoms further apart than neighbors. 
The higher value is, however, greater than the 
energy per atom in the diatomic molecule, while 
the lower is not. . 


17 Slater, Phys. Rev. 35, 509 (1930). 
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The Activated Complex in Chemical Reactions 


HENRY EyrinaG, Frick Chemical Laboratory, Princeton University 
(Received November 8, 1934) 


The calculation of absolute reaction rates is formulated 
in terms of quantities which are available from the potential 
surfaces which can be constructed at the present time. The 
probability of the activated state is calculated using ordi- 
nary statistical mechanics. This probability multiplied by 
the rate of decomposition gives the specific rate of reaction. 
The occurrence of quantized vibrations in the activated 


complex, in degrees of freedom which are unquantized in 
the original molecules, leads to relative reaction rates for 
isotopes quite different from the rates predicted using 
simple kinetic theory. The necessary conditions for the 
general statistical treatment to reduce to the usual kinetic 
treatment are given. 


INTRODUCTION 


HE customary procedure for calculating bi- 
molecular reaction rates has been to esti- 
mate the number of collisions between reacting 
molecules by using a cross-sectional area taken 
from measurements on momentum transfer. Such 
a cross section bears no very clear relationship to 
the area within which two molecules must fall 
in order to permit exchange of partners, i.e., to 
transfer mass. The violence of the collisions is of 
different orders of magnitude for one thing, and 
it is quite clear that many collisions which might 
result in momentum transfer are not oriented 
properly to permit exchange of atoms. This last 
difficulty is ordinarily met by introducing an 
empirical steric or orientation factor to take 
care of whatever discrepancy may arise between 
the observed and assumed collision area. This 
factor is often between 1 and 10 but may be 
as small as 10-*.! We propose here to obtain 
explicit expressions for the reaction rates. 

The ideas underlying the present calculations 
are the following ones. The forces between atoms 
are due to the motion and distribution of 
electrons and must be calculated, therefore, 
using quantum mechanics. However, after this is 
done the nuclei themselves can be assumed to 
move under the influence of these forces accord- 
ing to classical mechanics. It must be possible, 
therefore, to calculate the reaction rates by the 
methods of statistical mechanics (or kinetic 
theory), if one assumes the aforementioned forces 
to be known. This is what is done in the present 
paper using a modification of the schemes de- 

' Moelwyn-Hughes, Kinetics of Reactions in Solution, 


Oxford Press (1933). W. E. Vaughan, J. Am. Chem. Soc. 
55, 4115 (1933). 


veloped by Herzfeld,?, Tolman* and Fowler‘ 
among others and more recently applied in a 
very interesting way to the ortho-para hydrogen 
conversion by Pelzer and Wigner.® 

Cases occur when classical mechanics does not 
apply to the motion of the nuclei. Zero point 
energy may be present for some vibrations, and 
it will be necessary to deal with quantized 
vibrations in a semi-classical way. Tunneling may 
occasionally play some role in the motion. In 
other cases, probably also of very rare occurrence, 
there may be jumps from one energy level to 
another. The latter factors may also change the 
results calculated by neglecting them by orders 
of magnitude as, e.g., in the case of N.O.* We 
are not concerned here with reactions in which 
the last two effects are important. 

We now consider in more detail the nature of 
the surfaces and the motion which corresponds 
to a reaction. A group of atoms may of course 
arrange themselves in an infinitely large number 
of ways. If the energy of such a system of atoms 
for the lowest quantum state of the electrons is 
plotted against the various distances between the 
nuclei, we obtain the potential surface which 
governs (except in the aforementioned cases) the 
motion of the nuclei. 

Now a system moving on this surface will 
have kinetic energy which may be quantized 


for the different degrees of freedom in a variety 


K. F. Herzfeld, Kinetische Theorie der Warme ;(Miiller- 
Pouillets Lehrbuch der Physik) 1925. 
Tolman, Statistical Mechanics, Chemical Catalog 

0., A 

*R. H. Fowler, Statistical Mechanics, Cambridge Univ. 
Press, 1929. 

5H. Pelzer and E. Wigner, Zeits. f. physik. Chemie B15, 
445 (1932). 

® Volmer and Kummerow, Zeits. f. physik. Chemie B9, 
141 (1930). 
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of ways, consistent with the particular energy 
and the particular position on the surface. Low 
places in the potential surfaces correspond to 
compounds. If a particular low lying region is 
separated from all other low places by regions 
higher than about 23 kilocalories, the compound 
will be stable at and below room temperature. 
The higher the lowest pass the higher is the 
temperature at which the compound is still 
stable. A reaction corresponds to a system 
passing from one low region to another. In 
thermal reactions the Boltzmann factor makes it 
certain that the reaction will proceed by way of 
the lowest pass. The activated state is the highest 
point along this lowest pass. Before considering 
the activated state further we discuss the general 
problem of constructing our partition functions 
for a given surface. 


FORMULATION OF PARTITION FUNCTIONS 


We here simply sketch the procedure for the 
gas phase and indicate where the modifications 
for reactions in solution will come. If a system of 
atoms is represented by a point on a potential 
surface such that for the motion in every direc- 
tion the potential increases or remains constant, 
we may apply the well-known method of small 
vibrations to obtain the normal -coordinates.’ 
Now our system will show translational, vibra- 
tional and rotational degrees of freedom. The 
weighted number of states associated with a 
_ degree of translation are 


—p? _ 
aa: [exp 


2m ar)” h 
where q;, p; and m; are the generalized length, 
momentum and mass for this degree of freedom, 
respectively. For a solution or highly compressed 
gas the integration over the length /; would in- 
volve a potential function in an important way. 
Also the reacting system would now include in 
an essential way the degrees of freedom of a 
certain number of the solvent molecules. If we 
set the length /;=1 we have the number of unit 
cells per cm of length, a quantity frequently 
used in what follows. The other universal con- 


7See, for example, Whittaker, Analytic Dynamics, 
Cambridge Univ. Press, 1927. J. H. Van Vleck and P. C. 
Cross, J. Chem. Phys. 1, 357 (1933). 


108 HENRY EYRING 


stants require no definition. Similarly for a 
vibrational degree of freedom we obtain for the 
weighted number of states using classical theory 


d ( —faqi? kT 
ex i ex =— 


where (1/27) (fi/m;)*. 


Here f; and v; are the force constant and fre- 
quency, respectively. However, summing over 
the quantized states, we obtain the familiar 
quantum theoretical expression (1—exp (—hv;/ 
kT))— for the weighted number of states for the 
harmonic oscillator. For hvi«<kT of course (1 
—exp (—hv;/kT))“*=kT/hv;. The weighted sum 
of states for the two degrees of rotation of a 
linear molecule is 


oo —j(j+i)h? 
27+1) ex ( ): 
2. 


Here J is the moment of inertia. Now there are 
cases for which the position of the activated 
state on the potential surface depends on the 
quantum number j. Certain cases of this type 
will be treated in another place. If the critical 
configuration does not depend in an important 
way on j and if j(j+1)h?/8°I<kT, we can in 
the well-known way substitute for the summa- 
tion the integral 


1) (a IkT 


For a non-linear molecule we can write for 
the weighted sum of states the expression 87° 
(If the moments of in- 
ertia are small we should sum over the energy 
levels for a symmetric or asymmetric top.) 
The expression exp (—j(j+1)h?/(8r° 
X(ABC)*kT)) is satisfactory if the moments of 
inertia differ only a little and the corresponding 
integral 


yields the classical weighted sum given above. 
We now suppose that we have a great number 

of systems in thermal equilibrium and want to 

know the relative probability of any system 
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being at points 1 and 2 on a potential energy 
surface. We assume points 1 and 2 are of such a 
kind that we know the corresponding normal 
coordinates and that the energy of point 2 is 
greater by Ep» than that of point 1. The ratio 
of the corresponding weighted number of states 
of points 2 to 1 (which is the relative proba- 
bility) is obtained by multiplying exp (— Eo/kT) 
into the ratio of the product of the weighted sum 
for each normal mode for point 2 to the corre- 
sponding product for point 1. This very brief 
discussion of the equilibrium constant is sufficient 
for our purposes. For a rigorous deduction one 
may consult the work of Fowler previously 
referred to. The next conception we will require 
is the rate at which the systems at point 2 are 
moving along a particular normal coordinate. 
In our considerations point 2 will in general 
correspond to the activated state which we now 
consider in more detail. We call a system at the 
activated point an activated complex. 


THE ACTIVATED COMPLEX 


The activated state is because of its definition 
always a saddle point with positive curvature in 
all degrees of freedom except the one which 
_ corresponds to crossing the barrier for which it 
is of course negative. Further, the barriers are 
so flat near the top that tunneling may be 
neglected without appreciable error. A con- 
figuration of atoms corresponding to the acti- 
vated state thus has all the properties of a 
stable compound except in the normal mode 
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corresponding to decomposition and this mode 
because of the small curvature can be treated 
statistically as a translational degree of freedom. 
Thus a non-linear activated complex with n 
atoms (m=3) has three regular translational 
degrees of freedom corresponding to motion of 
the center of mass in addition to the one corre- 
sponding to passage over the top of the barrier. 
It also has three rotational degrees of freedom 
for the molecule as a whole, and the remaining 
(3n—7) degrees of freedom correspond to in- 
ternal rotations or vibrations. A linear molecule 
differs from this in that one of the degrees of 
freedom which was a rotation is instead a 
bending vibration. Now the calculation of the 
concentration of activated complexes is a straight 
forward statistical problem, given the moments 
of inertia of the complex and the vibration 
frequencies. This information is given with 
sufficient accuracy, even by our very approxi- 
mate potential surfaces, to give good values for 
the partition functions. 

The procedure for calculating the specific rate 
is the following: One first calculates the con- 
centration of activated complexes per unit length 
and with momentum / lying between and 
p+dp, both these quantities taken for the degree 
of freedom corresponding to decomposition. This 
is then multiplied by the associated velocity 
p/m* and summed for all values of momenta 
which correspond to passing over the barrier in 
the forward direction, i.e., for p=0 to ~. We 
now formulate the particular expressions for the 
various cases. 


THREE ATOM REACTIONS 


Consider the reaction 


A+BC~A—B—CAB+C, 


(1) 


where the activated complex A —B—C is linear. We write for the specific reaction rate constant 


C£ ala 


o3h? 


(2am3kT)} 3 —hv,* 


)) 


Lala m3 
( 


(2rme2kT)! —hv2\ \— 827 I2kT 
(FF) 


kT 


h? (1 —exp 


3 | 
(1—exp (—hv,*/kT)) 
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The subscript 1 refers to atom A, 2 to BC and 3 to the activated complex. Ep is the difference in 
energy between the initial substances and the activated state at the absolute zero. The quantities 
m;, I;o; are the mass, moment of inertia and symmetry number respectively of the particle 7. The 
three frequencies v;* (for i= 1, 2 and 3) are for the activated complex and 2 is the vibration frequency 
of AB. The letter g.i, and g,7, are the products of weights arising from electronic states and nuclear 
spin for the activated and ground state respectively. Frequently g.=g, and 7, and 7, need only be 
considered for reactions in which there is a change in ratio of ortho and para forms. m* and are 
the reduced mass and average momentum along the normal coordinate corresponding to decom- 
position. The quantity in Eq. (2) which multiplies p/m* (the average velocity of activated complex 
along the normal coordinate corresponding to passing over the barrier) is of course the concentration 
of activated complex per cm of length normal to the barrier when there is unit concentration of 
reactants per cc. Now 


p/m*= ( f / f exp 


Thus we see that the terms for the activated complex associated with the normal coordinate along 
which decomposition occurs give simply ((27m*kT)}/h)p/m* =kT/h. This factor will of course come 
in in the same way for every type of reaction. 

For some reactions it will happen that the same activated complex may cross the barrier and 
return without decomposing. This fact reduces the actual reaction rate. It is taken care of by the 
factor c which is the reciprocal of the average numbers of crossings required for each complex which 
reacts. It will generally be about unity. There are methods available for estimating the factor c. 
The other symbols appearing in Eq. (2) are well known and if we introduce our constants in grams 
cm sec. units, &; is given in cc molecules“ sec.—! If we want &; in cc moles~! sec.-! we must multiply 
the k of Eq. (2) by Avogadro’s number JN. 

A form frequently used for writing experimental specific reaction rates is 


ko= BT} exp (—E/kT)=BT' exp (—(Td log k2/dT—})). (3) 


The activation energy E in (3) is written as RT*d log ke/dT —4kT as this is the way it is calculated 
from the temperature increment. 
Putting (2) in the same form as (3) we find for comparison 


E=Ey+¥ (hos*/exp —1)—(hvs/exp (hos/kT) —1)—BT, 


i=1 


and 


BT'=c— 
LZntn 


ms )= h? (1—exp (—hv2/kT))N 


(24) (RT)! 3 


mym 
II (1—exp 


Xexp (exp hv;*/RT) —1)—'—hve/kT X (exp (hve/kT) —1)-!—-1 . (5) 


In (5) and in the preceding expressions the masses m; refer to the masses of a single atom or 
molecule. If the masses are taken in atomic weight units and written with primes we get: 
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BT = 1.41 (—hv2/kT)) 


x (1—exp (—hv,*/kT))— exp (exp (ho,;*/kT) 
i=1 i=1 
— (hve/kT) (exp (hve/kT) — cc mole~! sec.~! units. (5’) 


For many simple reactions all the terms in BT! except 1.41 X10” are of the order of magnitude of 1 
so that this factor multiplied by the activation energy term gives approximately the rate of the 
reactions and actually agrees to this approximation with the known experimental values. We here 
have an exact theoretical collision diameter for reaction rates which replaces the rough kinetic theory 
value and provides in addition a theory for predicting and explaining divergences. 

The inverse dependence on temperature of B is to be noted. This will be the case for molecules for 
which there are two more fairly stiff quantized vibrations in the activated complex than in the un- 
activated particles. Partition functions for fairly stiff vibrations depend only slightly on temperature. 
It is difficult to test this dependence of B upon T experimentally but for this same case there is a 
predicted dependence of isotopic reactions on mass which may be more readily detected. For the 
moment of inertia of the activated complex we have J; =m a?+-m,c? — (mya — mac)?/ms where m, and 
m, are the masses of atoms A and C, respectively, and m; is the sum of the masses of the three atoms; 
a is the distance between A and B, and c the distance between B and C. Now I2.=m4m;d?/(m4+-ms) 
where d is the distance between B and C whose masses are, respectively, m; and my. If all three atoms 
are alike [3/I2=(a+c)*/d’, i.e., the ratio is independent of the mass of the atoms so that for protium 
or deuterium mass enters explicitly into B only in the term (m;'/m,'m.’)!. Thus from this cause alone 
the protium reaction should go faster by a factor of 2! instead of 2! as is found if one assumes the only 
difference lies in the relative velocities of colliding particles. If the two extra bending frequencies in 
the activated complex are small, the term (1—exp (hv;/kT))—' approaches the value kT /hv;, and 
remembering that v;=(1/27)(f;/m;)' where f; is the corresponding force constant and m; the re- 
duced mass, we see that B will be dependent upon temperature and depend on mass in the way pre- 
viously supposed, i.e., B~ 7 /m!'. The actual dependence will lie between this extreme and B~1/Tm!. 
The results of Topley and the present author, soon to be published, indicate that actually the system 
should more closely approach the latter dependence. These authors have calculated values of B for 
both the ortho-para-hydrogen conversion and the reaction Br+He: which agree with experiment. 

It should be emphasized that in the particular formulation (2) of our specific reaction rate constant 
all the quantities may be calculated from the appropriate potential surface which can always be 
constructed at least approximately. However, in cases where any part of the partition functions is 
more accurately known from some other source such information can of course be incorporated. 

If one realizes that the term (27mkT)'/h for a translational degree of freedom has a value of the 
order of 10° for a light atom at ordinary temperature, one sees that the replacement of two terms 
like this in the initial products by two bending vibrations in the activated complex introduces a 
factor in k of the order 10~-"*. This factor is of course the analog of the collision area in the rough 
kinetic picture and explains why the kinetic picture works approximately, since the other terms in 
(2), except the average velocity, j/m*, are not very different from unity for many reactions. A dis- 
cussion showing when the general statistical treatment reduces to the usual kinetic theory treatment 
is given farther on. 

In the case where the least activation energy corresponds to a non-linear activated complex 
the terms 


/(osh*) (1 exp 


| 
| 
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in (2) are replaced by the quantity 


In this case we see a vibrational degree of freedom has become a rotational one. The letters A *B*C* 
now correspond to the principal moments of inertia of the activated complex. The calculation of 
the reaction rate and of B and E then proceeds exactly as before. 


Four ATOMS 
Consider the reaction: 


C—D 


(0) 


2°A 3B3C3) (RT) 3037 (1 —exp (—/v,/kT)) ORT /h Xexp (— Eo/kT) 
X [gntn(2 kT) —exp —exp (—hve/kT))~! 


K 10!2(300/T7) (A 3’B3’ C3’) 


—exp (—/hv, (—hv2/kT)) exp (— Eo/kT) 
in cc mole"! sec.~! units. (7) 


The subscripts 1, 2 and 3 refer to the species A—B, C—D and the activated complex respectively. 
Primed symbols as before mean that atomic weights are to be used as masses, and atomic weights 
x (Angstroms)? are the units in which moments of inertia are to be given. The significance of the 
other symbols will be clear from the definitions in connection with Eq. (2). Since all the quantities in 
(7) are of the order of unity for most reactions except the numerical factor and e~*»/*” we see again 
why the kinetic theory picture agrees approximately. 

Clearly, for any bimolecular reaction, we can immediately write down the expression corresponding 
to ks. If A, B, C and D instead of being atoms are radicals there will simply be additional vibration 
and internal rotation terms with the expression for moments of inertia of the initial substances in 
the appropriate cases, taking the form for non-linear molecules. The same reasons for approximate 
agreement with kinetic theory will remain. Now the activated complex for the type of reactions we 
are considering is the same for the forward or the reverse reaction. So that in calculating the specific 
reaction rate constant k for a unimolecular reaction which is bimolecular in the reverse direction we 
simply modify the denominator of the reverse reaction constant to correspond to the new initial 
reactants. 


UNIMOLECULAR REACTIONS 


Suppose we have a non-linear molecule of n atoms decomposing unimolecularly. We then write, 
cancelling out factors common to the initial and activated states: 


an—7 3n—6 


kg=c*a/o*(A*B*C*/ABC)! II (1 —exp (—hv*/kT))™ (1—exp (—hv;/kT)) 
(kT /h) exp (—Eo/kT). (8) 
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Quantities referring to the activated state in (8) are starred. Now in the particular case where 
hv;<«kT, i.e., all vibrational degrees of freedom approach a classical behavior we have (1—exp 
x (—hv;/kT))'=kT/hv; and (7) takes the form: 


3n—6 3n—7 


ky=c¥(o/o*) I (A*B*C*)(A BC) exp (—Eo/kT). (9) 


c* has the same meaning as the c defined in connection with (2). We of course come to this same 
result (9) directly if we integrate the appropriate classical expressions for vibration over phase space. 
Thus for each vibrational degree of freedom : 


(1/h) exp [exp RT 


if we use the relationship v;= }2(f;/m,)'. In using (8) it must be remembered that for certain mole- 
cules some of the degrees of freedom treated as vibrations can better be treated as internal rotations. 
In any particular case there is no particular difficulty in doing this. Eq. (9) is sufficiently near to 
that found for unimolecular reactions at high pressures that there seems no doubt of the wide applica- 
bility of both (8) and (9). A formula very similar to (9) was obtained by an approximate method in a 
paper by Polanyi and Wigner.*® 


General case 


Cases could of course be multiplied almost indefinitely but enough examples have been given to 
leave no doubt of the proposed method of procedure in a particular case. We may write for the 
specific reaction rate constant for a reaction of any order 


(10) 


where F, is simply the partition function (or Zustandsumme) for the activated state and F, is the 
same quantity for the normal state. F,’ is the partition function for the activated complex for all the 
normal coordinates except the one in which decomposition is occurring. The partition function for 
this normal coordinate is included in the factor (kT/h)e~*:'*”. The other quantities have been 
defined. 

The frequently observed negative temperature coefficient of trimolecular reactions has a ready 
explanation from the point of view presented here. Since the formation of an activated complex from 
three molecules involves a great loss in entropy, a reaction which goes with a reasonable rate at 
ordinary temperatures will necessarily have a low activation energy. Further in forming the complex 
a number of translational and rotational terms with direct dependence on temperature are converted 
into vibrational terms with very low temperature dependence. Thus the rate will vary inversely with 
the temperature to a comparatively high power so that if the activation energy is low enough, the k 
will have a negative temperature coefficient. Kassel® has discussed such reactions also. 

The extremely low rates, as compared with expectations from kinetic theory, observed in solu- 
tions' are to be thought of as associated with a change of translational or rotational degrees of 
freedom of the original molecules into vibrational or oscillational states of the activated complex. 
These changes may of course be in the enveloping solvent molecule. Any advantage of this formula- 
tion of the problem for solutions over any other consideration of entropy and heat content must 
come from a happy choice of the mechanism of reaction. 

When the rate determining step shifts to the collision process as it does for ‘‘unimolecular re- 
actions” at low enough pressure we again use well-known statistical methods, but our slow process is 
now connected with energy transfer in collision. 


®M. Polanyi and E. Wigner, Zeits. f. physik. Chemie A (Haber Band), 439 (1928). 
*L. S. Kassel, J. Phys. Chem. 34, 1777 (1930). 
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KINETIC THEORY DIAMETERS 


It becomes a matter of considerable interest to show under what circumstances the preceding 
general statistical method reduces to the simple kinetic theory scheme as ordinarily applied. We first 
calculate the number of collisions between two kinds of hard spheres A and B with radii 7; and r2 and 
masses m,; and mz, respectively. We use our general method. The respective concentrations per cc 
of A and B are N; and Nz. Our procedure is to calculate the number of pairs of molecules per cc per 
second which come closer to each other than (7:+72+ 6). We then let ¢ approach zero. Our collision 
complex then has three degrees of freedom associated with translational motion of the center of 
gravity ; one degree corresponding to relative translation along the line of centers; and two degrees 
corresponding to motion perpendicular to the line of centers, i.e., two degrees of rotation. Before 
collision there are six translational degrees of freedom, i.e., three for each sphere. The expression for 
the number of collisions when there is one molecule of each kind per cc may then be written at once: 


ku= (11) 


The significance of each term will be clear from our previous discussion. 
Now if the temperature is not too low we have j(j+1)h?/8x°I<kT ; so that we can make the usual 
approximation for the two rotational degrees of freedom, i.e., 


(2T-+1) exp = 80° LRT 


Also kT /h is just the term [(217m*kT)!/h ]p/m* of course. The moment of inertia J = (mym2/m,+mz2) 
X(n+1r2)?, so that we have after simplification 


(12) 
The number of collisions per cc per second is then 
Z= (mi (13) 


which is the usual expression for the number of collisions. Our method of treatment of collisions 
neglects certain of the refined features arising from the wave nature of the atoms. These are not of 
interest to us in our present treatment of reaction rates since here we make no explicit use of kinetic 
theory diameters. For an exposition of these features see a series of papers by Massey and Mohr.” 
For identical colliding systems a symmetry number should be included in (11) to (13). 

It is now easy to see when we are justified in using the simple kinetic picture. If the two colliding 
molecules have (a) none of their internal frequencies appreciably modified in the activated state and 
(b) if the two degrees of freedom replacing translation, which are not themselves translation, corre- 
spond to a rotation (as in the very special case of two colliding atoms) or if they are bending fre- 
quencies with force constants of practically zero, then we are justified in applying the simple kinetic 
theory. Even then there will be some difference arising from the fact that (71+72) for transfer of 
momentum is in extreme cases as much as 2.5 times as large as for the corresponding activated 
complex. Thus approximate agreement with simple kinetic theory will occur in ne cases, but 
much lower as well as higher values are to be expected in other reactions. 

In general it does not seem useful to separate our formulas into a collision factor and a sterie factor, 
but if this is to be done we should associate the kinetic theory diameter with the changes occurring 
in the particular six degrees of freedom which correspond to translation before the molecules collide. 
The changes in the other degrees of freedom would then be interpreted as the steric factor. It is 


10H. S. W. Massey and C. B. O. Mohr, Proc. Roy. Soc. A141, 434 (1933); A144, 188 (1934) and subsequent papers. 
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interesting to note that if the two bending frequencies arising from translational terms are stiff 
enough so that the system lies almost entirely in the corresponding lowest states, the reaction diam- 
eter as just defined will be less than the kinetic diameter by the factor 82*JkT/h? which may reach 
value of the order of 100 for heavy atoms and moderately high temperatures. 

The present formulation of the calculation of absolute rates of chemical reactions has certain 
features in common with a number of more intuitive previous treatments," but has more in common 
with the treatment of Pelzer and Wigner. It goes beyond these in formulating the general problem 
in a way susceptible to treatment with our present potential energy surfaces and in pointing out the 
consequences of quantization on the temperature coefficient and the difference in rate for isotopes. 
The fact that the activated complex is much like any other molecule except in the degree of freedom 
in which it is flying to pieces makes possible our comparatively simple formulation. A number of 
investigations are now in progress in which the absolute rate of reaction is being calculated. 

I want particularly to thank Dr. Bryan Topley for valuable discussions as it was with him the 
present calculations of absolute rates were begun. I also want to thank Professors Taylor and Webb 
for helpful discussions. 


1 W. H. Rodebush, J. Chem. Phys. 1, 440 (1933); V. K. La Mer, ibid. 1, 289 (1933); O. KX. Rice and H. Gershino- 
witz, ibid. 2, 853 (1934). 
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Thermodynamic Functions of Tetramethylmethane' 


Louis S. KasseEL,? Pittsburgh Experiment Station, U. S. Bureau of Mines, Pittsburgh, Pennsylvania 
(Received November 20, 1934) 


The Q-sum corresponding to La Coste’s eigenvalues for the rotational levels of tetramethyl- 
methane has been expressed by the usual asymptotic expansion. The relatively large coupling 
terms between the various rotations cancel almost exactly in the Q-sum, so that the thermo- 
dynamic functions are scarcely affected by the coupling. 


A COSTE* has given an exact solution for 

the rotational wave equation of tetramethyl- 
methane with zero potential energy. The energy 
levels are 


E= 
— 2gm,], (1) 


where J is the moment of inertia of the molecule 
with methyl groups fixed, J, the moment of 
a methyl group about the C—C axis, Ip=J 
—(4/3)I, and 


(2/3) (mime 


(2) 


1 Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to oo 

Associate physical chemist, U. S. Bureau of Mines, 
Pittsburgh Station, Pittsburgh, Pa. 

* La Coste, Phys. Rev. 46, 718 (1934). 


The quantum numbers m, m2, m3 and m, of the 
single methyl groups take integral positive and 


’ negative values, while 7, m, and m, are the ordi- 


nary quantum numbers for a symmetrical top. 
The term gq’? within the bracket represents a 
coupling between the various internal rotations 
which raises the energy levels and thus decreases 
the Q-sum. The term —2gm, splits the levels 
without shift in their center of gravity and 
hence raises the Q-sum. Each of these effects 
is quite large, but the two will be found to cancel 
almost exactly. 
With the notation 


(3); p=Lo/Is, (4) 


the Q-sum is 


Q=C 


27 my mM mi me ms m 


— (5) 
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The term e??"¢" may then be expanded as a power 
series, and each term replaced by its mean value 
over the range m,= —j to +j, by the exact re- 
lations 


(6) 
(7) 


m?= (1/3) j(j+1), 

(1/5) (1/15) G+1), 

+(1/21)j(j+1). (8) 


Averages over odd powers naturally vanish. The 
sums over j, m, and m, may now be executed at 
once, using Eqs. (34), (36), (39) and (43) of a 
previous paper.’ At this stage the result is 


Q=L LUD 1 (a — 0/6)? 


+705 

1 — (1/6) 


mi, m2 m3 m4 


+ (7/360) o4g*— (31/15120) 
Xexp [— ]. (9) 


The final sums are quite readily obtained, using 
the asymptotic expansion 


en am? — qil2q-1/2 
—oo 


and its derivatives with respect to a. All sums of 
the type > m?"+1e-«”’ vanish identically. The 


result is 


+ (7/54) — (62/945) +--+]. (11) 


4 Kassel, J. Chem. Phys. 1, 576 (1933). 


(10) 


KASSEL 


Since p=31.54,° the factor in brackets is ap- 
proximately 1—0.01¢+0.0010? and is thus much 
less important than e’/*, the correction term 
which carries over from the ordinary rigid top. 
Both factors may be discarded above 100°K, 
where o = 0.00220. The limiting high-temperature 
value is then 


and, apart from symmetry and nuclear spin 
corrections, 

S= RE(7/2)+1n RT/h?) 7/72 
H=(7/2)RT 


Cp=(7/2)R. 


(13) 
(14) 
(15) 


Also, of course, 
and 


When symmetry and nuclear spin are taken into 
account, the limiting high-temperature value of 
Q is multiplied by 2/34 12=1024/243.5 The 
effect on the asymptotic expansion is not in all 
cases quite this simple and cannot very easily 
be calculated. It is quite certain that the use of 
this factor is justified above 100°K. Deviations 
of the rotational specific heat from (7/2)R will 
then be due to a potential energy associated 
with the internal rotation, and might be used as 
an experimental approach to the problem of 
determining this potential. 

The limiting equation (13) agrees with the 
general equation tentatively proposed by Hal- 
ford’ except for the occurrence of Jy instead of J, 
the slightly large total moment of the rigid 
molecule. For ethane, Halford’s equation agrees 
with the calculation of Mayer, Brunauer and 
Mayer® only if the external moments are taken 
to be A, A and C, where the moments of the 
rigid molecule are A, A and 2C. A correction of 
this nature to the external moments will always 
be necessary, but ethane is perhaps the only 
case in which it is not small. 


5 This is calculated with dogg =1.54X1078 cm and dcx 
= 1.08 X cm. 

6 Cf. Ludloff, Zeits. f. Physik 57, 227 (1929). 

7 Halford, J. Chem. Phys. 2, 694 (1934). 

8 Mayer, Brunauer and Mayer, J. A. C. S. 55, 37 (1934). 
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The Structure of Caesium Enneachlordiarsenite, Cs;As.Cl, 


J. L. Hoarp Anp LEONARD GOLDsTEIN, Department of Chemistry, Stanford University 
(Received November 28, 1934) 


Cs;As2Cly is shown to have a hexagonal unit of structure 
with d9=7.37A, c@ =8.91A, containing one molecule. The 
atomic arrangement is based upon closest-packing of 
caesium and chlorine, with D;?— C 3 2 as the most probable 
space-group, but the amount of distortion from the closest- 
packed positions is not large enough to permit of the un- 
ambiguous elimination of D3u3—C 3 m. A structure based 
upon D;? accounts well for observed intensity data. Arsenic 
is shown definitely by the intensity data to have only a 
pseudo-coordination number of six, and the structure may 


be qualitatively described as containing AsCl; molecules 
embedded in an environment of caesium and chloride ions 
in such a way that all large atoms fit into an approximately 
closest-packed arrangement. This is the first structure of 
this type to be discovered. The ratio of 3Cl per Cs is a 
particularly happy one for structures of this kind, and a 
study of this structure gives much insight into the exis- 
tence and stability of the related group of compounds dis- 
cussed in the Introduction. 


INTRODUCTION 


NUMBER of compounds with the type 

formula A;B2Xy in which A may represent 
caesium or rubidium, B arsenic or antimony, and 
X chlorine, bromine, or iodine were originally 
prepared and described by Wheeler.' These com- 
pounds crystallize with some form of three- 
fold or sixfold symmetry, or in some cases with 
only a close approximation thereto. The crys- 
tallographic data indicate that two or more 
essentially distinct but probably closely related 
structural types are represented. Two points 
seem of particular interest: (1) How can we 
account for the existence and stability of these 
apparently complicated compounds, and (2) 
what sort of bonding and coordination properties 
should we ascribe to arsenic (or antimony)? 
In order to answer these questions, we have made 
an x-ray study of caesium enneachlordiarsenite, 
Cs;AseClo, resulting in a determination of the 
atomic arrangement existing within the crystal. 
We have found from x-ray data that two com- 
pounds of similar formula, Cs3TleCly and 
Cs;Tl,Bro, first prepared by Pratt,? belong to 
quite a different structural type.* 


DETERMINATION OF THE STRUCTURE 


Following the directions of Wheeler, crystals 
of CssAseCly were prepared in the form of thin, 
six-sided, optically uniaxial plates somewhat less 


' Wheeler, Am. J. Sci. 146, 88 (1893); 146, 269. See also 
Groth, Chemische Krystallographie, Engelmann, Leipzig, 
1, 434 (1908). 

* Pratt, Am. J. Sci. 149, 397 (1895). 


Unpublished material. 


than one mm in diameter. Laue and oscillation 
photographs were prepared using, respectively, 
general radiation from a tube with tungsten 
target operated at a peak voltage of about 50 kv, 
and K-series radiation from a tube with water- 
cooled molybdenum target. 

Application of the Bragg formula to the ob- 
served spacings of the equatorial reflections from 
the basal plane appearing upon oscillation photo- 
graphs gives a value for dqpo.1)/m: of 8.910 
+0.020A, spectra from a cleavage face of calcite 
being used as a comparison standard. Application 
of the Polanyi equation, s\=/ sin u, to observed 
spacings of layer lines on these photographs 
leads to identity distances along the axes of 
oscillation of about 7.4A and 12.8A (= ¥3-7.4A) 
accordingly as the axis of crystallographic de- 
scription (using hexagonal axes) or a direction 
making an angle of 30° with it is chosen as the 
axis of oscillation. The data from three com- 
pletely indexed Laue photographs taken with 
the x-ray beam inclined at various small angles 
to c are completely accounted for on the basis 
of a hexagonal unit with a9=7.37A, co=8.91A; 
for no observed reflection gives a calculated value 
of md less than 0.24A, the short wavelength 
limit of the radiation employed. On one photo- 
graph, for example, a total of about 135 spots 
appear with more than 40 forms reflecting in 
the first order ; and in no case is a spot of doubtful 
nature observed, so that this unit may be con- 
fidently accepted as the true one. Reflections 
from (10-0) and (11-0) obtained by transmission 
through the very thin plate-like crystals lead to 
values for a) of about 7.4A. Inasmuch as these 
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reflections were too wide for accurate measure- 
ment, we have preferred to use Wheeler’s ap- 
parently reliable value of 1.209 for co/ao in 
conjunction with our value of 8.91A for co, thus 
leading to a9=7.37A. 

With one molecule of Cs;AseCly within the 
unit cell, the calculated density is 3.42 g/cc, 
somewhat greater than that of methylene iodide. 
As we had only a small amount of material, no 
convenient method of determining the density 
experimentally was available, but inasmuch as 
the calculated value lies within the limits to be 
anticipated, we conclude that we have to deal 
with a one-molecule unit. As we shall show, all 
other evidence supports this conclusion. 

A Laue photograph with the x-ray beam 
accurately perpendicular to the basal plane 
shows no deviation from the threefold axis and 
three planes of symmetry of Ds. A number of 
first-order reflections of the types {mo-1} with / 
odd and {mm-l} with 1 odd appearing on Laue 
photographs together with the observation that 
{00-1} is not thirded eliminates all space-groups 
with special criteria. As ao is found to lie along 
a twofold axis and 30° from a symmetry plane, 
there remain as possible space-groups C;,' 
—C3m, D?—C32, and D3#—C3m, with no 
conclusive evidence at this stage to enable us to 
decide upon one of them. We are able, however, 
to proceed with a structural determination based 
upon convincing premises, and leading to the 
selection of D;? as the probable space-group. 

The radius® of Cl- (for coordination number 
six) is 1.81A. A hexagonal unit of three layers 
made up of 12 spheres of this size in closest- 
packing would have ap>=7.24A, co=8.86A. The 
radius of Cst is 1.69A, but little less than that 
of Cl-, and our unit cell contains 3Cst+ and 9CI-. 
Comparing our lattice constants of a@9=7.37A, 
Co= 8.91A with those above, the conclusion seems 
inescapable that the structure of Cs;As2Cly is 
based upon a nearly closest-packed arrangement 
of caesium and chloride ions. Furthermore, if 
we make the additional and quite reasonable 
requirement that two caesium ions shall not be 
in contact with each other, it is not difficult to 


“See Astbury and Yardley, Tabulated Data for the Ex- 
amination of the 230 Space-groups by Homogeneous X-Rays, 
Phil. Trans. 224, 221-57 (1924). 

5 Linus Pauling, J. Am. Chem. Soc. 49, 765 (1927). 
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see that we are led to a unique distribution 
(excepting small distortions) of the 3Cs+ and 
9CI- within the unit cell. This ideal or un- 
distorted arrangement is made up of three 
identical layers, and may be conveniently sym- 
bolized as follows: Let A represent a closest- 
packed layer with Cs* at 00, B a layer with 
Cs+ at 32, and C a layer with Cst at 24. The 
sequence of layers along c is then ABCABC--.-, 
corresponding to cubic closest-packing. (See 
Figs. 1-3.) The face-centered unit cube, con- 
taining 4 spheres, has Cl- replaced by Cs* at 
the corners. We note that the rhombohedral 
angle® a for Cs3AseCly is 90° 30’, which is indica- 
tive of the small amount of distortion of the 
structure from cubic closest- packing. 

Taking layer A at z=0, we first place ar- 
senic in 43%, 33%, at the centers of octahedra of 
chlorine atoms in positions such that arsenic 
atoms remain far apart one from another, and 
also well-removed from caesium ions. The result- 
ing arrangement is obtainable from each of the 
three space-groups which must be considered, in 
each case with highly specialized values of a 
number of parameters. Moreover, a comparison 
of calculated structure factors for this structure 
with observed intensities for the various orders 
of (00-1), (10-0), and (11-0) make it appear 


Fic. 1. Photograph of a model of the unit cell of Cs;As2Cls, 
with the chloride ions of layer A and all caesium ions shown 
as spheres in closest-packing. Positions of arsenic and of the 
chlorine atoms of layers B and C are indicated as the ver- 
tices of the triangular pyramids, which have the dimensions 
of structure II. 


6 The value of the crystallographic rhombohedral angle « 
. <9 given by Groth, reference 1, as 89° 3’ instea 
of 90° 30’. 
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Fic. 2. Packing drawing of three layers parallel to (00-1) 
showing the undistorted closest-packed structure made up 
of spheres of the same size throughout. Caesium may be 
distinguished from chlorine atoms by the type of shading. 


INK 
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Fic. 3. Packing drawing of the unit cell normal to (10-0) 
with caesium and chlorine atoms shown as spheres in 
closest-packing. The positions of arsenic are indicated by 
the apices of the pyramids. 


altogether likely that a suitable distortion of the 
ideal structure will account for all of the x-ray 
data. 

Each chlorine of layer 4 is shared between two 
AsCl, octahedra, while each chlorine of layers 
B and C is a part of but one such group. Regard- 
ing the structure as a purely ionic aggregate of 
Cst+, Ast+#, and Cl-, we see that the strength of 
the electrostatic valence (including the con- 
tribution of Cs+) reaching each Cl- of layer A 
is 23, while to each Cl- of layers B and C it is 12. 
Actually we expect the arsenic to have a strong 
tendency to form three normal covalent bonds 
with chlorine, and the ionic model shows clearly 
that such bonds would be formed between 
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arsenic and the chlorines of layers B and C; 
for the formation of such bonds means that 
arsenic layers will be found somewhat closer to 
layers B and C than they are to layer A. It 
appears also that the twofold axes are to be 
retained and that the actual structure is based 
upon D;? or rather than C3,'. 

With either D;? or D34*, caesium and arsenic 
are placed in the following positions: Cs in 
000, 43 u, 33 a, As in 332, 34 3. We expect 
u =}, 3<v< Moreover, the complete structure 
factor including chlorine for (00-/) reflections is 
identical for the two space-groups, namely 


S{ 00.1) = fes(1+2 cos 2mlu) 
+2f,4, cos 2xlv+3f¢:(1+2 cos 27/z), 


in which z=}. We have intensity data for the 
first twelve orders of reflections from (00-1) 
with which to compare theoretical values for the 
proposed structure as calculated from the for- 


mula 
1+cos? 


I=C 
2 sin 26 


No corrections are included for extinction or for 
the decrease in intensity arising from the thermal 
vibrations of the atoms within the lattice. The 
f-values of Pauling and Sherman’ are employed 
throughout this paper. The factor / is required 
since we are considering reflections at small 
angles from a thin plate-like but strongly absorb- 
ing crystal which was at all times completely 
bathed in radiation ; the effective cross section of 
the beam was proportional then to sin @ and 
therefore to /. Estimating the mass absorption 
coefficient as about 40 for Mo Ka-radiation, we 
calculate that there would be only a quite 
negligible transmission except, perhaps, for the 
last two or three orders. Experimental intensities 
were estimated by visual comparison of re- 
flections from a number of photographs with 
widely varying times of exposure with an in- 
tensity scale prepared for this use. The reflections 
become more diffuse as the order increases, and, 
as our considerations have to do with integrated 
intensities, it may be that we have somewhat 
underestimated the intensities of the last two 
or three reflections. 


7 Pauling and Sherman, Zeits. f. Krist. 81, 1 (1932). 
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With the values of the parameters corre- 
sponding to closest-packing of Cst and Cl-, 
neither caesium nor chlorine would contribute to 
reflections with //3 not integral. We do not 
find it possible to account for such reflections by 
considering them as due to arsenic alone, but 
rather that with 0.795 <v<0.815 and with both 
u and z slightly less than 3, a general agreement 
is obtained between calculated and observed 
intensities. Table I contains these data for 


TABLE I. Comparison of calculated with observed intensities 
for (00-1) reflections. 
v=0.805, 


u=0.320, s=0.325 


Observed Calculated 
intensity* intensity* 


5 
28 
93 
41 

4 
32 


Observed Calculated 
intensity intensity* l 


* Temperature factor not included. 


v=0.805, u=0.320, z=0.325, which values seem 
to be most satisfactory. As 3f¢: is nearly equal 
to fcs, we see that the contributions of these 
two kinds of atoms are nearly identical in case 
u=z. It seems improbable that either wu or z is 
less than 0.315, or that both together can be 
appreciably less than 0.320. However, it seems 
most satisfactory to have both uw and z some- 
what less than 4. Making the caesium parameter 
slightly less than that for chlorine appears 
reasonable from a consideration of dimensional 
relations, but is otherwise quite arbitrary. 

We now proceed to compare calculated in- 
tensities for (40-0) and (hh-0) reflections with 
the experimental data. These data are the result 
of transmission photographs taken through the 
thin plates, which, however, were of such small 
diameter as to be completely bathed in radiation 
at all times. The intensity formula employed was 
identical in form with that already given except 
for the omission of the factor J. In this case the 
correction for absorption increases with the order, 
but calculation shows that the increase is not 
large enough at most to be very important. 

The form of the structure factor depends upon 
the space-group of the structure. If it be based 
upon D;,, all chlorine atoms lie in planes 
of symmetry and the structure has but one 
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parameter at right angles to c; in particular, all 
atoms of layer A occupy no-parameter positions, 
and each caesium ion is surrounded by six 
equidistant chloride ions. If on the other hand 
the true space-group be D,?, the planes of 
symmetry are no longer retained, and three of 
these six chlorides could be nearer than the others 
to each caesium ion. This second possibility 
seems to us the more probable, especially in 
view of the fact that Rb;AseCl, apparently is 
isomorphous with our compound, and the con- 
siderably smaller radius of 1.48A for Rb* as 
compared with 1.69A for Cs+ and 1.81A for Cl- 
would seem to require such a closer anion- 
cation approach. Six chlorine atoms would then 
occupy the general positions of D;’, and the 
chlorines of all layers could be rotated somewhat 
about the threefold axes away from the closest- 
packed positions for identical spheres, with 
probably a better distribution of structural 
strains than is obtainable from D;,°. We have 
then, 
3 Cl in w00; 0w0; @HO with w=0.480, 
6 Cl in x, y, 2; y—x, %, 2; 9, x—y, 2; 
X, y—xX, Z;x—y, 2, 

with x=}, z=0.325. If the structure were 
based on x=2y, and w=}. 

The amounts of the deviations of x, y and w 
from the closest-packed positions seems to be 
too small to have any large effects upon the 
calculated intensities. In Table II we give calcu- 
TABLE II. Comparison of calculated with observed intensities 

for (hO-0) and (hh-0) reflections. 
Structure I: x=}, y= %. Structure II: x=0.365, y=0.200. 


Calculated intensity* 


Observed Structure Structure 
intensity I II 


dR 
ow 


* Temperature factor not included. 


lated intensities for structure I for which «=}, 
y=, w=} as well as for structure II based 
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upon D;? for which x= 0.365, y= 0.200, w= 0.480. 
The agreement is generally satisfactory in both 
cases, although somewhat better for the second. 
The data indicate that there is indeed some dis- 
tortion from structure I, and that structure II 
is one variation which meets the intensity re- 
quirements. Table III contains some intensity 


TasBLe IL]. Comparison for structure II of calculated with 
observed intensities for (Ok-1) reflections. ‘ 


Observed Calculated 
intensity intensity* 


120 140 
0.3 
2 
12 
12 
125 
55 
1.5 
12 
10 
0.3 
50 
$3 
17 


Observed Calculated 
intensity intensity* 


3.0 
47 
1.5 
0.1 
80 
20 
0.5 
5.5 
13 
5.5 
5.5 
0.5 
12 
9.2 
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01- 
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* Temperature factor not included. 


comparisons which are representative of a more 
general type of reflection. The data are given 
only for structure II, but again the comparisons 
are appreciably better than for the undistorted 
structure. That we have arrived at essentially 
the correct structural type is scarcely to be 
doubted, but it is primarily upon the basis of 
dimensional and structural rather than intensity 
considerations that we decide that distortions 
from the ideal structure are to be introduced in 
accordance with D,? rather than D;,’. 


DISCUSSION OF THE STRUCTURE 


A photograph of a model of the structure is 
reproduced in Fig. 1. In this model the chloride 
ions of layer A and all caesium ions are shown 
as spheres of equal size in closest-packing, while 
the positions of arsenic atoms and of the chlorine 
atoms of layers B and C are indicated as the 
vertices of the triangular pyramids. Figs. 2 and 3 
show packing drawings of the unit cell, again 
with caesium and chlorine atoms indicated as 
spheres of equal size in true closest-packed 
array, and with the positions of arsenic marked 
by the apices of the pyramids. Study of these 
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drawings indicates the essential nature of the 
structure, and the process of going over to the 
parameters of structure II may then be visualized 
quite easily. 

There are two factors which show that the 
plane z=} between layers B and C should be a 
natural cleavage plane: (1) the absence of arsenic 
atoms in this region, and (2) the fact that both 
chlorine and caesium parameters along c are 
slightly less than }. 

Since the basal plane intensities require v 
=0.805+0.010 rather than 2, we see, as is 
indicated in the various figures, that arsenic 
trichloride molecules would seem to have their 
identity preserved within the structure, and that 
the process of crystallization has consisted in 
embedding AsCl; in an environment of com- 
pletely ionized caesium chloride in such a 
manner that the large atoms or ions are in nearly 
closest-packing. The ratio of 3Cl per Cs, leading 
to the kind of layers already described, is par- 
ticularly well suited to the requirements of this 
general type of structure. This ratio is also 
shown by the widely occurring perovskite or 
KMgF; structure. 

We believe, moreover, that the size and shape 
of the AsCl; molecules are probably described 
well by structure II. Arsenic is expected to form 
three normal covalent bonds at angles somewhat 
greater than 90° but less than the tetrahedral 
angle.’ Thus in elementary arsenic three bonds 
are formed at 97°, while electron diffraction 


TABLE IV. Interatomic distances calculated from structures 
I and II, (The subscripts refer to the layers.) 


Separation 


Separation 
structure I] 


Neighbors structure | 
3Cle 
3Cla 
3Cla 
C3lz, 3Cle 


8 Linus Pauling, J. Am. Chem. Soc. 53, 1367 (1931). 


O1 
01 
01 
Atom 
As, 2.25A 
As, 2.76 
Cs4 3.54 
3.83 
Csa 3.71 
Csp 3Cl4 3.56 
Csp 3Cle 3.69 
Csp 3.57 
Csp 3.83 
Cls 4Cl, 3.70 
Cla Cla, Cle 3.55 
Cla Cla, Cle 3.55 
Cle 2Clp 3.35 
Cle 2Clg 4.04 
Cle Cle 3.78 
Cle Cle 4.08 
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studies® of AsCl; lead to 101+4° bond angles. 
An As—Cl separation of 2.18+0.03A is given 
by the electron diffraction study, while a value 
of 2.20A is obtained from Pauling and Huggins’ 
covalent radii.!° Structure II, which has v= 0.805, 
x=0.365, y=0.200, 2:=0.325, contains AsCl; 
molecules with bond angle 95°50’, As—Cl 
=2.25A, Cl—Cl=3.34A. Table IV lists all of 
the pertinent interatomic distances for the two 
structures for which we have calculated in- 
tensities. The As—Cl separation of 2.25A given 
by structure II is nearly that expected for a 
covalent bond, whereas the precise significance 


® Brockway and Wall, J. Am. Chem. Soc. 56, 2373 (1934). 
10 Pauling and Huggins, Zeits. f. Krist. 87, 205 (1934). 
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of the corresponding value of 2.42A given by 
structure I would remain doubtful. Moreover 
the relations of caesium with chlorine are 
probably more satisfactory in general for struct- 
ure II than for structure I. 

Possible variations in the parameter values 
would include having z about 0.005 smaller, 
v from 0.005—0.010 smaller, and x and y a little 
nearer the closest-packed positions. Such changes 
would not materially alter our concept of co- 
valent As—Cl bonds, but would lead to bond 
angles a little greater than 96°. They would 
result also in a somewhat smaller value than 
3.71A for the Cs, — Clg distance, which is perhaps 
to be anticipated. 
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A wave function in the form of a series with coefficients 
to be determined by the variational method is used to 
calculate the dissociation energies of the hydrogen molecule 
in two of its excited states: and 1se2po 'Z,*. 
In the former case an 11-term wave function gives an energy 
less than 0.02 e.v. from the most probable experimental 
value and lying within the limits of experimental error. 
For the second state an energy less than 0.08 e.v. from 


INTRODUCTION 


HE treatment by wave mechanics of the 
electronic states of a diatomic molecule has 
led in general to results which are unsatisfactory 
from a quantitative point of view. The very 
simplest problems, including the excited states 
of the hydrogen molecule, have not been carried 
to the point where quantitative agreement with 
experiment could be obtained. The successful 
treatment of the normal state of the hydrogen 
molecule by James and Coolidge! furnishes a 
method whose possibilities it seems worth while 
to explore more thoroughly; and with this in 
mind a calculation has been made on several 
of the two quantum states of hydrogen,? and the 
1 James and Coolidge, J. Chem. Phys. 1, 825 (1933). 
2 Two-quantum molecular states dissociate upon adia- 


batic separation of the nuclei into one normal atom and one 
atom in a two-quantum state. 
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experiment is obtained by means of a function comprising 
18 terms. The interelectronic distance has been introduced 
explicitly into the wave function and the role that it plays 
in reducing the energy is discussed at length. The limits of 
accuracy of the calculation with reference to convergence 
and nuclear motion are considered. Explicit wave functions 
are given for the two states. 


results for two of these states are reported in this 
paper.’ 

Previous work on these states has been based 
on trial wave functions that have been con- 
structed from hydrogen atom functions?‘ (cf. the 
Heitler-London method), from the wave func- 
tions for the hydrogen molecular ion,® or from a 
combination of the two.* In a molecule as closely 


3 Work has also been done on the repulsive state 1so2pc 
3y,,* and will be published shortly in a joint paper by Drs. 
James and Coolidge with the writer. 

4 Kemble and Zener, Phys. Rev. 33, 512 (1929); Zener 
and Guillemin, Phys. Rev. 34, 999 (1929). 

5E. A. Hylleraas, Zeits. f. Physik 71, 739 (1931). 
Hylleraas fails to treat the state 1s02sa *2,*+, but mentions 
results obtained by a rough extrapolation which appears to 
have little theoretical justification. 

6 J. K. L. MacDonald, Proc. Roy. Soc. A136, 528 (1932). 
MacDonald treats all of the two-quantum states for the 
same internuclear distance, 2 Bohr radii. Both Hylleraas 
and MacDonald obtain fairly good agreement with exper!- 
ment on the II states, but their results on the = states are 
unsatisfactory. 
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bound as hydrogen it is to be expected that the 
molecular ion function of Hylleraas and the 
hybrid function of MacDonald will be more 
effective than a Heitler-London function, as is 
confirmed by their results. However these wave 
functions, together with all others previously 
employed, suffer a twofold disadvantage which 
renders them unsatisfactory if accurate results 
are desired: (1) they are insufficiently flexible, 
if flexible at all; (2) they take no adequate ac- 
count of electronic interaction; e.g., for 2 states 
they are functions of only four of the five es- 
sential coordinates, the interelectronic distance 
or its equivalent having been omitted. 

The importance of taking account of elec- 
tronic interaction in the wave function is to 
ensure that there will be a small probability for 
those regions of configuration space in which the 
electrons are close together. For example, the 
introduction of positive powers of the inter- 
electronic distance (ri2-terms) into the wave 
function will serve to keep the electrons apart, 
and where that service is partly accounted for 
by other means, 71:2-terms will be less important. 
Thus 7i2-terms should be less important for 
triplet than for singlet states, since the wave 
function for a triplet state is antisymmetric in 
the space coordinates of the electrons and will 
vanish when they coincide. Also the explicit 
appearance of azimuthal angles in the wave 
function for a II state will compensate in some 
measure for the absence of 7;2-terms, if the state 
happens to be a triplet. In consequence we should 
expect that functions of the type previously 
used would have given somewhat better agree- 
ment for triplet than for singlet states and better 
for *II than for *2. This appears to be in agree- 
ment with the results of MacDonald, although 
the former effect is not so marked as one might 
suppose. A more detailed discussion of some of 
these points is reserved till later when a compari- 
son can be made with the results of the present 
work. 


METHOD 


The observed two-quantum © states are 
1sc2so *Z,+ (the most probable upper state for 
the continuous spectrum), 1so2pc1Z,+ (the B 
state) and (the Heitler-London 
repulsive state). The symmetry characteristics 
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impose twe types of conditions on the wave 
functions: 

(1) Symmetry or antisymmetry to interchange 
of the space coordinates of the electrons (singlet 
or triplet states). 

(2) Symmetry or antisymmetry to interchange 
of the space coordinates of the nuclei (even (g) or 
odd (u) positive electronic states).? Furthermore 
the wave function of a = state must have rota- 
tional symmetry about the nuclear axis and can 
contain only the difference between the azimuthal 
angles of the electrons. 

It is proposed to expand the wave function 
in terms of a power series in the five variables 
with an exponential factor in the elliptical co- 
ordinate \ to ensure convergence of the integrals 
involved. If the nuclei are designated by a and b 
and the electrons by 1 and 2, and the distances 
between the particles expressed in Bohr radii 
(B.R.) are denoted by the customary symbols, 
then the coordinates are: 


(fiatrw)/R, 
p=2ri2/R. 
The wave function is: | 


mnjkp 


where [ménéejkp | stands for 


( 1/2 1) { 


No spin function is introduced, but the spin is 
determined implicitly by the symmetry char- 
acteristics of the wave function y in accordance 
with Pauli’s principle. For singlet states the 
plus sign is to be used, for triplet states the 
minus sign. For the even (g) states 7+ must be 
even, for the odd (u) states odd. Subject to this 
last restriction the summation is to be extended 
over positive or zero values of the indices and 
as many terms taken as result in an effective 
lowering of the energy. The quantities 5; and 62 
should preferably be taken as arbitrary param- 
eters to be varied, but since there appeared to 


symmetric 
antisymmetric 
even(g) 


odd() 


7 + terms are in the space coordinates of 


the nuclei if they are 
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be no simple way of performing this variation, 
suitable values were selected at the outset of the 
calculation by a preliminary computation de- 
scribed below. A small error in the choice of 6; 
and 6: should be compensated for by a slight 
alteration in the coefficients of terms involving 
powers of \; and dz. The best values of these 
constants will depend on the internuclear dis- 
tance R, and a particular choice of 6’s will give 
reasonably good results over a small range in R.! 
The electronic energy plotted as a function of R 
will be the envelope of a family of curves char- 
acterized by different pairs of values of the param- 
eters 6; and 62 such that each pair of values 
will be the best for that particular R at which its 
curve touches the envelope. 

The variational method was applied to de- 
termine a set of values for the parameters 
Cmnjkp Which would give the lowest energy ob- 
tainable with the type of function used. Since 
only the lowest states of a particular symmetry 
are considered, the computed value of the energy 
for any trial function must always lie above the 
true energy. The actual computation falls 
naturally into three parts. In the first, extensive 
tables of fundamental integrals, related by 
convenient recursion formulas, are prepared for 
use in evaluating the matrix elements of the 
successive approximations. The principal labor 
occurs in the integrations over 1/712, which is 
expanded by the Neumann formula, and in the 
evaluation of integrals between 7j.-terms. In 
the second part the matrix elements correspond- 
ing to combinations of the simplest terms in 
the series, are evaluated and introduced into the 
secular equation. If the parameters 6; and 62 are 
approximately correct, it is found that the series 
converges rapidly and that only low powers 
(0, 1, 2 but never more than 3) are required. 
Since in the mth approximation therearen(n+1)/2 
matrix elements to be computed and a secular 
equation of the mth degree to be solved, it is not 
feasible to test every term in combination with 
all the others nor to carry along any terms 
which give inappreciable contributions. As the 
importance of any term was generally found to 
decrease as the number of terms combined with 
it was increased, terms which yielded negligible 
improvement at any stage of the calculation 
were usually dropped. However, exceptions have 
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been found to this rule, notably in the case of 
certain pairs of functions, e.g. [06,06.000] and 
[16,16,000 ] in treating *2,+ which gave in com- 
bination a better contribution to the energy than 
the sum of their individual contributions when 
inserted independently. This phenomenon, which 
is fortunately infrequent in occurrence, seems 
to indicate that the 6’s have not exactly the right 
values to ensure the closest convergence. 

The third part of the work consists in the 
solution of the secular equations which arise in 
each approximation. The procedure is neces- 
sarily one of trial and error and involves the 
evaluation of large numbers of symmetric 
determinants.' Here there were no neighboring 
electronic states of the same symmetry, i.e., no 
nearby roots of the secular equation, and it was 
found sufficiently accurate to interpolate linearly 
for \, the value of the root. The coefficients of 
the wave function are readily obtained from a 
secular determinant evaluated at the _ root. 
Before working out this determinant it is de- 
sirable to rearrange the terms by putting an 
important function at the end to ensure that the 
last factor of the determinant will reduce to a 
vanishingly small quantity. 


THE STATE 1so2so *Z,* 


In order to make use of tables already pre- 
pared! in treating the normal state of the mole- 
cule, computations were begun on the state 
1so2so0*Z,+ with a wave function for which 
5:=62=0.75. It was soon found that the con- 
vergence was extremely poor, and this was at- 
tributed to an unfavorable choice of shielding 
constants. Thereupon a preliminary investiga- 
tion for determining a suitable pair of 6’s was 
made upon the single term [161062000 ] which in 
the previous work had proved to be the most 
important single term of the series; this led to 
the selection of 6:=0.5 62:=1 for treating this 
state of the molecule. The calculation was then 
carried out for the internuclear distance R= 1.87 
B.R. at which occurs the minimum of the ex- 
perimentally determined potential energy curve.* 


8 [t would have been more convincing to have established 
this also as the theoretical minimum. However it is well 
known that the experimentally determined Morse curves 
are quantitatively accurate in the neighborhood of the 
minimum. 
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QUANTUM STATES OF Hz 


The six terms [10020 ] [10200] [10000] [00000 ] 
[10110] [11000] were found to give an energy 
much better than any previously reported,’® less 
than 0.08 e.v. from experiment. Upon the in- 
sertion of 712 terms the energy was still further 
reduced and the final wave function involving 
[10020] [10200] [10000] [00000] [10110] 
[11000] [11020] [10001] [00001] [10021] 
[11001 ] gave a total energy of —0.7361 in units 
of twice the ionization energy of the normal 
hydrogen atom: 2R,jic (27.074 e.v.). A number of 
other terms were investigated at various stages of 
the calculation including: [20000]* [01000] 
[11200] [10130] [12000] [11110] [10111]; 
these were all found to yield 0.001 e.v. or less. 
The 712 terms contributed 0.057 e.v. in marked 
contrast to the improvement of 0.44’ e.v. ob- 
served in the ground state (singlet) by James and 
Coolidge; this difference is easily understandable 
from the orbital antisymmetry of the triplet 
state. Because of the extremely rapid conver- 
gence it was possible to make an estimate of the 
probable contribution of remaining terms, from 
which the -absolute minimum may be said 
to occur at A= — 0.7367 +0.0004(— 19.945+0.01 
e.v.). Subtracting the zero point energy of 0.162 
e.v. and the energy — 16.921 e.v. of a normal and 
a two-quantum hydrogen atom (with fixed 
nuclei), we obtain 2.862+0.01 e.v. for the disso- 
ciation energy of the hydrogen molecule in the 
state 1so2so*D,*+. The experimental value is 
given by Richardson as 2.867 e.v., by Mulliken 
as 2.86 e.v. and by Jevons as 2.85 e.v." 

The unusual degree of accuracy which has been 
obtained for this state makes it necessary for 
us to consider closely the approximation which 
has been made in the customary assumption of 
“clamped nuclei.”” The correction to the energy 
for © states may be roughly estimated by means 
of the limiting case of the united atom. The 
electronic energy” is then effectively raised by 
an amount:!*: 14 
mbindsjkp | is abbreviated to [mnjkp] for 6,=0.5, 


10 MacDonald’s result is 0.40 e.v. from experiment. 
1Q. W. Richardson, Molecular Hydrogen and Its Spec- 
trum, p. 324; R.S. Mulliken, Rev. Mod. Phys. 4, 78 (1932); 
W — Report on Band Spectra of Diatomic Molecules, 
p. 


- “Electronic” energy includes energy of nuclear re- 
pulsion. 

L. Kronig, Band Spectra and Molecular Structure, 
pp. 

*R. de L. Kronig, Physica 1, 7 (1934). 
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pL(L+1) ? 
f (—) | 
spin OR 
where yu is the reduced mass of the molecule, 
L the resultant orbital angular momentum of the 
united atom and ® the electronic part of the 
wave function. Since this expression is always 
positive it appears that the uncorrected elec- 
tronic energy may lie below the experimental 
value. A consideration of the effect of nuclear 
motion in Hylleraas’ very exact calculation on 
parahelium® led to a similar result. For the state 
under consideration *2,+) the first term 
in U will vanish. It would seem to be a matter of 
considerable difficulty to appraise the second 
term; fortunately this has been done by Mac- 
Donald® who finds it to be worth about 0.01 e.v. 
The minimum of the corrected potential energy 
curve will therefore occur at A= —19.935+0.01 
e.v. However, since corrections have now been 
made for nuclear motion, we must subtract the 
true value of the two-quantum asymptote 
(—16.912) and this leaves us with nearly the 
same dissociation energy as before: 2.861+0.01 
e.v. The approximations made in assuming the 
minimum of the experimentally-determined po- 
tential energy curve to be the true minimum, 
in estimating the convergence of the series wave 
function and in correcting for the effect of nu- 
clear motion take us into a region of accuracy 
where spin and relativistic effects can no longer 
be neglected. 


THE STATE 1so2po 


Next to be treated was the B state 1sa2po 'Z,,*. 
In order to secure the best results, it would have 
been advisable to have adopted the same pro- 
cedure for determining 6; and 42 as that employed 
in the preceding work, but this would have 
necessitated the entire recalculation of all the 
tables of fundamental integrals from the very 
beginning. The calculation was therefore carried 
out with 6,=0.5, 6.=1 at the Morse curve 
minimum R=2.42 B.R. A great many terms 
were tried and it was found that few could be 
rejected. A wave function comprising the terms 
[10010] [10100] [11100] [00100] [11010] 
[10120] [20100] [01100] [00300] [02100] 
[01120] [01300] [11120] [00101] [10101] 


15 E, A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
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[11101] [01101] [10121] gave the best energy 
of —0.7522 (—20.37 e.v.) which is approxi- 
mately 0.08 e.v. from the experimental value 
given by Mulliken.’ The addition of 712 terms 
to the wave function lowered the energy by 
0.29 e.v. The slowness of convergence is un- 
doubtedly to be attributed primarily to the fact 
that 6; and 6, have not been given their best 
values. Since the r12 terms are quite important 
and since they were inserted last in order, it 
might be thought that if they had been put in 
somewhat earlier, several of the non-rj2 terms 
would have proved to be less important and 
might even have been dropped altogether. The 
contribution of some of these terms was notice- 
ably decreased, but it appears that none of the 
18 could have been omitted without sacrificing 
0.01 e.v. or more, e.g., the best 15-term wave 
function, whose coefficients are given in Table I, 
gave an energy worse by 0.06 e.v. In estimating 
the convergence of remaining terms, the greatest 
error appears to be due to the absence of those 
with higher powers of 712. In the treatment of 
the normal state these terms proved to be worth 
about one-sixth of the energy contributed by 
the first power terms. Allowing from 0.04 to 
0.06 e.v. for this purpose and, upon examina- 
tion, from 0.01 to 0.05 e.v. for convergence of 
other terms, one obtains 3.44+0.03 e.v. for the 
dissociation energy of the B state, as compared 
with the values 3.42, 3.44 and 3.47 e.v. given by 
Richardson, Mulliken and Jevons, respectively. 
This estimate is not quite so reliable as that made 
on the state *Z,*. 

Relative to the total energy, the 712 terms in 
the B state are of the same order of importance 
as in the ground state. Thus it is clear that the 
greater importance of the ri2 terms for these 
two states is to be attributed almost entirely 
to the fact that they are both singlets. One 
might suppose that there would be less overlap- 
ping of the electron orbits in a two-quantum 
state than in the normal state, and that 712. terms 
would be correspondingly less important in the 
former case. Thus Hylleraas and Undheim! 
find that for two-quantum para-helium the rie 
terms contribute only one-tenth as much to the 


16 The energies calculated by Hylleraas and MacDonald 


are, respectively, 0.25 and 0.50 e.v. from experiment. 
17 Hylleraas and Undheim, Zeits. f. Physik 65, 759 (1930). 
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energy as in the ground state. The situation in 
the molecule is sufficiently different from that in 
the atom to make understandable the fact that 
the 712 terms do not drop off in importance in 
the molecule to the extent that they do in the 
atom as one passes from the ground state to the 
excited levels.'® 


CONCLUDING REMARKS 


The best wave function for the state *Z,+ and 
the best 15-term wave function for the B state 
are given in Table I. A recalculation of the latter 
state with a specially selected pair of 6’s would 
probably reduce the number of terms consider- 
ably and provide a somewhat lower value of the 
energy. A preliminary examination of the term 
[06,062100] for a range of values of 6; and 6, 
should enable one to select a suitable pair, or 
else one might determine 6; and 52 by requiring 
a simpler function to approximate the function 
given in Table I or a certain number of its more 
important terms. Tables of matrix elements for 
combinations of numerous terms have been pre- 


TABLE I. Coefficients of the normalized* wave functions 
(Cmnjkp)- 


= —0.73609, R=1.87) *+(A= —0.74986, R=2.42) 


1.441073 
—0.231754 
0.002751 
—0.560549 
0.007094 
—0.009393 
0.019232 
0.207040 
0.077989 
—0.106513 
0.096702 
0.210583 
—0.016744 
—0.020524 
—0.014764 


— 1.161959 
0.499368 
— 0.077834 
—0.014422 
0.071975 
0.034815 
—0.018709 
0.091777 
—C.004742 
— 0.000503 
0.009030 


10121] 


* Note the factor 1/2 in the wave function. 


18 In the two-quantum atom the inner electron will be 
moving rapidly with respect to the outer electron in a 
smaller orbit. Classically speaking, its motion could then 
be represented by a spherical cloud of charge shielding the 
nucleus (the shielding for two-quantum para-helium is 80 
percent effective). The interaction of the outer electron 
with this charge cloud is a function only of its distance from 
the center, i.e., the nucleus, so that 712-terms will not be 
needed to. describe the motion. On the other hand in the 
excited molecule the electrons will have lost their one-and 
two-quantum atomic characteristics to a large extent (all 
the more so for a closely bound molecule such as hydrogen), 
and will move with more nearly equal velocities in much 
the same region of space between the nuclei. Thus 7:2-terms 
will be nearly as important as in the ground state. 
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pared for treating both even (g) and odd (x) 
states. These tables, although unpublished, are 
accessible and may be used for calculating any of 
the 3-quantum = states or the unobserved 
state 1so2so'1X,+. These states should appear 
as higher roots of the secular determinants, for 
to each higher root there will correspond a state 
of the system such that the uth lowest root 
will be an upper limit for the energy of the mth 
lowest state of the given symmetry."” 

In possible extensions of the James-Coolidge 
method to other molecules, the present work 
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would indicate that where it becomes a practical 
impossibility to take into account 712 terms or 
their equivalent, quantitative results for = 
states can be expected only if the states are 
antisymmetrical in the space coordinates of all 
the electrons.'® 

In conclusion I wish to thank Dr. H. M. 
James and Dr. A. S. Coolidge for the use of their 
tables, and I am especially indebted to Dr. 
James for valuable advice and encouragement. 
I am also obliged to Professor Kemble for help- 
ful comment on the manuscript. 


The procedure adopted for evaluating the matrix elements follows closely that described in the 
appendix of James and Coolidge, and will consequently not be reproduced in detail. Modifications 
of most of the formulas are needed because the previous work was restricted to the special case 


j= bo. 


The required quantities can all be found in terms of the integrals: 


1 
X f f f f f f (Aa? — poy?) 9" gid 


1 
X f f f f f f (Aa? — pr?) 28241-28102) gid go, 


x f f f ff f (Ag? — poy?) 9 gid go, 


where a2=252, a= 


The indices are all positive or zero except p which may have the value —1. 
Consider the matrix component of unity between the two terms [m.6\Mabejakapa] and 


[m51246ej0keps ]. Since each term is the sum or difference of two parts which we may designate 
Va’, Wa’’ and Wo’, o’’, respectively, the matrix component of unity will have four parts: 


Sa= f f ret f f Vel f f f 


where the plus sign is for singlet and the minus sign for triplet states. Remembering that 


R 
dridt2= 1? — 1”) (Ao? — we”) dA gid 


we find: 


Ma tM t2, a2, jot jo, kate, pat po) 
—X(matm, a1, Nats, Jato, katk,+2, Pat pr) 
a, natms+2, a, Ratjo, Pat pr) 


—X(matm, a, Natm, Qa, jath, Ratjot2, Pat fr) 


"Or at least in the space coordinates of the binding electrons. 
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The matrix component of the Hamiltonian energy will similarly consist of four parts: 
Ro jor Ro joPo “Myb jokoPo 
The J’s may be evaluated in terms of the X’s by means of the following general formula: 
I (” 
byedysfhpr 
—4{(d?+d—h?—h—-,4?)X (00000) — 2y4(1+d)X (01000) + y42X (02000) + 2y4dX (0 — 1000) 
—d(d—1)X(0— 2000) +A(h—1)X (000 — 20) } —4p,{ (d—h)X (00000) — y4X (01000) 
+(d+h+ p,+1)£X (0200-2) — X (0002 — 2) ]+ (h—d)X (2000 — 2) +X (0020 — 2) ] 
(2100 — 2) — X(0300— 2) +X (0120-2) +X (0102 —2) —2X(1011—2) ] 
+2dX(1—111—2) —2hX(111—1-—2)}}. 


R4 
) R{ X (02000) — X (00020) — 8X (01000) +2.X (0200— 1) — 2X (0002 — 1)} 


The abbreviation X(rstuv) has here been used for the integral 
X(atb+r, yitye, C+dts, et+ftt, gthtu, patpt?). 


Fortunately, by virtue of the unsymmetrical nature of this formula, the bracket multiplied by p, 
will drop out for all combinations of non-r,z terms and for all combinations of an 72-term with a 
non-7j2 term. 

The X’s may all be evaluated in terms of generalized Z’s by the following formulae: 


X (manarjkp) = Z°(m +2, ai, n, a2, j, k, p) —Z°(m, ar, n, a2, k, p) 
1 
M?= (di? — 1) (Ao? — 1) (1 — i?) (1 — 


For the calculation of the Z’s there are needed the generalized integrals H,’(maynaz) and S(manaz), 
already familiar in molecular problems for the special case a:= az. The details of these generalizations 
are more or less obvious and the relevant formulae” and recurrence relations may be obtained from 
the appendix of the James-Coolidge paper if one makes the following alteration in their formulae: 
replace the single a by a and a2 in such a way that m always goes with a, and m with az. For ex- 
ample: 


Z(mayna2jkO) = 4A m(a1)A n(a2)/(J+1)(k+1), 
Z(maynarjk—1)= (27+ 1)R,(7) RAR) 
7=0 
H,(manaz) = Ho(m+1, a1, n+1, a2) — S(m, a1, n+1, ae) — S(n, Qe, m-+1, a). 


20 Some of these formulae are contained in a recent paper by H. M. James, J. Chem. Phys. 2, 794 (1934). 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


The Fractionation of the Oxygen Isotopes in an 
Exchange Reaction 


Recently Urey and Greiff! have presented theoretical 
calculations on the equilibrium constants of exchange re- 
actions involving the isotopes of the lighter elements and 
among others the exchange reaction, 


2H,0'* +CO,"* = 2H,0* 


At 0°C the equilibrium constant for this reaction is found 
to be 1.097, which gives for the fractionation factor for this 
reaction 1.047, the O'* being concentrated in the carbon 
dioxide. This reaction is well suited to the separation of 
the oxygen isotopes by using counter-current scrubbing 
apparatus similar to fractionation columns in construction. 

In this note we wish to report the experimental con- 
firmation of the fractionation factor calculated from theory. 
The method consisted in passing carbon dioxide gas from a 
tank through three spiral wash bottles in series containing 
water from which the deuterium had been largely removed. 
The gas was then mixed with deuterium-free hydrogen gas 
prepared electrolytically and the mixture passed over a 
supported nickel catalyst to convert the carbon dioxide to 
methane and water. The specific gravity of the water in 
the wash bottles and the water produced catalytically from 
the carbon dioxide were compared. The deuterium-free 
water used in these experiments was prepared by electro- 
lyzing a large volume of ordinary water, collecting the 
first 10 percent electrolyzed, then re-electrolyzing this and 
catching the first 10 percent re-electrolyzed. The fractiona- 
tion factor of these cells runs between 9 and 10 and hence 
most of the deuterium should have been removed. The 
specific gravity of the water at 25° was 0.999979, de- 
termined with a pycnometer.* The specific gravities of the 
water from the wash bottles and from the carbon dioxide 
were compared by a pressure float method described by 
Gilfillan.2 The determination of the specific gravity of the 
water from the wash bottles by this method agreed with 
the pycnometer value given above. Three determinations 
of the specific gravity of the water from the carbon dioxide 
gave a mean value for the specific gravity of 0.9999892, 
giving a difference in the two specific gravities of 10.2 
parts per million with an error of not over 1 part per 
million. The calculated difference in density is 10.3 parts per 
million using 1.047 as the fractionation factor. The differ- 
ence between these two values is within the experimental 
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issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


limit of error and also the probable error in the calculations 
as estimated by Urey and Greiff. 

We are at present constructing a counter-current 
scrubbing apparatus with the object of separating the 
oxygen isotopes by the method outlined by these authors. 


L. A. WEBSTER 
M. H. Want 
H. C. Urey 
Department of Chemistry, 
Columbia University, 
January 16, 1935. 


1 Urey and Greiff, J. Am. Chem. Soc. 57, in press (January, 1935). 

2 Gilfillan, J. Am. Chem. Soc. 56, 406 (1934). 

* This water was kindly prepared for us by Dr. W. G. Brown and 
A. F. Daggett, who will describe their plant in a forthcoming publi- 
cation. 


A Correction and Addition to the Discussion of the 
Ground State of H; 


Our recent paper! on the ground state of H, contains a 
minute error, detected on the eve of publication, about 
which we wish to make a more complete statement. When 
correctly computed the energies corresponding to the wave 
functions defined in that paper come out lower than the 
values there given by 0.003 e.v. for the 11-term function 
and by 0.001 e.v. for the 13-term function at the equi- 
librium distance. Furthermore, the functions are not 
accurately normalized; they should be reduced by approxi- 
mately 0.05 percent, or the corresponding correction 
applied to the result of any calculation based on their 
published form. The 5-term function and the corresponding 
energy value were not affected by the error. In principle, 
it would be possible to improve the functions (aside from 
the question of normalization) by revising the ratios of the 
coefficients, but we have not undertaken this, as it is 
clear that the resulting improvement would be quite 
negligible. 

The largest binding energy computed was —4.698 e.v. 
Correcting this by 0.006+0.003 e.v. for the effect of 
nuclear motion,? and making an allowance of —0.03+0.01 
e.v.3 for the improvement available from the addition of 
other terms, we get —4.722+0.013 e.v. as the binding 
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energy of H; at the potential minimum. Correcting for the 
zero point vibration we then obtain for the dissociation 
energy of HzDy,=4.454+0.013 e.v., as the probable result 
of a complete theoretical treatment of the problem. 

We should like to mention also an unsuccessful attempt 
to improve our method by the introduction of the Wang 
function as a term in the series. For larger separations of 
the nuclei this would certainly be helpful, but for those 
under consideration this was not the case. The computa- 
tions became somewhat more complicated through loss of 
uniformity in character, and no results were obtained on 
introduction of this function which could not be duplicated 
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by the addition in its place of one or two terms of the 
standard form. 
HusBert M. JAMEs 
ALBERT SPRAGUE COOLIDGE 
Research Laboratory of Physics, 
Mallinckrodt Chemical Laboratory, 
Harvard University, 
January 14, 1935. 


1H. M. James and A. S. Coolidge, J. Chem. Phys. 1, 825 (1933). 

2 A first order correction for this effect entered our work through the 
comparison of the fixed nucleus problem for both molecule and atoms. 

3 This differs from our previous estimate, which was —0.03 +0.02 e.v. 
An example of the method used in estimating convergence is contained 
in a paper on the 1se2pe 22, state, now in preparation. 
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